Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,136,026 G→T 5.3% V155L (GTA→TTA) ‡ 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,0260GT 5.3% 219.1 / 2.8 76V155L (GTA→TTA) ‡flgHflagellar protein of basal‑body outer‑membrane L ring
Reads supporting (aligned to +/- strand):  ref base G (26/46);  new base T (2/2);  total (28/48)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01

CTTTGATACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGC  >  NC_000913/1135878‑1136171
                                                                                                                                                    |                                                                                                                                                 
cTTTGATACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGt                                                                                                                                                  <  2:1195341/150‑1 (MQ=255)
cTTTGATACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGt                                                                                                                                                  <  1:757641/150‑1 (MQ=255)
  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:954446/1‑150 (MQ=255)
                                                                                 cAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTgg                                                                 <  2:407280/150‑1 (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:847152/1‑150 (MQ=255)
                                                                                                                 ccTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGgcagca                                 <  2:444311/150‑1 (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         >  1:1235458/3‑150 (MQ=255)
                                                                                                                                         ccgtTGACCAGTACCTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTgg         >  2:383930/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)
                                                                                                                                                    |                                                                                                                                                 
CTTTGATACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGC  >  NC_000913/1135878‑1136171

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: