Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 749,648 A→T 5.1% A20A (GCT→GCA ybgP ← putative fimbrial chaperone YbgP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913749,6480AT 5.1% 507.0 / 11.1 156A20A (GCT→GCAybgPputative fimbrial chaperone YbgP
Reads supporting (aligned to +/- strand):  ref base A (70/78);  new base T (6/2);  total (76/80)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

Only 100 of 160 total aligned reads displayed.

CAGAAGAGCATCGCTTTTCTCCCCTTTTTCATTTTCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCTGGCAGATGAATAATGCAGCTTTGCTC  >  NC_000913/749501‑749796
                                                                                                                                                   |                                                                                                                                                    
 agaagaGCATCGCTTTTCTCCCCTTTTTCATTTTCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCgcagca                                                                                                                                                   >  2:1934827/1‑150 (MQ=255)
     gagCATCGCTTTTCTCCCCTTTTTCATTTTCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGc                                                                                                                                               <  1:2019261/150‑1 (MQ=255)
     gagCATCGCTTTTCTCCCCTTTTTCATTTTCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGc                                                                                                                                               <  2:164896/150‑1 (MQ=255)
       gCATCGCTTTTCTCCCCTTTTTCATTTTCTATCCAGGAATAAGCAAGATATGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGGTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCaa                                                                                                                                             <  1:2381060/150‑1 (MQ=255)
       gCATCGCTTTTCTCCCCTTTTTCATTTTCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCaa                                                                                                                                             >  2:929926/1‑150 (MQ=255)
        cATCGCTTTTCTCCCCTTTTTCATTTTCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAAc                                                                                                                                            >  1:988953/1‑150 (MQ=255)
        cATCGCTTTTCTCCCCTTTTTCATTTTCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAAc                                                                                                                                            <  1:845958/150‑1 (MQ=255)
        cATCGCTTTTCTCCCCTTTTTCATTTTCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAAc                                                                                                                                            >  2:1831727/1‑150 (MQ=255)
          tCGCTTTTCTCCCCTTTTTCATTTTCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCAATGCTGCGGTGCAAccc                                                                                                                                          <  2:934973/150‑1 (MQ=255)
          tCGCTTTTCTCCCCTTTTTCATTTTCTACCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAAccc                                                                                                                                          <  1:85256/150‑1 (MQ=255)
             cTTTTCTCCCCTTTTTCATTTTCTATCCTGGACTAAGCAAGATAAGGCAGTTTGGCACTATGGTTTTCGATCCGCAGGCTGGTGGCTTTGTCAATGGCGTTAAACACAATGCGTGTGGTATCGGTTTGCAATGCTGCGGTGCAAcccgat                                                                                                                                       <  1:2073819/150‑4 (MQ=255)
              ttttCTCCCCTTTTTCATTTTCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTTGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGAAGCATTGCAACCCAAGc                                                                                                                                      <  2:799621/150‑1 (MQ=255)
               tttCTCCCCTTTTTCATTTTCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCt                                                                                                                                     <  1:1538708/150‑1 (MQ=255)
               tttCTCCCCTTTTTCATTTTCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCt                                                                                                                                     >  1:2261294/1‑150 (MQ=255)
               tttCTCCCCTTTTTCATTTTCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCt                                                                                                                                     <  1:1502366/150‑1 (MQ=255)
                ttCTCCCCTTTTTCATTTTCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTg                                                                                                                                    <  2:1491228/150‑1 (MQ=255)
                   tcCCCTTTTTCATTTTCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATg                                                                                                                                 >  2:642127/1‑150 (MQ=255)
                   tcCCCTTTTTCATTTTCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATg                                                                                                                                 >  1:578083/1‑150 (MQ=255)
                     cccTTTTTCATTTTCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTCCAACCCAAGCTGATGGt                                                                                                                               >  2:1338200/1‑150 (MQ=255)
                        tttttCATTTTCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTaa                                                                                                                            >  2:129017/1‑150 (MQ=255)
                          tttCATTTTCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTaaca                                                                                                                          <  2:1402459/150‑1 (MQ=255)
                            tCATTTTCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTaacaac                                                                                                                        <  2:480394/150‑1 (MQ=255)
                               ttttCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAAGAtt                                                                                                                     >  1:319995/1‑150 (MQ=255)
                               ttttCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACAtt                                                                                                                     <  2:610005/150‑1 (MQ=255)
                                tttCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTa                                                                                                                    <  1:1028020/150‑1 (MQ=255)
                                  tCTATCCAGTAATAAGCAAGATAAGTCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAg                                                                                                                  <  1:2439679/150‑1 (MQ=255)
                                  tCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAg                                                                                                                  <  2:133481/150‑1 (MQ=255)
                                   cTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAgg                                                                                                                 <  2:1015639/150‑1 (MQ=255)
                                    tATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAggg                                                                                                                <  2:418416/150‑1 (MQ=255)
                                    tATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAggg                                                                                                                <  2:565320/150‑1 (MQ=255)
                                     aTCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGt                                                                                                               >  2:1228213/1‑150 (MQ=255)
                                      tCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTGAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTa                                                                                                              >  1:1955885/1‑150 (MQ=255)
                                          ggAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATcc                                                                                                          <  2:1418158/150‑1 (MQ=255)
                                            aaTAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCtt                                                                                                        >  2:1769351/1‑150 (MQ=255)
                                               aaGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTtaa                                                                                                     >  2:2356643/1‑150 (MQ=255)
                                                aGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTtaat                                                                                                    >  1:336397/1‑150 (MQ=255)
                                                aGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTtaat                                                                                                    <  2:988953/150‑1 (MQ=255)
                                                 gCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTtaata                                                                                                   <  1:1116923/150‑1 (MQ=255)
                                                 gCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTtaata                                                                                                   >  1:94425/1‑150 (MQ=255)
                                                     gATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATg                                                                                               <  1:1684158/150‑1 (MQ=255)
                                                         aGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAg                                                                                                                  <  2:2402624/127‑1 (MQ=255)
                                                          ggCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAAcc                                                                                                                                           >  1:46689/1‑101 (MQ=255)
                                                          ggCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAAcc                                                                                                                                           <  2:46689/101‑1 (MQ=255)
                                                           gCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATtc                                                                                         <  1:2198328/150‑1 (MQ=255)
                                                           gCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATtc                                                                                         <  2:794807/150‑1 (MQ=255)
                                                                ttGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCttt                                                                                    >  1:1876797/1‑150 (MQ=255)
                                                                   tCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTcc                                                                                 >  2:388655/1‑150 (MQ=255)
                                                                   tCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTcc                                                                                 >  1:746305/1‑150 (MQ=255)
                                                                     aCTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCtt                                                                               <  2:2251877/150‑1 (MQ=255)
                                                                     aCTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCtt                                                                               >  2:2170769/1‑150 (MQ=255)
                                                                          ggTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCttaatta                                                                          <  1:1228213/150‑1 (MQ=255)
                                                                            ttttCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTtaattaat                                                                        <  1:929926/150‑1 (MQ=255)
                                                                             tttCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATg                                                                       <  2:1375343/150‑1 (MQ=255)
                                                                             tttCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATg                                                                       <  2:1714024/150‑1 (MQ=255)
                                                                             tttCGATACGCAGGCTCGTGCCTTTGTCATTGGCGTAAAACACAATGCGAGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAAAAACAATAAGGGTAATCCTTTAATAAATGTCATTCACTTTACCTTAATTAATg                                                                       >  2:2414744/1‑150 (MQ=255)
                                                                               tCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGAc                                                                     <  2:2089272/150‑1 (MQ=255)
                                                                                  aTACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGTGTTGCACCGCAGCATTGCAACCGAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAg                                                                  <  1:1630061/150‑1 (MQ=255)
                                                                                   tACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAgg                                                                 <  2:2261294/150‑1 (MQ=255)
                                                                                     cGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTa                                                               >  2:636989/1‑150 (MQ=255)
                                                                                        aGGCGGTTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTaaaa                                                            <  2:1053656/150‑1 (MQ=255)
                                                                                         ggCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAAt                                                           >  1:1001824/1‑150 (MQ=255)
                                                                                         ggCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAAt                                                           >  1:795124/1‑150 (MQ=255)
                                                                                             ggtggCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAAGTAATGACAAGGTAAAATCAgg                                                       >  2:1754179/1‑150 (MQ=255)
                                                                                                  ctatctCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCCCAGCATTGCAACCAAAGCTCATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATAAATGACAAGGTAAAATCAGGCGttt                                                  >  1:1063701/6‑150 (MQ=255)
                                                                                                   tttGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTg                                                 <  1:2400718/150‑1 (MQ=255)
                                                                                                   tttGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTg                                                 <  1:349883/150‑1 (MQ=255)
                                                                                                    ttGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTgg                                                <  1:221117/150‑1 (MQ=255)
                                                                                                     tGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGt                                               >  2:2196858/1‑150 (MQ=255)
                                                                                                        cATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCgt                                            >  1:2153350/1‑150 (MQ=255)
                                                                                                        cATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGACgt                                            >  2:328547/1‑150 (MQ=255)
                                                                                                               gTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTcc                                     <  1:1934827/150‑1 (MQ=255)
                                                                                                                  aaaCACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGa                                  <  2:2078695/150‑1 (MQ=255)
                                                                                                                    acacAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAgg                                                       >  2:1953011/1‑127 (MQ=255)
                                                                                                                    acacAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCGGACg                                >  1:293153/1‑150 (MQ=255)
                                                                                                                    acacAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACg                                <  1:642127/150‑1 (MQ=255)
                                                                                                                       caATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTc                             <  1:1831727/150‑1 (MQ=255)
                                                                                                                        aaTGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCt                            <  1:1183813/150‑1 (MQ=255)
                                                                                                                        aaTGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCt                            <  2:613263/150‑1 (MQ=255)
                                                                                                                         aTGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCTg                           <  2:2031435/150‑1 (MQ=255)
                                                                                                                         aTGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCTg                           >  2:1839552/1‑150 (MQ=255)
                                                                                                                          tGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCTgg                          <  2:2141333/150‑1 (MQ=255)
                                                                                                                          tGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCTgg                          <  1:1404443/150‑1 (MQ=255)
                                                                                                                           gCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCTGGc                         <  1:2170769/150‑1 (MQ=255)
                                                                                                                               gtgCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCTGGCAGAt                     <  1:1224719/150‑1 (MQ=255)
                                                                                                                               gtgCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCTGGCAGAt                     <  2:53416/150‑1 (MQ=255)
                                                                                                                               gtgCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCTGGCAGAt                     >  1:1401447/1‑150 (MQ=255)
                                                                                                                                 gCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCTGGCAGATGa                   >  2:1260738/1‑150 (MQ=255)
                                                                                                                                 gCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCTGGCAGATGa                   >  1:1493672/1‑150 (MQ=255)
                                                                                                                                  cGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGCCGTTTGGTCGTGTCTTCCAGACGTTCTGGCAGATGaa                  >  2:1933576/1‑150 (MQ=255)
                                                                                                                                  agcttgggTTGCAATGCTGCGGTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCTGGCAGATGaa                  >  1:565646/6‑150 (MQ=255)
                                                                                                                                   gcttgggTTGCAATGCTGCGGTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCTGGCAGATGaat                 >  2:304747/5‑150 (MQ=255)
                                                                                                                                    ttCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCTGGCAGATGaata                <  1:2414744/149‑1 (MQ=255)
                                                                                                                                    aTCGGGGTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGGCGTGGCTTCCAGACGTTCTGGCAGATGaata                <  2:1644237/150‑1 (MQ=255)
                                                                                                                                     ttgggTTGCAATGCTGCGGTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCTGGCAGATGaataa               >  2:1215023/3‑150 (MQ=255)
                                                                                                                                     tCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCTGGCAGATGaataa               >  2:1084023/1‑150 (MQ=255)
                                                                                                                                        ggTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCTGGCAGATGAATAATGc            >  2:683301/1‑150 (MQ=255)
                                                                                                                                        ggTTGCAATGCTGCGGTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCTGGCAGATGAATAATGc            >  1:321485/1‑150 (MQ=255)
                                                                                                                                           tGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCTGGCAGATGAATAATGCAGc         >  2:2010975/1‑150 (MQ=255)
                                                                                                                                            gCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAAAATTAAGGGTAATCCTTTAATAAATGTCATTATCTTTTCCTTAATTAATCACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGAAGTTCTGGCAGATGAATAATGAAGCt        >  2:1243851/1‑150 (MQ=255)
                                                                                                                                               ccGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCTGGCAGATGAATAATGCAGCTTTg     >  2:1101357/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CAGAAGAGCATCGCTTTTCTCCCCTTTTTCATTTTCTATCCAGGAATAAGCAAGATAAGGCAGTTTGTCACTTTGGTTTTCGATACGCAGGCTGGTGGCTTTGTCATTGGCGTTAAACACAATGCGTGTGCGATCGGGTTGCACCGCAGCATTGCAACCCAAGCTGATGGTTAACAACATTAAGGGTAATCCTTTAATAAATGTCATTCTCTTTTCCTTAATTAATGACAAGGTAAAATCAGGCGTTTGGTCGTGTCTTCCAGACGTTCTGGCAGATGAATAATGCAGCTTTGCTC  >  NC_000913/749501‑749796

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