Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,101,543 T→G 8.2% K48Q (AAA→CAA)  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,5430TG 8.2% 333.7 / 17.1 120K48Q (AAA→CAA) csgGcurli production assembly/transport outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (62/47);  new base G (6/4);  total (68/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 = 2.96e-01

Only 100 of 122 total aligned reads displayed.

CCAGCGGTATAAACCAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATT  >  NC_000913/1101396‑1101691
                                                                                                                                                   |                                                                                                                                                    
ccAGCGGTATACAGCAGCGAGAATCTTTCAGGGCCGTGACCCGCATTGCCGTGGCGATTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTAc                                                                                                                                                    <  1:393231/150‑1 (MQ=255)
ccAGCGGTATAAACCAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTAc                                                                                                                                                    <  2:1543212/150‑1 (MQ=255)
ccAGCGGTATAAACCAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTAc                                                                                                                                                    >  1:1323538/1‑150 (MQ=255)
ccAGCGGAATAAACCAGCGAGAATCTTTCAGTGCCGTGACCAGAATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGTGTTTAAATAGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTAc                                                                                                                                                    <  1:1234233/150‑1 (MQ=255)
 cAGCGGTATAAACCAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACTGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTAcc                                                                                                                                                   <  1:645452/150‑1 (MQ=255)
 cAGCGGTATAAACCAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTGAcc                                                                                                                                                   >  2:1402738/1‑150 (MQ=255)
   gCGGTATAAACCAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCg                                                                                                                                                 <  2:1806945/150‑1 (MQ=255)
      gTATAAACCAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCg                                                                                                                                              >  1:2041460/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)
              cAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGc                                                                                                                                      >  1:1397659/1‑150 (MQ=255)
              cAGCGAGAAGCTTTCAGAGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGc                                                                                                                                      <  2:523565/150‑1 (MQ=255)
               aGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCa                                                                                                                                     >  1:983113/1‑150 (MQ=255)
                gCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAg                                                                                                                                    >  2:1589047/1‑150 (MQ=255)
                gCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAg                                                                                                                                    <  2:1545153/150‑1 (MQ=255)
                gCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCACGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACACGTCGGCGCTGGCAg                                                                                                                                    >  1:1729005/1‑150 (MQ=255)
                   agaATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATg                                                                                                                                 <  2:1417309/150‑1 (MQ=255)
                   agaATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATg                                                                                                                                 >  2:1968770/1‑150 (MQ=255)
                     aaTCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATggg                                                                                                                               <  2:958233/150‑1 (MQ=255)
                       tCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTc                                                                                                                             <  1:689709/150‑1 (MQ=255)
                         tttCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGGCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCaa                                                                                                                           <  2:784672/150‑1 (MQ=255)
                           tCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAAt                                                                                                                         <  1:508065/150‑1 (MQ=255)
                           tCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAAt                                                                                                                         <  2:1378656/150‑1 (MQ=255)
                              gTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCtt                                                                                                                      >  1:181962/1‑150 (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)
                                        cAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTg                                                                                                            <  2:1095288/150‑1 (MQ=255)
                                               gCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGa                                                                                                     >  1:2023702/1‑150 (MQ=255)
                                               gCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGa                                                                                                     >  1:441307/1‑150 (MQ=255)
                                                ccGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAg                                                                                                    >  2:1707392/1‑150 (MQ=255)
                                                       gcTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTat                                                                                             >  1:374604/1‑149 (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                                                                                     <  2:618142/150‑1 (MQ=255)
                                                               gAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTc                                                                                     >  1:938030/1‑150 (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)
                                                                     cagTTGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCTAAACAAAGATTTTACCCGTCGGCGCTGTCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTggcgg                                                                               <  2:1836022/150‑1 (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:937872/1‑150 (MQ=255)
                                                                              aGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAg                                                                      >  1:1958343/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:589239/150‑1 (MQ=255)
                                                                                     tGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAggcggggc                                                               >  2:287667/1‑150 (MQ=255)
                                                                                     tGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAggcggggc                                                               >  1:432180/1‑150 (MQ=255)
                                                                                     tGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCGTTTggcggggc                                                               >  2:160210/1‑150 (MQ=255)
                                                                                       ccGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAggcggggcgg                                                             >  2:594922/1‑150 (MQ=255)
                                                                                        cGGGTAGGGTTTAATTTTCCCGGTGTCGTCCTGAATGTTGTATACCGAAACAAAGATATTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGt                                                            <  2:386377/150‑1 (MQ=255)
                                                                                        cGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGt                                                            >  2:1130069/1‑150 (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)
                                                                                                                        gAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCa                            <  2:355849/150‑1 (MQ=255)
                                                                                                                            gttgtATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAg                        >  2:1858903/1‑150 (MQ=255)
                                                                                                                            gttgtATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCATTGTAGCTCAGAGCACGAGGCATTAATGACGGACTGGCGGCTTCATAAGGAGGGGCGGTTAAGCATCCGCACACTTAGATAACGGCACCCAAAAg                        >  2:1329262/1‑150 (MQ=255)
                                                                                                                                 ataCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATa                   <  1:1238306/150‑1 (MQ=255)
                                                                                                                                  taCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATaa                  >  2:1963456/1‑150 (MQ=255)
                                                                                                                                    gagctaCAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAgc                >  1:602617/6‑150 (MQ=25)
                                                                                                                                    gagctaCAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAgc                >  1:1069073/6‑150 (MQ=25)
                                                                                                                                    ccGAAACAAAGATTTTACCCGTCGGCGCTGGCTGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTTGGCGGGGCGTTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAgc                >  2:302918/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          >  2:1298134/1‑150 (MQ=255)
                                                                                                                                          cAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCa          >  1:916351/1‑150 (MQ=255)
                                                                                                                                          cAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCa          >  1:875977/1‑150 (MQ=255)
                                                                                                                                          cAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCTGGGCGCTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCa          >  2:1086376/1‑150 (MQ=37)
                                                                                                                                           caaGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTGCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCAt         <  2:354810/149‑1 (MQ=255)
                                                                                                                                                  ttACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGattatt  <  1:1655583/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCAGCGGTATAAACCAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATT  >  NC_000913/1101396‑1101691

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: