Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,104,1750TG28.2% 147.0 / 13.4 71G75G (GGT→GGG) csgBcurlin, minor subunit
Reads supporting (aligned to +/- strand):  ref base T (16/35);  new base G (19/1);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 6.95e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGACACAACGTTAATTTCCATTCG  >  NC_000913/1103928‑1104416
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
aCGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TAGTTGCGCAAGAAGgtag                                                                                                                                                                                                                                                    <  1:218859/250‑1 (MQ=255)
    ccAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCa                                                                                                                                                                                                                                               <  1:764775/251‑1 (MQ=255)
      aGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCATAATATAACTTCGCGTTAAATGAATTGAGTAATTCTTCATTTAATCAGGCAGTCATAATTGTTCAAGTTTGGACTAATAATAGTGCTCAGTTAGGGCAGGGAGGCTCAAAACTGTTGGCGT‑TTGTTGCGCAGGAAGGGAGTTGCAAc                                                                                                                                                                                                                                             <  2:684437/251‑1 (MQ=255)
               gTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTTGCGG‑TTGTTTCGCAAGAAGGTAGTAGCAACCGGGCaa                                                                                                                                                                                                                                      >  1:331487/1‑249 (MQ=255)
                   caacaTGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAACAAGGGAGGAGCAACCGGGCAAAAATTg                                                                                                                                                                                                                                >  2:809883/1‑251 (MQ=255)
                    aacaTGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTg                                                                                                                                                                                                                                >  1:175409/1‑250 (MQ=255)
                     acaTGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGAc                                                                                                                                                                                                                              <  1:549078/251‑1 (MQ=255)
                     acaTGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGACTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAAATTTTATCGG‑GTGTTTCTCAAAGAGGGAGTAGGAACCGGGAAAAGATTggc                                                                                                                                                                                                                              >  2:174477/1‑249 (MQ=255)
                         gAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAg                                                                                                                                                                                                                           <  1:184762/250‑1 (MQ=255)
                           aaaaCAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACa                                                                                                                                                                                                                        <  1:416020/251‑1 (MQ=255)
                                       tatttatGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTACAAACTTTTGGGGG‑ATGTTTCGCCAGGAGGTAGTAGCAACAGGGCAAAAATTGAACAGAAAGGGGAtttttac                                                                                                                                                                                                            >  2:531472/1‑246 (MQ=255)
                                           tatGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCtt                                                                                                                                                                                                        <  2:838840/250‑1 (MQ=255)
                                            atgatgTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGc                                                                                                                                                                                                      <  1:614814/251‑1 (MQ=255)
                                               atgTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGGAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCat                                                                                                                                                                                                    >  1:164761/1‑250 (MQ=255)
                                                    aaCATTACTGGGGCCGCCTGGGATTGCAGCCTCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTGCGCGTTAAATGAATTGAGTAATTCTTCATTTAAGCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAACCTTTTGGCGG‑TTGTTGCGCAAGAATGTAGTAGCAACCGGGCAAATATTGTCCAGTCAGGAGATTATAT‑CCTTGCATATATTg                                                                                                                                                                                              <  2:212048/251‑1 (MQ=255)
                                                     aCAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAAATTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAAGAGATTATAA‑ACTTTCCTATAtttt                                                                                                                                                                                             >  2:273859/1‑249 (MQ=255)
                                                      cAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGGAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGAt                                                                                                                                                                                            >  1:305318/1‑251 (MQ=255)
                                                       aaTACTGGGTGCGCCTGGGATTGCAGACGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGGTCAAAAATTTTGGCGG‑TTGTTGCGCAAGAAGGGAGTAGCAACCGGGCCAAGAATGACCAGACCGGGGATTATAA‑CCTTGCATATATTGAtt                                                                                                                                                                                           >  1:885759/1‑250 (MQ=255)
                                                              ggTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCggg                                                                                                                                                                                    <  2:175409/251‑1 (MQ=255)
                                                                  cgcCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGt                                                                                                                                                                                <  1:200103/251‑1 (MQ=255)
                                                                   gcCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑ATTTTGCGCAAGAAGGGAGTAGCAACCGGGCAAAGATTGACCACACGGGAGGTTATAA‑ACTTTCATTTTTTGATCAGGCGGGCAggg                                                                                                                                                                               >  2:100127/1‑249 (MQ=255)
                                                                    ccTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGAACAGACAGGAGATTATAA‑CCTTGCATATATTGATAAGGCGGGCAGTGc                                                                                                                                                                              >  1:274144/1‑251 (MQ=255)
                                                                       gggATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCa                                                                                                                                                                            >  1:535418/1‑250 (MQ=255)
                                                                               aGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCa                                                                                                                                                                   <  1:273859/251‑1 (MQ=255)
                                                                                  cGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTa                                                                                                                                                                <  2:9625/251‑1 (MQ=255)
                                                                                     agcagGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGGAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTAttt                                                                                                                                                             >  1:642080/1‑251 (MQ=255)
                                                                                     agcagGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGGAGTAGCAACCGGGCAAAGATTGACCAGACAAGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTAttt                                                                                                                                                             >  1:444322/1‑251 (MQ=255)
                                                                                      gcagGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTc                                                                                                                                                            <  1:100127/251‑1 (MQ=255)
                                                                                          gTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTTGGCAAAAAGGGAGTAGCAAACGGGCAAAAAATTACAAGAAAAGAGACTATAA‑ACTTGCATATATTTATCAAGCGGGCAgggaccacgaatccagtttaggta                                                                                                                                                          >  2:491826/1‑226 (MQ=255)
                                                                                          gTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTAt                                                                                                                                                               <  1:491826/244‑1 (MQ=255)
                                                                                                    aGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGGAGGAGGAAACGGGCAAAGAATGACCAGACAGGGGGTTATAA‑CCTTGCCTATATTTAACAGGCGGGGAGTGGCAACGAAGCCAGGATTTTGCaaaggggttttt                                                                                                                                              >  2:90925/1‑241 (MQ=255)
                                                                                                     gCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAAGCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTTAGTTACGGCAGGGAGGCTCAAAACTTGTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTTACCAGACAGGAGATTATAA‑CCTTGCATATATTGAGCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATgg                                                                                                                                             <  2:201161/251‑1 (MQ=255)
                                                                                                        aaTTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGTCAGTATTTCGCAAGGTGCTTATGGTaa                                                                                                                                          <  2:670770/251‑1 (MQ=255)
                                                                                                        aaTTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGt                                                                                                                                            <  1:248084/249‑1 (MQ=255)
                                                                                                            cAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACt                                                                                                                                      <  1:121188/251‑1 (MQ=255)
                                                                                                                 tataACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAgg                                                                                                                                                                                                                      <  2:30094/167‑1 (MQ=255)
                                                                                                                 tataACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGGAGTAGCAACCGGGCAAAGATTGACCAGACAgg                                                                                                                                                                                                                      >  1:30094/1‑167 (MQ=255)
                                                                                                                     aCTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑Att                                                                                                                            <  1:426635/251‑1 (MQ=255)
                                                                                                                          gcgGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCa                                                                                                                       <  1:614088/251‑1 (MQ=255)
                                                                                                                                 aTGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGt                                                                                                                <  2:305318/251‑1 (MQ=255)
                                                                                                                                   gAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGtt                                                                                                               <  1:533501/250‑1 (MQ=255)
                                                                                                                                   gAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCCCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGGAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTTATCAGGCGGGCAGTGACAACGATGCCAGGATTTCGCAAGGTGCTTATGGGAATACTGCGATG‑ATTATCCAGAAAGGttt                                                                                                              >  1:640308/1‑250 (MQ=255)
                                                                                                                                       tGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCCAGGTGCTTATGGTAATACTGCGGTG‑ATTATCCAGAAAGGTTCT‑GGt                                                                                                         >  1:2723/1‑251 (MQ=255)
                                                                                                                                            aaGTCTTCATTTAATCAGGCAGCCATAATTTGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGtaataa                                                                                                    <  1:91222/251‑1 (MQ=255)
                                                                                                                                              gTCTTCATTTAAGCAGGCAGCCAGAATTGGGCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGAGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAg                                                                                                  <  2:137659/251‑1 (MQ=255)
                                                                                                                                                   cATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAat                                                                                             <  1:703054/251‑1 (MQ=255)
                                                                                                                                                    aTTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAata                                                                                            <  1:90925/251‑1 (MQ=255)
                                                                                                                                                        aaTCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGGGGCTTAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGAAGCAACCCGGCAAAAAATGACCAGAAAGGAGATTATAAACCTTG‑ATAATTTGATCAGGCGGGCGGGGCCCAAGATGCCAGGTTTTTTCAAGGGGCTTATGGTAATACGTCGAGG‑ATTTTCCAGAAAGGTTTT‑GGGAATAAAgacaaaattac                                                                                        >  2:388580/1‑241 (MQ=255)
                                                                                                                                                                aGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGGGCAAGAAGGGAGTAGCAACCGGGCAAAGATTGACCCGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGGATTTCGCAAGGGGCTTATGGTAATAATTCGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAAAATTAAACAGTATg                                                                                >  1:514344/1‑251 (MQ=255)
                                                                                                                                                                 gCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTTCGCAAGAAGGTAGTAGCAAACGGGCCAAGATTGGCCAGACAGGGGATTATAA‑CCTTGCATATATTTATCAGGCGGGCAGGGCCCACGATGCCAGTATTTTGCAAGGTGCTTTTGGGAATACTTCGATG‑ATTTTCCCGAAAGGTTCTTGGTA‑TAAAGCAAAATTTTCACAGTAttg                                                                               >  2:82066/1‑249 (MQ=255)
                                                                                                                                                                 gCCATAATTGGACAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGATTTTTGC‑CAAAAAGGGAGTAGCAACCGGGCAAAAAATGACCAAACAGGAGGTTTTAA‑CCTTTCATATATTGATAAGGGGGGCAGTGCCAACCGTGCCAGTATTTCGCAAGGTGCTTATTGTAAAACTGCGATGTATT‑TCCCAAAAGGGTAT‑GGGAAAAAAGCAAATATTACAAAGTaagg                                                                               >  2:500118/1‑248 (MQ=255)
                                                                                                                                                                  ccATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGGAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTTATCAGGCGGGCAGGGGCAACCATGCCAGTATTTTGCAAGGGGCTTATGGTAAAACTGCGATG‑AGTATCCAGAAAGGTTTT‑GGTAATAAAGCAAATATTACACAGTAttgg                                                                              >  2:417299/1‑248 (MQ=255)
                                                                                                                                                                    aTAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAATATTACACAGTATGGTAc                                                                            <  1:468125/251‑1 (MQ=255)
                                                                                                                                                                                   gggACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCGTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCTATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCaa                                                              <  2:585498/250‑1 (MQ=255)
                                                                                                                                                                                   gggACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGGAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAAt                                                             >  1:739479/1‑251 (MQ=255)
                                                                                                                                                                                     gACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAAACGGGCAAAGGTTGACCAGACAGGAGATTATAA‑CCTTGCCTATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGGGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAATATTACACAGGATGGTACTCAAAAAAAGGCAATTg                                                           >  1:248078/1‑251 (MQ=255)
                                                                                                                                                                                        taataaTAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGGAGTAGGAACCGGGCAAAGATTGACCAGAAAGGAGGTTATAA‑CCTTGCATATATTGATCAGGCGGGCAGGGCCAACGATGCCAGTATTTAGCAAGGGGCTTATGGGAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGAAATAAAACAAATAATACAAAGTATGGTACTCAAAAAACGGCAATTgg                                                          >  1:86840/1‑248 (MQ=255)
                                                                                                                                                                                           taataGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGc                                                     >  1:416188/1‑251 (MQ=255)
                                                                                                                                                                                               aGTGCTCAGTTACGGCAGGGAAGCTCAAAACTTTTGGCGG‑GTGTTGGGCAAGAAGGGAGGAGCAACCGGGCAAAGATTGGACAGGCAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGGAATACTGCGATG‑ATTATCCAGAAAGGGTCT‑GGGAATAAAGCAAATATTAAACCGGATGGTAATCAAAAAAACGCAATTGTATTGCag                                                   >  1:12165/1‑249 (MQ=255)
                                                                                                                                                                                                     cAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑GTGTTGCGCAAGAAGGGAGGAGCAACCGGGCAAAGATTGGCCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGGGCCAACGATGCCAGTATTTCGCAAGGGGCTTATGGTAATAACGCGATG‑ATTATCCCGAAAGGGTCT‑GGGAATAAAGCAAATATTACACAGTATGGTACTCAAAAAAAGGCAATTGGAGTGCAGAGACAGTc                                           >  2:265707/1‑251 (MQ=255)
                                                                                                                                                                                                      aGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCg                                          <  1:17399/251‑1 (MQ=255)
                                                                                                                                                                                                               aGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCt                                 <  1:862537/251‑1 (MQ=255)
                                                                                                                                                                                                                 ggAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTAt                               <  2:535418/251‑1 (MQ=255)
                                                                                                                                                                                                                 ggAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTAt                               >  2:288003/1‑251 (MQ=255)
                                                                                                                                                                                                                 ggAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTAt                               >  1:329839/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ggCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGTCAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGACACAACGtt              <  1:871921/251‑1 (MQ=255)
                                                                                                                                                                                                                                      g‑ttgttGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGATATTATAA‑CCTTGCATATATGGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAATCGGCAATTGTAGTGCAGAGACAGTGGCAAATGGATATTCGCTTGACACAACGTTAAtt          <  2:755456/251‑1 (MQ=255)
                                                                                                                                                                                                                                      g‑ttgttGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGc                                                                >  2:515724/1‑197 (MQ=255)
                                                                                                                                                                                                                                         tgttgCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGACACAACGTTAAttt         <  2:2723/250‑1 (MQ=255)
                                                                                                                                                                                                                                             gcgcAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAC‑CCTTGCATATAGTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGACACAACGTTAATTTCCAtt    <  2:652681/251‑1 (MQ=255)
                                                                                                                                                                                                                                               gcAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGACACAACGTTAATTTCCATTCg  >  2:654237/1‑251 (MQ=255)
                                                                                                                                                                                                                                               gcAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTAGCGTGACACAACGTTAATTTCCATTCg  >  1:704090/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
ACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGG‑TTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAA‑CCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATG‑ATTATCCAGAAAGGTTCT‑GGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGACACAACGTTAATTTCCATTCG  >  NC_000913/1103928‑1104416

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑