Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,017,124 C→T 5.3% L191L (CTG→TTG)  fliI → flagellum‑specific ATP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,017,1240CT 5.3% 241.6 / 2.2 75L191L (CTG→TTG) fliIflagellum‑specific ATP synthase
Reads supporting (aligned to +/- strand):  ref base C (29/42);  new base T (2/2);  total (31/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

CTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGCACGCTCAGTGGTGATTGCCGCTCCGGC  >  NC_000913/2016977‑2017269
                                                                                                                                                   |                                                                                                                                                 
cTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGtgct                                                                                                                                                  >  7:314977/1‑149 (MQ=255)
cTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGtgct                                                                                                                                                  >  1:68406/1‑149 (MQ=255)
 tGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGtgctg                                                                                                                                                 >  5:100043/1‑149 (MQ=255)
   aTTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTTTTTGCCGGGTCCGGCGTTGGTAAAAGTGTgctgct                                                                                                                                               >  6:289174/1‑149 (MQ=255)
    ttACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTCCCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGTCGGGTCCGGCGTTGGTAAAAGTGTGctgctg                                                                                                                                              <  8:337692/149‑1 (MQ=255)
       ccccGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGctgctg                                                                                                                                              >  3:149806/1‑146 (MQ=255)
       ccccGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGctgctg                                                                                                                                              <  4:149806/146‑1 (MQ=255)
        cccGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGga                                                                                                                                          <  1:38716/149‑1 (MQ=255)
         ccGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTAGGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGgat                                                                                                                                         >  2:370915/1‑149 (MQ=255)
           gCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGgatga                                                                                                                                       <  3:77337/149‑1 (MQ=255)
            ccATTTAACCCGTGGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGgatgat                                                                                                                                      <  2:65533/149‑1 (MQ=255)
            ccACTTAACCCGTTGCAACGTCCACCGATTGACCATGGGCTGGACACCGGCGTGCGCCCATTCAATTCCCTGCTTCCCGCTGGGCGGGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCGGGGgatgat                                                                                                                                      <  2:259861/149‑1 (MQ=255)
              aTTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATgg                                                                                                                                    >  7:170487/1‑149 (MQ=255)
                     ccGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTa                                                                                                                             >  5:298600/1‑149 (MQ=255)
                      cGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTac                                                                                                                            <  1:352747/149‑1 (MQ=255)
                      cGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTa                                                                                                                             >  2:183772/1‑148 (MQ=255)
                       gTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTaca                                                                                                                           <  1:360920/149‑1 (MQ=255)
                          gCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACAcc                                                                                                                         <  7:217421/148‑1 (MQ=255)
                           cAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCcg                                                                                                                       <  8:6500/149‑1 (MQ=255)
                                 acacCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGa                                                                                                                 >  7:291531/1‑149 (MQ=255)
                                  cacCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGAt                                                                                                                <  1:183772/149‑1 (MQ=255)
                                   acCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATg                                                                                                               <  7:70028/149‑1 (MQ=255)
                                      gATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCa                                                                                                            >  8:8851/1‑149 (MQ=255)
                                       aTTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCAt                                                                                                           >  3:307640/1‑149 (MQ=255)
                                           aaCATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGACATTGTc                                                                                                       >  2:370837/1‑149 (MQ=255)
                                             cttgtgCTGGACCCCGGCGTGCGCCCAATCAAGGCCCTGCTTACCGTTGGGCGGGGGCAGCGTATGGGGCGGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGt                                                                                                     <  6:111155/147‑1 (MQ=255)
                                               tgttCTTGACACCGGCGTGCGCCAAATCAATGCCCTGCCTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTTGTAAAAGTGTGCTGCTGGTGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTgg                                                                                                   <  4:116197/149‑1 (MQ=255)
                                                   gtGGAGACCGGCGTGCGCCCAATCAATGCCCTGCGTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGttt                                                                                               <  5:360037/148‑1 (MQ=255)
                                                   cTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGttt                                                                                               >  3:137541/1‑149 (MQ=255)
                                                      gACACCGGCGTGCGCCCAATCACTGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGGAAAAGGGTGCTGCTGGGGACTATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGag                                                                                            >  4:282528/1‑148 (MQ=255)
                                                      gACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGAt                                                                                            >  3:160960/1‑149 (MQ=255)
                                                      gACACCGGCGTGCGACCAATCAATGCCCAGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGGTTGAt                                                                                            >  8:222669/1‑149 (MQ=255)
                                                          ccGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGt                                                                                        <  4:307640/149‑1 (MQ=255)
                                                           cGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGtt                                                                                                                              >  4:317014/1‑110 (MQ=255)
                                                           cGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGtt                                                                                                                              <  3:317014/110‑1 (MQ=255)
                                                             gCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGaa                                                                                     <  4:2126/149‑1 (MQ=255)
                                                                tGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGt                                                                                  <  8:291531/149‑1 (MQ=255)
                                                                   gcCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGtt                                                                                                >  4:158451/1‑132 (MQ=255)
                                                                   gcCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGtt                                                                                                <  3:158451/132‑1 (MQ=255)
                                                                    cccAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGgc                                                                              <  8:314977/149‑1 (MQ=255)
                                                                    cccAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGgc                                                                              >  6:134618/1‑149 (MQ=255)
                                                                     ccAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGTTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGaa                                                                                     <  7:333365/141‑1 (MQ=255)
                                                                     ccAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGTTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGaa                                                                                     <  7:327738/141‑1 (MQ=255)
                                                                     ccAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGTTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGaa                                                                                     >  8:327738/1‑141 (MQ=255)
                                                                     ccAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGTTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGaa                                                                                     >  8:333365/1‑141 (MQ=255)
                                                                     ccAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGgcg                                                                             >  7:393242/1‑149 (MQ=255)
                                                                      caatcaatGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTc                                                                                                       >  6:244879/1‑122 (MQ=255)
                                                                      caatcaatGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTc                                                                                                       <  5:244879/122‑1 (MQ=255)
                                                                           aatGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTGCACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGt                                                                       <  7:2138/149‑1 (MQ=255)
                                                                               cccTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTaaaa                                                                   <  4:137541/149‑1 (MQ=255)
                                                                               cccTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTaaa                                                                    <  7:8851/148‑1 (MQ=255)
                                                                                 cTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTaa                                                                     >  2:359331/1‑145 (MQ=255)
                                                                                    cTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGAttt                                                              <  1:370837/149‑1 (MQ=255)
                                                                                       aCCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTa                                                            <  8:170487/148‑1 (MQ=255)
                                                                                        ccGTTGGGCGTGGGCAGGGTATGGGGGTGTTTGCCGGGTCCGGCGTTGGTAATAGTGTGCTGCTGTGGATGATGGCAGGTTACACCCGCGCCGATGTCTTTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTTAAAGtttttttt                                                          <  3:217890/149‑1 (MQ=255)
                                                                                        ccGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTAtt                                                          <  8:372956/149‑1 (MQ=255)
                                                                                         cGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTAtt                                                          >  7:100173/1‑148 (MQ=255)
                                                                                          gTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTAtt                                                          >  1:289056/1‑147 (MQ=255)
                                                                                             gggCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTgaga                                                      <  7:130788/148‑1 (MQ=255)
                                                                                              ggCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGaa                                                     <  5:304143/148‑1 (MQ=255)
                                                                                                cGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTTGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACAt                                                  <  6:100043/149‑1 (MQ=255)
                                                                                                      cAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATTGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGtt                                                                                                <  5:51382/97‑1 (MQ=255)
                                                                                                      cAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATTGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGtt                                                                                                >  6:51382/1‑97 (MQ=255)
                                                                                                        gCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTg                                          >  3:176041/1‑149 (MQ=255)
                                                                                                         cGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGc                                         <  6:298600/149‑1 (MQ=255)
                                                                                                           tATGGGGCTGTTTGCTGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGcc                                        >  8:47941/1‑148 (MQ=255)
                                                                                                                     tttGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGt                               <  1:370915/147‑1 (MQ=255)
                                                                                                                     tttGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGc                             <  7:222669/149‑1 (MQ=255)
                                                                                                                           gggTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGc                             <  2:363988/143‑1 (MQ=255)
                                                                                                                           gggTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGc                             >  1:363988/1‑143 (MQ=255)
                                                                                                                            ggTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGCACGCTCa                      <  2:289056/149‑1 (MQ=255)
                                                                                                                                cGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGCACGCTCAgtgg                  >  5:222251/1‑149 (MQ=255)
                                                                                                                                cGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGCACGCTCAgtg                   <  7:47941/148‑1 (MQ=255)
                                                                                                                                cGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGCACGCTCAgtg                   <  8:110510/148‑1 (MQ=255)
                                                                                                                                   cGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGCACGCTCAGTGGTGa               <  2:68406/149‑1 (MQ=255)
                                                                                                                                     ttGGAAAAAGTGGGCTGCTGGGGACGATGGCACGTGACACCCGCGCCGATGTCATTGTCGTGGGTTTGCTTGGTGCACGGGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGCACGCTCAGTGGTGAtt             <  3:282528/149‑1 (MQ=255)
                                                                                                                                        gTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGCACGCTCAGTGGTGATTGcc          >  6:384574/1‑149 (MQ=255)
                                                                                                                                        gTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGCACGCTCAGTGGTGATTGc           <  7:359933/148‑1 (MQ=255)
                                                                                                                                         tAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGCACGCTCAGTGGTGATTGCCg         >  2:7206/1‑149 (MQ=255)
                                                                                                                                          aaaaGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTg                              >  2:337894/1‑127 (MQ=255)
                                                                                                                                          aaaaGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTg                              <  1:337894/127‑1 (MQ=255)
                                                                                                                                                gtgCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGCACGCTCAGTGGTGATTGCCGCTCcggc  >  7:296646/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGCACGCTCAGTGGTGATTGCCGCTCCGGC  >  NC_000913/2016977‑2017269

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: