Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,101,568 C→A 6.3% L39F (TTG→TTT)  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,5680CA 6.3% 293.9 / 3.2 94L39F (TTG→TTT) csgGcurli production assembly/transport outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (41/46);  new base A (2/5);  total (43/51)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGCTGGGGCTTAA  >  NC_000913/1101422‑1101714
                                                                                                                                                  |                                                                                                                                                  
ttCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGAT‑GGTCaaa                                                                                                                                                 <  2:854786/149‑1 (MQ=255)
     tGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCtt                                                                                                                                             <  2:1318344/149‑1 (MQ=37)
      gCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCttt                                                                                                                                            <  1:530566/149‑1 (MQ=37)
       ccGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTg                                                                                                                                           >  1:796738/1‑149 (MQ=255)
       ccGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTGTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTg                                                                                                                                           >  1:791461/1‑149 (MQ=255)
       ccGTGACCAGCATTGCCGTGGCGCTTTGAGTATCAGAATTTGAGAAGTTACTTGCCGTGTAGGGTTTAAATTGCCCGTTTTCGTACTGAATGTTGTATACCTAAACAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTg                                                                                                                                           <  2:1133068/149‑1 (MQ=16)
          tGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAg                                                                                                                                        >  1:888215/1‑149 (MQ=255)
            aCCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGct                                                                                                                                      >  2:1310940/1‑149 (MQ=255)
             ccAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTGtttcgg                                                                                                                                     <  2:316443/149‑6 (MQ=25)
              cAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGctct                                                                                                                                    >  2:238807/1‑149 (MQ=255)
              cAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGctct                                                                                                                                    >  2:825497/1‑149 (MQ=255)
              cAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCACGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGctct                                                                                                                                    >  1:1440911/1‑149 (MQ=255)
                gCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGa                                                                                                                                  <  2:701210/149‑1 (MQ=255)
                      ccGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGa                                                                                                                            >  2:573542/1‑149 (MQ=255)
                       cGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAg                                                                                                                           <  1:1494849/149‑1 (MQ=255)
                         tGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGTGTTAAAATTGCCCGGATTCGTCCTGAATGTTGTATACCGAAACAAAGAATTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCAGAGCACGAGGc                                                                                                                         >  1:1436090/1‑149 (MQ=255)
                            cgcTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCAtt                                                                                                                      <  2:536470/149‑1 (MQ=255)
                            cgcTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGGTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCAtt                                                                                                                      <  2:1289931/149‑1 (MQ=255)
                             gcTTTGCGGAACAGCAGTGGAGAAGTTACTTGTCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTa                                                                                                                     >  1:575127/1‑149 (MQ=255)
                               tttGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAAt                                                                                                                   <  1:927550/149‑1 (MQ=255)
                               tttGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAAt                                                                                                                   >  1:440670/1‑149 (MQ=255)
                                 tGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGt                                                                                                                 <  2:443652/149‑1 (MQ=255)
                                  gCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTc                                                                                                                <  2:1391983/149‑1 (MQ=255)
                                   cGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAg                                                                                                                           <  1:1157460/137‑1 (MQ=255)
                                   cGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAg                                                                                                                           >  2:1157460/1‑137 (MQ=255)
                                   cGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTATAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGTCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGAATTAATTTCg                                                                                                               >  1:966985/1‑149 (MQ=255)
                                                  taaGGTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACATAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGATCACGAGGCATTAATGTCGGTCTGGCGGcttctt                                                                                                <  1:967060/148‑1 (MQ=255)
                                                  gaAGTTACATGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGcttctt                                                                                                >  2:1532351/1‑149 (MQ=255)
                                                   aaGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCttt                                                                                               <  2:1312584/149‑1 (MQ=255)
                                                   aaGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCttt                                                                                               >  2:738446/1‑149 (MQ=255)
                                                    aGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTa                                                                                              <  1:665018/149‑1 (MQ=255)
                                                     gTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAg                                                                                             >  1:1242716/1‑149 (MQ=255)
                                                      ttACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTTAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAgg                                                                                            >  1:1110498/1‑149 (MQ=255)
                                                         cTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAggcgg                                                                                         <  1:899948/149‑1 (MQ=255)
                                                          ttGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAggcggg                                                                                        >  2:50406/1‑149 (MQ=255)
                                                               gtgtattgTTTAAATTTCCCGGTTTCGTCATGAATGTTGTACACCTAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGCGNNATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGt                                                                                   <  2:705264/142‑1 (MQ=255)
                                                                 gTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTa                                                                                 >  2:1275379/1‑149 (MQ=255)
                                                                  tAGTGTGTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTAACCGTCGGAGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                                <  1:295789/149‑1 (MQ=255)
                                                                  tAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                                <  1:672084/149‑1 (MQ=255)
                                                                  tAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                                <  1:764512/149‑1 (MQ=255)
                                                                  tAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                                <  1:573542/149‑1 (MQ=255)
                                                                  tAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                                <  1:1519687/149‑1 (MQ=255)
                                                                  tAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                                <  1:415099/149‑1 (MQ=255)
                                                                  tAGGGGTTAAATTGCCCGGTGTCGTCCTGACTGTTGTATACCGAAACAAAGATTTTAACCGTCGGCGCTGGCAGATGTGTCAAATCTTTGTAGCTCTGAGCAAGAGTCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                                <  1:262551/149‑1 (MQ=255)
                                                                      gTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCAt                                                                            <  2:888215/149‑1 (MQ=255)
                                                                      gTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCAt                                                                            >  2:709979/1‑149 (MQ=255)
                                                                      gTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCAt                                                                            >  2:293570/1‑149 (MQ=255)
                                                                       tttAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATc                                                                           >  1:197521/1‑149 (MQ=255)
                                                                        ttAAATTGCCCGGTTTCGTCCAGAATGTTGTATTCCGAAACAAAGATTTTACCAGTCGGCTCTTGCAGATGGGTCAAATCTTAGTAGCTCTGAGCACGAGGCATTAATGTAGGTCTGGCGGCTTCTTTAGGCGGGGCTGTTTAGAATcc                                                                          >  1:1390330/1‑149 (MQ=255)
                                                                          aaaTTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGc                                                                        <  2:1357662/149‑1 (MQ=255)
                                                                           aaTTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCCTCCGCt                                                                       >  1:1075117/1‑149 (MQ=255)
                                                                                 ccGGTTTCGCCCTGAATGTTGTATACCGAAACAAAGATTATACCCGGCGGCGCTGGCAGATGGGTCTAATCTTTGTAGATCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTaa                                                                 <  1:300232/149‑1 (MQ=255)
                                                                                        cGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACg                                                          <  2:1045149/149‑1 (MQ=255)
                                                                                         gTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACTAGGCATTAATGTCGGTCTGGCGCCTTCTTTAGGCGGGCCGGTTAAGCATCCGCTCAGTAACATGACgg                                                         >  2:1230366/1‑149 (MQ=255)
                                                                                         gTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACgg                                                         >  1:1335333/1‑149 (MQ=255)
                                                                                            cTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCaa                                                      >  1:38479/1‑149 (MQ=255)
                                                                                             tGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAAc                                                     <  2:791461/149‑1 (MQ=255)
                                                                                             tGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAAc                                                     >  2:249954/1‑149 (MQ=255)
                                                                                             tGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAAc                                                     <  2:796738/149‑1 (MQ=255)
                                                                                             tGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAAc                                                     >  1:1007928/1‑149 (MQ=255)
                                                                                                aTGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCCCGAGGCAGTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGCCGGCAACCaa                                                  >  2:190841/1‑149 (MQ=255)
                                                                                                      tataCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                                <  2:138511/113‑1 (MQ=255)
                                                                                                      tataCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                                >  1:138511/1‑113 (MQ=255)
                                                                                                       ataCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAAt                                           <  1:322013/149‑1 (MQ=255)
                                                                                                         aCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAATTCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTTAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAATTaa                                         <  2:1451993/149‑1 (MQ=255)
                                                                                                           cGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATTTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAgc                                       <  1:238807/149‑1 (MQ=255)
                                                                                                              taCAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCt                                    >  2:548515/2‑149 (MQ=37)
                                                                                                               aCAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTg                                   <  2:1563349/149‑1 (MQ=255)
                                                                                                               aCAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTg                                   <  2:1569502/149‑1 (MQ=255)
                                                                                                                 aaaGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTAGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCa                                 >  1:46539/1‑149 (MQ=255)
                                                                                                                 aaaGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCa                                 >  1:525367/1‑149 (MQ=37)
                                                                                                                  aaGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCAt                                >  1:1221277/1‑149 (MQ=255)
                                                                                                                         tACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGattatt                         <  1:577841/149‑1 (MQ=255)
                                                                                                                          aCCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTAttt                        <  2:147391/149‑1 (MQ=255)
                                                                                                                          aCCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTAttt                        <  1:825497/149‑1 (MQ=255)
                                                                                                                           accttCGGCGCTTGCAGATGGGTAAAATCTTTGTAGCTCTGAGCACGATGCATTAATGTCGGTATGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCATCCATAAGAAATAAGCGCTGCATGATTAtttt                       <  2:1390330/145‑1 (MQ=255)
                                                                                                                            ccGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTc                      <  2:923131/149‑1 (MQ=255)
                                                                                                                            ccGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTc                      <  2:38479/149‑1 (MQ=255)
                                                                                                                            ccGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTc                      <  2:1440911/149‑1 (MQ=255)
                                                                                                                            ccGTCGGCGCGGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTc                      <  1:288962/149‑1 (MQ=255)
                                                                                                                             cGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTcc                     >  1:1395016/1‑149 (MQ=255)
                                                                                                                              gTCTGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGGCACCAAAAGAAATAAGCGCTGCATGATTATTTTCCt                    <  1:25947/149‑1 (MQ=255)
                                                                                                                              gTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCGCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCt                    <  2:1007928/149‑1 (MQ=255)
                                                                                                                               tCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAAAGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCtt                   <  1:332031/149‑1 (MQ=255)
                                                                                                                                      tGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGc            >  1:1119161/1‑149 (MQ=255)
                                                                                                                                      tGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGc            >  1:253703/1‑149 (MQ=255)
                                                                                                                                        gCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGCTg          >  2:418796/1‑149 (MQ=255)
                                                                                                                                         cAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGCTgg         <  1:1275379/149‑1 (MQ=255)
                                                                                                                                          cgATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGCTggg        <  2:1168510/148‑1 (MQ=255)
                                                                                                                                           gATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGCTgggg       >  1:1332062/1‑149 (MQ=255)
                                                                                                                                            aTGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGCTGGGGc      >  2:766533/1‑149 (MQ=255)
                                                                                                                                            aTGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGCTGGGGc      >  2:762181/1‑149 (MQ=255)
                                                                                                                                               ggTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGCTGGGGCTTa   <  2:1335333/149‑1 (MQ=255)
                                                                                                                                                gTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGCTGGGGCTTaa  >  1:1437950/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGCTGGGGCTTAA  >  NC_000913/1101422‑1101714

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: ‑