Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,136,028 A→C 5.3% V155V (GTA→GTC) ‡ flgH → flagellar protein of basal‑body outer‑membrane L ring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,136,0280AC 5.3% 220.9 / 3.0 75V155V (GTA→GTC) ‡flgHflagellar protein of basal‑body outer‑membrane L ring
Reads supporting (aligned to +/- strand):  ref base A (27/44);  new base C (2/2);  total (29/46)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.91e-01

TTGATACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCAT  >  NC_000913/1135880‑1136177
                                                                                                                                                    |                                                                                                                                                     
ttGATACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTAc                                                                                                                                                      >  1:444311/1‑150 (MQ=255)
ttGATACTCTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTAGAAGCCTCCGGTGGTAACGCGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTaa                                                                                                                                                      >  2:1281335/1‑149 (MQ=255)
 tGATACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACt                                                                                                                                                     <  2:433023/150‑1 (MQ=255)
   aTACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTgg                                                                                                                                                   <  2:938461/150‑1 (MQ=255)
       tgtgCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCaa                                                                                                                                               >  2:92471/1‑150 (MQ=255)
        gtgCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGTACCTGGTCAAc                                                                                                                                              <  1:478828/150‑1 (MQ=255)
         tgCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGTACCTGGTCAAcc                                                                                                                                             <  2:982288/150‑2 (MQ=255)
          gCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCGCGTGCCGATGTCGACGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACgg                                                                                                                                            <  1:927770/150‑1 (MQ=255)
            cgcgCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGGGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCa                                                                                                                                          <  1:794145/150‑1 (MQ=255)
             gcgcTATTTGCAGGGGCTGGTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCaa                                                                                                                                         <  1:784506/150‑1 (MQ=255)
                cTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAAccg                                                                                                                                      >  2:345062/1‑149 (MQ=255)
                  aTTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGc                                                                                                                                    <  1:593848/150‑1 (MQ=255)
                   tttGCCGGGGCGGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCa                                                                                                                                   <  2:491394/150‑1 (MQ=255)
                    ttGCAGGGGCTGTTTGGTAACGCTCGGGCCGATGTCGAAGCCTCCGGTGGTAACACGCTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCAt                                                                                                                                  <  1:94688/150‑1 (MQ=255)
                     tGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCAtg                                                                                                                                 <  2:780694/150‑1 (MQ=255)
                         ggggCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGCCGGTTGACCAGGTACTGGTCAACGGCAACCTGCATgtggt                                                                                                                             >  1:866339/1‑150 (MQ=255)
                         ggggCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTGCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGGTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATgtggt                                                                                                                             <  1:790724/150‑1 (MQ=255)
                         ggggCTGTTTGGTAACGCGCGTGCCGATGTCGCAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATgtggt                                                                                                                             <  2:45803/150‑1 (MQ=255)
                          gggCTGGTTGGTAACGCTCGGGCCGAGGTCGAAGCCTCCGGTGGTACCACGGTCACCGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATgtggtg                                                                                                                            <  1:231281/150‑1 (MQ=255)
                           ggCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGtggtgg                                                                                                                           >  1:1313256/1‑150 (MQ=255)
                           gacTGTTTGGTAACGCGCGTGCCGATGTCGACGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGtggtgg                                                                                                                           <  1:800690/148‑1 (MQ=255)
                            gCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTggtggg                                                                                                                          <  2:437483/150‑1 (MQ=255)
                                  tGGTAACGCTCGTGCCGATGACGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCACCCTGCATGTGGTGGGTGaaaa                                                                                                                    >  1:769114/1‑150 (MQ=255)
                                   ggTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTACCACGTTCAACGGAAAGGGCGGGGCCAAAGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGCCCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGaaaaa                                                                                                                   >  1:1061066/1‑150 (MQ=255)
                                   ggTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCACCCTGCATGTGGTGGGTGaaaaa                                                                                                                   >  2:1130824/1‑150 (MQ=255)
                                    gTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAAc                                                                                                                  <  1:1021578/150‑1 (MQ=255)
                                    gTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGGGAAAAAc                                                                                                                  <  1:678877/150‑1 (MQ=255)
                                     tAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACa                                                                                                                 <  2:334728/150‑1 (MQ=255)
                                      aaCGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGCCGGTTGACCAGGTACTGGTCAACGGCAACCTGCAGGTGGTGGGTGAAAAACAg                                                                                                                >  2:545026/1‑150 (MQ=255)
                                         gagcGTGCCGATGTCGAAGCGTCCGGTGGTTACACTTTCAAGGGTAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGAtt                                                                                                             <  2:977379/147‑1 (MQ=255)
                                           tCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGc                                                                                                           <  2:462598/150‑1 (MQ=255)
                                            cGTGCCGATGACGAAGCCTCCGGAGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCCGCAATACCTTTAGCGGCACGTTGACGGTGCCGGTTGACCACGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGcc                                                                                                          >  1:181805/1‑150 (MQ=255)
                                             gggCCGATGGCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAAGGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCa                                                                                                         <  1:127697/148‑1 (MQ=255)
                                               gCCGATGTCGAAGCCTCCGGTGGTAACACGTTCCACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCAtt                                                                                                       <  1:976620/150‑1 (MQ=255)
                                                    tGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCa                                                                                                  >  1:439336/1‑150 (MQ=255)
                                                     gTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAAGCAg                                                                                                 >  1:878411/1‑150 (MQ=255)
                                                      tggAAGCCGCCGGTGGTAACACGTTCTACGGAAAGGTCGGTGAGAATGCCATCATTACCTTAAGCGGCACGTTGACGGGGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAgg                                                                                                <  2:208704/148‑1 (MQ=255)
                                                      tCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCAAGCAATACCTTTAGCGGCCCGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAACCAGATTGCCATTAATCAgg                                                                                                >  1:347164/1‑150 (MQ=255)
                                                       cGAAGCCTCCGGTGGAAACACGTTCAACGGAAAGGGCGGGGCAAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAggg                                                                                               <  1:828132/150‑1 (MQ=255)
                                                        gAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGt                                                                                              <  1:306502/150‑1 (MQ=255)
                                                          aGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTAc                                                                                            <  2:863479/150‑1 (MQ=255)
                                                            ccTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCg                                                                                          >  1:68431/1‑150 (MQ=255)
                                                                 ggtggtAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAAttt                                                                                     <  1:960719/150‑1 (MQ=255)
                                                                      tAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCg                                                                                >  2:407796/1‑150 (MQ=255)
                                                                       aaCACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGc                                                                               <  1:858851/150‑1 (MQ=255)
                                                                             ttCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTctct                                                                         <  2:331383/150‑1 (MQ=255)
                                                                               cAACGGAAAGGGCGGGGCCAATGCCATCAATACCTTTAGCGGCACGTTGACGGGTACGGATGAACATGTACTTGTCAACGGCAACCTGCATGTTGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTgg                                                                       <  1:328943/150‑1 (MQ=255)
                                                                               cAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTgg                                                                       <  2:407280/150‑1 (MQ=255)
                                                                               cAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTgg                                                                       >  2:954446/1‑150 (MQ=255)
                                                                                    gAAAGGGCGGGGCCAATGCCAGCAATACCTGTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGGCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTgg                                                                  >  1:503581/1‑150 (MQ=255)
                                                                                         ggCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTCCCGAATTTATTCGCTTCTCTGGCGTGGTGAAt                                                             >  1:963469/1‑150 (MQ=255)
                                                                                              ggCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGc                                                       <  1:608949/150‑1 (MQ=255)
                                                                                                ccAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACGGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCa                                                      <  1:603193/150‑1 (MQ=255)
                                                                                                  aaTGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACt                                                    <  1:177429/150‑1 (MQ=255)
                                                                                                      ccAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCa                                                >  1:442562/1‑150 (MQ=255)
                                                                                                         gCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATATATTCGCTTCTCAGGCGTGGTTAATCCACGCACTATCAGCg                                             >  2:300984/1‑150 (MQ=255)
                                                                                                           aaTACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGgc                                           <  2:612426/150‑1 (MQ=255)
                                                                                                               ccTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGgcagca                                       >  1:1215858/1‑150 (MQ=255)
                                                                                                               ccTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGgcagca                                       <  2:444311/150‑1 (MQ=255)
                                                                                                               ccTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGgcagca                                       >  2:847152/1‑150 (MQ=255)
                                                                                                               ccTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGACAACGGCAACCTGAATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGTCGTGGTTAATCCACGCACTATCAGCGgcagca                                       >  1:1044808/1‑150 (MQ=255)
                                                                                                                  ttAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAAta                                    <  1:539437/150‑1 (MQ=255)
                                                                                                                     gCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTCGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAAtaccg                                 <  2:717955/150‑1 (MQ=255)
                                                                                                                      cGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAAtaccgt                                <  2:318199/150‑1 (MQ=255)
                                                                                                                        gCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGGGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAAtaccgtac                              <  2:1143538/150‑1 (MQ=255)
                                                                                                                             ttGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCt                         <  1:945873/150‑1 (MQ=255)
                                                                                                                                     gacggtTGACCAGGTACTTGTCAACGGCAACCTGCATGTTGTGGGTGAAAAACAGATTGTCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGt                 <  2:239679/150‑1 (MQ=255)
                                                                                                                                      acggtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTg                <  1:287904/150‑1 (MQ=255)
                                                                                                                                       ccgtTGACCAGTACCTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTgg               >  2:383930/3‑150 (MQ=255)
                                                                                                                                       ccgtTGACCAGTACCTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTgg               >  1:1235458/3‑150 (MQ=255)
                                                                                                                                        ggtTGACCATGTACTGGTCAACGGCAACCTTCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTTGc              <  2:981159/150‑1 (MQ=255)
                                                                                                                                         gtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCg             <  1:92471/150‑1 (MQ=255)
                                                                                                                                              ccAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACAAAGGTGGCTGATgc        >  1:640607/1‑150 (MQ=255)
                                                                                                                                              ccAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCAACGCACGATCAGCGGCAGCAATACCGTACCGTCGACTCAGGTGGCGGATgc        >  1:922123/1‑150 (MQ=255)
                                                                                                                                                   tACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCa   >  1:429436/1‑150 (MQ=255)
                                                                                                                                                   tACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCa   >  1:392609/1‑150 (MQ=255)
                                                                                                                                                    aCTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCAt  <  2:181805/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTGATACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCGCGCAT  >  NC_000913/1135880‑1136177

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