Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,101,551 G→T 7.7% P45Q (CCG→CAG)  csgG ← curli production assembly/transport outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,101,5510GT 7.7% 328.3 / 15.8 117P45Q (CCG→CAG) csgGcurli production assembly/transport outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (59/48);  new base T (5/4);  total (65/52)
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.72e-01

Only 100 of 122 total aligned reads displayed.

GTATAAACCAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATG  >  NC_000913/1101402‑1101700
                                                                                                                                                     |                                                                                                                                                     
gTATAAACCAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCg                                                                                                                                                       >  1:2041460/1‑150 (MQ=255)
gTATAAACCAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCg                                                                                                                                                       >  2:186564/1‑150 (MQ=255)
 tataAACCCGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCACGGTTTCGTCCTGAATTATGTATACAGAAACAAAGATTTTACCCGTCgg                                                                                                                                                      >  1:1011259/1‑150 (MQ=255)
 tataAACCAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCgg                                                                                                                                                      >  2:949853/1‑150 (MQ=255)
    aaaCCAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGgcgc                                                                                                                                                   >  1:1595185/1‑150 (MQ=255)
      aCCAGCGAGAATCTTTCAGTGCCGGGACCAGCATTGCCGTGGCGCGTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGGAGGGTTTAACTTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTg                                                                                                                                                 <  1:161199/150‑1 (MQ=255)
        cAGCGAGAAGCTTTCAGAGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGc                                                                                                                                               <  2:523565/150‑1 (MQ=255)
         aGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCa                                                                                                                                              >  1:983113/1‑150 (MQ=255)
          gCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAg                                                                                                                                             >  2:1589047/1‑150 (MQ=255)
          gCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCACGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACACGTCGGCGCTGGCAg                                                                                                                                             >  1:1729005/1‑150 (MQ=255)
             agaATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATg                                                                                                                                          >  2:1968770/1‑150 (MQ=255)
             agaATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATg                                                                                                                                          <  2:1417309/150‑1 (MQ=255)
               aaTCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATggg                                                                                                                                        <  2:958233/150‑1 (MQ=255)
                   tttCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGGCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCaa                                                                                                                                    <  2:784672/150‑1 (MQ=255)
                     tCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAAt                                                                                                                                  <  2:1378656/150‑1 (MQ=255)
                     tCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAAt                                                                                                                                  <  1:508065/150‑1 (MQ=255)
                         tGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCCGCGCTGGCAGATGGGTCAAAGCttt                                                                                                                              >  1:1044762/1‑150 (MQ=255)
                          gCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTg                                                                                                                             <  2:1191621/150‑1 (MQ=255)
                          gCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTg                                                                                                                             <  1:438878/150‑1 (MQ=255)
                           ccctGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGGAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTGt                                                                                                                            <  1:1292387/147‑1 (MQ=25)
                           ccGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTGt                                                                                                                            <  1:1402739/150‑1 (MQ=37)
                             gTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTGttt                                                                                                                          <  2:2033213/150‑3 (MQ=37)
                             gTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCAGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTCCCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAg                                                                                                                          >  2:1690248/1‑150 (MQ=255)
                             gTGACCAGCATTGCCGTGGCGCTTTGCGGAAAAGCAGTGGACAAGTTACTTGCCGGGTAGGGTTGAAATTGCCCTGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTAAAATCTTTGTAg                                                                                                                          <  2:1594093/150‑1 (MQ=255)
                               gACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAACCAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTGtttcg                                                                                                                        <  1:211758/150‑5 (MQ=25)
                               gACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGct                                                                                                                        >  2:1544034/1‑150 (MQ=255)
                                  cAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTg                                                                                                                     <  1:1817278/150‑1 (MQ=255)
                                  cAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTg                                                                                                                     <  2:1095288/150‑1 (MQ=255)
                                         gCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGa                                                                                                              >  1:441307/1‑150 (MQ=255)
                                          ccGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAg                                                                                                             >  2:1707392/1‑150 (MQ=255)
                                            gTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCTCGAGGc                                                                                                           >  1:1348093/1‑150 (MQ=255)
                                                 gcTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCAAGAGGCATTaa                                                                                                      >  2:1124613/1‑150 (MQ=255)
                                                  cTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGTGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAAt                                                                                                     <  1:1708981/150‑1 (MQ=255)
                                                    ttGCGGAACAGCAGTGGAGAAGTTACATGCCGGGTAGGGTTTAAATTGCCCGGAATCGTCCTGAATGTTGTATACCGAAACAAAGATTATACACGACGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGt                                                                                                   >  1:1127304/1‑150 (MQ=255)
                                                     tGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTc                                                                                                  >  2:1500816/1‑150 (MQ=255)
                                                      gCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCg                                                                                                 >  1:1871699/1‑150 (MQ=255)
                                                       cGGAACAGCAGTGGAGAAGTTACTTGCCGGGGAGGGTTTAAATTGCCCGGGTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCgg                                                                                                <  2:184693/150‑1 (MQ=255)
                                                         gAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGTTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTc                                                                                              <  2:1620974/150‑1 (MQ=255)
                                                         gAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTc                                                                                              >  1:938030/1‑150 (MQ=255)
                                                         gAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTc                                                                                              <  2:618142/150‑1 (MQ=255)
                                                            cagcagTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTgg                                                                                           >  1:642843/1‑150 (MQ=255)
                                                            cagcagTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTgg                                                                                           >  1:1825547/1‑150 (MQ=255)
                                                             agcagTGGAGAAGTTACTTGCCGGGTTGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGTTTTTACCCGTCGGCGCTGTCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTggc                                                                                          <  1:943151/150‑1 (MQ=255)
                                                              gcagTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTggcg                                                                                         >  1:1080154/1‑150 (MQ=255)
                                                                agTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTggctga                                                                                       >  2:386172/1‑147 (MQ=255)
                                                                    gagaAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCTTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGcttct                                                                                   <  2:1697113/150‑1 (MQ=255)
                                                                    gagaAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGcttct                                                                                   <  1:52760/150‑1 (MQ=255)
                                                                    gagaAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACACGTCGGCGCTTGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGcttct                                                                                   >  1:1243451/1‑150 (MQ=255)
                                                                     agaAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTCTATACCGAAACAAAGATTTTACACGTCGGCGCTGGCAGATGGGTCAAATATTTGTAGCTCTGAGCACGCGGCATTAATGTCGGTCTGGCGGcttctt                                                                                  >  1:1847342/1‑150 (MQ=255)
                                                                      gaAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCttt                                                                                 >  1:445475/1‑150 (MQ=255)
                                                                        aGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAg                                                                               >  1:1958343/1‑150 (MQ=255)
                                                                        aGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAg                                                                               >  1:937872/1‑150 (MQ=255)
                                                                         gTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAgg                                                                              <  1:1994034/150‑1 (MQ=255)
                                                                         gTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAgg                                                                              >  1:466280/1‑150 (MQ=255)
                                                                           tACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTTTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAggcg                                                                            >  1:1801953/1‑150 (MQ=255)
                                                                           tACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAggcg                                                                            >  2:926845/1‑150 (MQ=255)
                                                                           tACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAggcg                                                                            >  2:1580397/1‑150 (MQ=255)
                                                                           tACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCAGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAggcg                                                                            >  2:1113475/1‑150 (MQ=255)
                                                                             cTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAggcggg                                                                          >  1:1989918/1‑150 (MQ=255)
                                                                             cTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAggcggg                                                                          >  2:1998209/1‑150 (MQ=255)
                                                                              ttGCCGGGGAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGAATTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAggcgggg                                                                         <  1:453229/150‑1 (MQ=255)
                                                                               tGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAggcggggc                                                                        >  1:432180/1‑150 (MQ=255)
                                                                               tGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAggcggggc                                                                        <  1:589239/150‑1 (MQ=255)
                                                                               tGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAggcggggc                                                                        <  1:680953/150‑1 (MQ=255)
                                                                                 ccGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAggcggggcgg                                                                      >  2:594922/1‑150 (MQ=255)
                                                                                  cGGGTAGGGTTTAAATTGCCCTGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCTGCTTCTTTAGGCGGGGCGGt                                                                     <  2:1709909/150‑1 (MQ=255)
                                                                                  cGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGt                                                                     >  2:1130069/1‑150 (MQ=255)
                                                                                   gggTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGtt                                                                    <  2:1331293/150‑1 (MQ=255)
                                                                                   gggTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGtt                                                                    >  1:1850029/1‑150 (MQ=255)
                                                                                      tAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGCGGCGGTTAAg                                                                 >  1:1800914/1‑150 (MQ=255)
                                                                                      tAGGGGTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAg                                                                 <  1:791000/150‑1 (MQ=255)
                                                                                        gggTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCa                                                               <  1:868770/150‑1 (MQ=255)
                                                                                          gTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATc                                                             >  2:686828/1‑150 (MQ=255)
                                                                                                      ctggttagTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCTCTGGCAGATGGTTCAAATCTTTGTAGCTCTGAGCACGCGGCATTAATGTCGGTCTGGAGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACa                                                 <  2:191114/143‑1 (MQ=255)
                                                                                                              tCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCa                                         >  2:1305079/1‑150 (MQ=255)
                                                                                                                 tGAATGTTGTATACCGAAACAAAGATTTAACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAAcc                                      <  1:391695/150‑1 (MQ=255)
                                                                                                                  gAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCa                                     <  2:355849/150‑1 (MQ=255)
                                                                                                                      gttgtATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAg                                 >  2:1858903/1‑150 (MQ=255)
                                                                                                                      gttgtATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAg                                 <  2:1595184/150‑1 (MQ=255)
                                                                                                                           ataCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATa                            <  1:1238306/150‑1 (MQ=255)
                                                                                                                              gagctaCAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAgc                         >  1:1069073/6‑150 (MQ=25)
                                                                                                                              gagctaCAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAgc                         >  1:602617/6‑150 (MQ=25)
                                                                                                                              ccGAAACAAAGATTTTACCCGTCGGCGCTGGCTGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTTGGCGGGGCGTTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAgc                         >  2:302918/1‑150 (MQ=255)
                                                                                                                              ccGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAgc                         >  1:2026471/1‑150 (MQ=255)
                                                                                                                              ccGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGAAGCTCAGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAgc                         >  2:279098/1‑150 (MQ=255)
                                                                                                                                gctaCAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAgcgc                       >  1:486181/4‑150 (MQ=25)
                                                                                                                                 ctaCAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCt                      >  1:1432893/3‑150 (MQ=37)
                                                                                                                                  aaCAAAGATTTTACCCGTCGGCGCTGGCAGATGCGTCAAATCTTTGTAGCTCTGATCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGATAAGCATCAGCACAGTAAGATCACGGCAACCAAATGAAATAAGCTCTg                     >  1:313573/1‑150 (MQ=255)
                                                                                                                                    cAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCa                   >  1:875977/1‑150 (MQ=255)
                                                                                                                                    cAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCa                   >  1:916351/1‑150 (MQ=255)
                                                                                                                                    cAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCa                   >  2:1298134/1‑150 (MQ=255)
                                                                                                                                    cAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCTGGGCGCTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCa                   >  2:1086376/1‑150 (MQ=37)
                                                                                                                                     caaGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTGCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCAt                  <  2:354810/149‑1 (MQ=255)
                                                                                                                                              aCCCGGCGGCGCTGGCAGCTGGGTCAAATCTTTGTAGATCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTAtttt         <  2:1243892/150‑1 (MQ=255)
                                                                                                                                               cccGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTc        >  1:1172295/1‑150 (MQ=255)
                                                                                                                                                ccGGCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTcc       <  1:1842195/150‑1 (MQ=255)
                                                                                                                                                 cGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCt      >  1:896188/1‑150 (MQ=255)
                                                                                                                                                  gTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCtt     <  2:441307/150‑1 (MQ=255)
                                                                                                                                                   tCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTa    <  1:1998210/150‑1 (MQ=255)
                                                                                                                                                     ggCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATg  <  2:1989917/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTATAAACCAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATG  >  NC_000913/1101402‑1101700

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