Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,965,176 A→T 6.6% V339E (GTG→GAG)  flhB ← flagellin export apparatus, substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,965,1760AT 6.6% 174.0 / 3.7 61V339E (GTG→GAG) flhBflagellin export apparatus, substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base A (31/25);  new base T (2/2);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

GCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAGCGCGACCAGCCCTGCACCTT  >  NC_000913/1965037‑1965320
                                                                                                                                           |                                                                                                                                                
gCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCCCCAGCCAGACGCCAGCGGTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCgcggcg                                                                                                                                          >  6:31150/1‑148 (MQ=255)
    gATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGACGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAg                                                                                                                                     <  3:146120/149‑1 (MQ=255)
            aTGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCg                                                                                                                             <  1:357162/149‑1 (MQ=255)
            aTGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCg                                                                                                                             >  2:357162/1‑149 (MQ=255)
               ggTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGa                                                                                                                          <  8:183798/149‑1 (MQ=255)
                   ggTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATttg                                                                                                                      <  3:250626/149‑1 (MQ=255)
                    gTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATttgt                                                                                                                     >  1:52254/1‑149 (MQ=255)
                    gTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTACAGTTGCCAGACCCAGGCCCCCACTTCCGCCACCGCGGCGTACCGTTGACCCGGGATttgt                                                                                                                     >  6:273560/1‑149 (MQ=255)
                     tttCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATttgtt                                                                                                                    <  2:300386/149‑1 (MQ=255)
                        ctcGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGa                                                                                                                  <  2:51475/148‑1 (MQ=255)
                           gTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCAAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCAGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACaaa                                                                                                              >  6:354661/1‑149 (MQ=255)
                           gTTAATAAAATCCAGCGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGAcc                                                                                                                <  4:317183/147‑1 (MQ=255)
                              aataaAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCTCCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAAtct                                                                                                           >  4:91866/1‑149 (MQ=255)
                               ataaAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGGCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAAtctc                                                                                                          <  2:379799/149‑1 (MQ=255)
                                taaAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTcc                                                                                                         >  2:178534/1‑149 (MQ=255)
                                       cAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGGCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTc                                                                                                  <  5:243462/149‑1 (MQ=255)
                                        aGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTc                                                                                                  >  6:295507/1‑148 (MQ=255)
                                            cTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATAc                                                                                             <  2:5532/149‑1 (MQ=255)
                                             ttCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATa                                                                                              >  5:54703/1‑147 (MQ=255)
                                              tCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGATGCCAGACCCAGGCCAGCACTTCCGCCCCCGCGCCGTACAGTTGACCAGGGATTTGTTGACCAATCTCCGCAGGTGGATACAg                                                                                           >  6:316062/1‑149 (MQ=255)
                                                 ggCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTcc                                                                                                         <  6:330555/132‑1 (MQ=255)
                                                 ggCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTcc                                                                                                         >  5:330555/1‑132 (MQ=255)
                                                   cACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTc                                                                                      <  5:295507/149‑1 (MQ=255)
                                                   cACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTc                                                                                      <  7:311283/149‑1 (MQ=255)
                                                      cGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTcgcg                                                                                   >  1:170768/1‑149 (MQ=255)
                                                        gAAGATGAGTAGGTTGTACAGGGCGCTGTCAACCAGCCAGACGACAGCGTTTCAGTTGCCAGACACAGGCCAGCACTTCCGCCACCGCGGGGTAAAGTTGACCCGGGATTAGTTGACCAATCTCCGCATGTCGATACAGCGCTCgcgc                                                                                  >  8:180402/1‑148 (MQ=255)
                                                           gATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCCGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGGTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAgc                                                                              <  6:268799/149‑1 (MQ=255)
                                                                gTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAgcggcgg                                                                         >  2:390155/1‑149 (MQ=255)
                                                                   ggTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAgcggcggcg                                                                       <  5:80185/148‑1 (MQ=255)
                                                                    gTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCt                                                                     >  1:104789/1‑149 (MQ=255)
                                                                     ttGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCtt                                                                    >  2:80023/1‑149 (MQ=255)
                                                                       gTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGGTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTc                                                                   <  6:286995/148‑1 (MQ=255)
                                                                          cAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAg                                                                <  6:54703/148‑1 (MQ=255)
                                                                           aGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGCCCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCg                                                              >  7:126182/1‑149 (MQ=255)
                                                                             ggCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTc                                                            >  4:220347/1‑149 (MQ=255)
                                                                             ggCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTc                                                            >  6:48694/1‑149 (MQ=255)
                                                                                gcTGTCCACCAGCCAGACGCCAGCGGTTCAGGTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCggg                                                         <  2:357586/149‑1 (MQ=255)
                                                                                       accaGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTAt                                                  >  6:64444/1‑149 (MQ=255)
                                                                                        ccagccagACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTAtt                                                 <  5:64444/149‑1 (MQ=255)
                                                                                           gccagACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTAtttt                                               >  5:240509/1‑148 (MQ=255)
                                                                                                  gCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGcc                                        >  3:94875/1‑148 (MQ=255)
                                                                                                    cAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGc                                                                                                       >  2:63425/1‑83 (MQ=255)
                                                                                                    cAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGc                                                                                                       <  1:63425/83‑1 (MQ=255)
                                                                                                         tttCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCAc                                <  4:105450/149‑1 (MQ=255)
                                                                                                           tCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACg                               >  4:294051/1‑148 (MQ=255)
                                                                                                               ttGCCAGACCCAGGCCAGCACTTCCGCCTCCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAAt                            <  7:36693/147‑1 (MQ=255)
                                                                                                               ttGCCAGACCCAGGCCAGCACTTCCGCCTCCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAAt                            <  7:42669/147‑1 (MQ=255)
                                                                                                               ttGCCAGACCCAGGCCAGCACTTCCGCCTCCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAAt                            >  8:36693/1‑147 (MQ=255)
                                                                                                               ttGCCAGACCCAGGCCAGCACTTCCGCCTCCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAAt                            >  8:42669/1‑147 (MQ=255)
                                                                                                                   cAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAg                      >  3:338115/1‑149 (MQ=255)
                                                                                                                        ccAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAGCGCGa                 >  5:252895/1‑149 (MQ=255)
                                                                                                                        ccAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAGCGCGa                 >  3:352463/1‑149 (MQ=255)
                                                                                                                          aGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGCCCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAGCGCGAc                >  1:349951/1‑148 (MQ=255)
                                                                                                                           ggCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAGCGCGACCa              <  5:48694/149‑1 (MQ=255)
                                                                                                                              cAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAGCGCGACCAGcc           >  2:144615/1‑149 (MQ=255)
                                                                                                                                  aCTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTAtt                                                 <  4:210923/107‑1 (MQ=255)
                                                                                                                                  aCTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTAtt                                                 >  3:210923/1‑107 (MQ=255)
                                                                                                                                  aCTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAGCGCGACCAGCCCTGc       <  1:144615/149‑1 (MQ=255)
                                                                                                                                   cTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAGCGCGACCAGCCCTGCa      >  5:101806/1‑149 (MQ=255)
                                                                                                                                   cTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAGCGCGACCAGCCCTGCa      <  2:170768/149‑1 (MQ=255)
                                                                                                                                        gCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAGCGCGACCAGCCCTGCACCtt  >  6:352851/1‑148 (MQ=255)
                                                                                                                                           |                                                                                                                                                
GCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAGCGCGACCAGCCCTGCACCTT  >  NC_000913/1965037‑1965320

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑