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→GGGcsgBcurlin, 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.

ACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGACACAACGTTAATTTCCATTCG  >  NC_000913/1103928‑1104416
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
aCGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTAGTTGCGCAAGAAGgtag                                                                                                                                                                                                                                                    <  1:218859/250‑1 (MQ=255)
    ccAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCa                                                                                                                                                                                                                                               <  1:764775/251‑1 (MQ=255)
      aGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCATAATATAACTTCGCGTTAAATGAATTGAGTAATTCTTCATTTAATCAGGCAGTCATAATTGTTCAAGTTTGGACTAATAATAGTGCTCAGTTAGGGCAGGGAGGCTCAAAACTGTTGGCGTTTGTTGCGCAGGAAGGGAGTTGCAAc                                                                                                                                                                                                                                             <  2:684437/251‑1 (MQ=255)
               gTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTTGCGGTTGTTTCGCAAGAAGGTAGTAGCAACCGGGCaa                                                                                                                                                                                                                                      >  1:331487/1‑249 (MQ=255)
                   caacaTGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAACAAGGGAGGAGCAACCGGGCAAAAATTg                                                                                                                                                                                                                                >  2:809883/1‑251 (MQ=255)
                    aacaTGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTg                                                                                                                                                                                                                                >  1:175409/1‑250 (MQ=255)
                     acaTGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGAc                                                                                                                                                                                                                              <  1:549078/251‑1 (MQ=255)
                     acaTGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGACTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAAATTTTATCGGGTGTTTCTCAAAGAGGGAGTAGGAACCGGGAAAAGATTggc                                                                                                                                                                                                                              >  2:174477/1‑249 (MQ=255)
                         gAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAg                                                                                                                                                                                                                           <  1:184762/250‑1 (MQ=255)
                           aaaaCAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACa                                                                                                                                                                                                                        <  1:416020/251‑1 (MQ=255)
                                       tatttatGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTACAAACTTTTGGGGGATGTTTCGCCAGGAGGTAGTAGCAACAGGGCAAAAATTGAACAGAAAGGGGAtttttac                                                                                                                                                                                                            >  2:531472/1‑246 (MQ=255)
                                           tatGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCtt                                                                                                                                                                                                        <  2:838840/250‑1 (MQ=255)
                                            atgatgTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGc                                                                                                                                                                                                      <  1:614814/251‑1 (MQ=255)
                                               atgTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGGAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCat                                                                                                                                                                                                    >  1:164761/1‑250 (MQ=255)
                                                    aaCATTACTGGGGCCGCCTGGGATTGCAGCCTCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTGCGCGTTAAATGAATTGAGTAATTCTTCATTTAAGCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAACCTTTTGGCGGTTGTTGCGCAAGAATGTAGTAGCAACCGGGCAAATATTGTCCAGTCAGGAGATTATATCCTTGCATATATTg                                                                                                                                                                                              <  2:212048/251‑1 (MQ=255)
                                                     aCAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAAATTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAAGAGATTATAAACTTTCCTATAtttt                                                                                                                                                                                             >  2:273859/1‑249 (MQ=255)
                                                      cAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGGAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGAt                                                                                                                                                                                            >  1:305318/1‑251 (MQ=255)
                                                       aaTACTGGGTGCGCCTGGGATTGCAGACGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGGTCAAAAATTTTGGCGGTTGTTGCGCAAGAAGGGAGTAGCAACCGGGCCAAGAATGACCAGACCGGGGATTATAACCTTGCATATATTGAtt                                                                                                                                                                                           >  1:885759/1‑250 (MQ=255)
                                                              ggTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCggg                                                                                                                                                                                    <  2:175409/251‑1 (MQ=255)
                                                                  cgcCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGt                                                                                                                                                                                <  1:200103/251‑1 (MQ=255)
                                                                   gcCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGATTTTGCGCAAGAAGGGAGTAGCAACCGGGCAAAGATTGACCACACGGGAGGTTATAAACTTTCATTTTTTGATCAGGCGGGCAggg                                                                                                                                                                               >  2:100127/1‑249 (MQ=255)
                                                                    ccTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGAACAGACAGGAGATTATAACCTTGCATATATTGATAAGGCGGGCAGTGc                                                                                                                                                                              >  1:274144/1‑251 (MQ=255)
                                                                       gggATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCa                                                                                                                                                                            >  1:535418/1‑250 (MQ=255)
                                                                               aGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCa                                                                                                                                                                   <  1:273859/251‑1 (MQ=255)
                                                                                  cGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTa                                                                                                                                                                <  2:9625/251‑1 (MQ=255)
                                                                                     agcagGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGGAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTAttt                                                                                                                                                             >  1:642080/1‑251 (MQ=255)
                                                                                     agcagGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGGAGTAGCAACCGGGCAAAGATTGACCAGACAAGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTAttt                                                                                                                                                             >  1:444322/1‑251 (MQ=255)
                                                                                      gcagGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTc                                                                                                                                                            <  1:100127/251‑1 (MQ=255)
                                                                                          gTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTTGGCAAAAAGGGAGTAGCAAACGGGCAAAAAATTACAAGAAAAGAGACTATAAACTTGCATATATTTATCAAGCGGGCAgggaccacgaatccagtttaggta                                                                                                                                                          >  2:491826/1‑226 (MQ=255)
                                                                                          gTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTAt                                                                                                                                                               <  1:491826/244‑1 (MQ=255)
                                                                                                    aGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGGAGGAGGAAACGGGCAAAGAATGACCAGACAGGGGGTTATAACCTTGCCTATATTTAACAGGCGGGGAGTGGCAACGAAGCCAGGATTTTGCaaaggggttttt                                                                                                                                              >  2:90925/1‑241 (MQ=255)
                                                                                                     gCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAAGCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTTAGTTACGGCAGGGAGGCTCAAAACTTGTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTTACCAGACAGGAGATTATAACCTTGCATATATTGAGCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATgg                                                                                                                                             <  2:201161/251‑1 (MQ=255)
                                                                                                        aaTTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGTCAGTATTTCGCAAGGTGCTTATGGTaa                                                                                                                                          <  2:670770/251‑1 (MQ=255)
                                                                                                        aaTTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGt                                                                                                                                            <  1:248084/249‑1 (MQ=255)
                                                                                                            cAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACt                                                                                                                                      <  1:121188/251‑1 (MQ=255)
                                                                                                                 tataACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAgg                                                                                                                                                                                                                      <  2:30094/167‑1 (MQ=255)
                                                                                                                 tataACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGGAGTAGCAACCGGGCAAAGATTGACCAGACAgg                                                                                                                                                                                                                      >  1:30094/1‑167 (MQ=255)
                                                                                                                     aCTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGAtt                                                                                                                            <  1:426635/251‑1 (MQ=255)
                                                                                                                          gcgGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCa                                                                                                                       <  1:614088/251‑1 (MQ=255)
                                                                                                                                 aTGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGt                                                                                                                <  2:305318/251‑1 (MQ=255)
                                                                                                                                   gAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGtt                                                                                                               <  1:533501/250‑1 (MQ=255)
                                                                                                                                   gAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCCCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGGAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTTATCAGGCGGGCAGTGACAACGATGCCAGGATTTCGCAAGGTGCTTATGGGAATACTGCGATGATTATCCAGAAAGGttt                                                                                                              >  1:640308/1‑250 (MQ=255)
                                                                                                                                       tGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCCAGGTGCTTATGGTAATACTGCGGTGATTATCCAGAAAGGTTCTGGt                                                                                                         >  1:2723/1‑251 (MQ=255)
                                                                                                                                            aaGTCTTCATTTAATCAGGCAGCCATAATTTGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGtaataa                                                                                                    <  1:91222/251‑1 (MQ=255)
                                                                                                                                              gTCTTCATTTAAGCAGGCAGCCAGAATTGGGCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGAGATGATTATCCAGAAAGGTTCTGGTAATAAAg                                                                                                  <  2:137659/251‑1 (MQ=255)
                                                                                                                                                   cATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAat                                                                                             <  1:703054/251‑1 (MQ=255)
                                                                                                                                                    aTTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAata                                                                                            <  1:90925/251‑1 (MQ=255)
                                                                                                                                                        aaTCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGGGGCTTAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGAAGCAACCCGGCAAAAAATGACCAGAAAGGAGATTATAAACCTTGATAATTTGATCAGGCGGGCGGGGCCCAAGATGCCAGGTTTTTTCAAGGGGCTTATGGTAATACGTCGAGGATTTTCCAGAAAGGTTTTGGGAATAAAgacaaaattac                                                                                        >  2:388580/1‑241 (MQ=255)
                                                                                                                                                                aGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGGGCAAGAAGGGAGTAGCAACCGGGCAAAGATTGACCCGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGGATTTCGCAAGGGGCTTATGGTAATAATTCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAAAATTAAACAGTATg                                                                                >  1:514344/1‑251 (MQ=255)
                                                                                                                                                                 gCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTTCGCAAGAAGGTAGTAGCAAACGGGCCAAGATTGGCCAGACAGGGGATTATAACCTTGCATATATTTATCAGGCGGGCAGGGCCCACGATGCCAGTATTTTGCAAGGTGCTTTTGGGAATACTTCGATGATTTTCCCGAAAGGTTCTTGGTATAAAGCAAAATTTTCACAGTAttg                                                                               >  2:82066/1‑249 (MQ=255)
                                                                                                                                                                 gCCATAATTGGACAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGATTTTTGCCAAAAAGGGAGTAGCAACCGGGCAAAAAATGACCAAACAGGAGGTTTTAACCTTTCATATATTGATAAGGGGGGCAGTGCCAACCGTGCCAGTATTTCGCAAGGTGCTTATTGTAAAACTGCGATGTATTTCCCAAAAGGGTATGGGAAAAAAGCAAATATTACAAAGTaagg                                                                               >  2:500118/1‑248 (MQ=255)
                                                                                                                                                                  ccATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGGAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTTATCAGGCGGGCAGGGGCAACCATGCCAGTATTTTGCAAGGGGCTTATGGTAAAACTGCGATGAGTATCCAGAAAGGTTTTGGTAATAAAGCAAATATTACACAGTAttgg                                                                              >  2:417299/1‑248 (MQ=255)
                                                                                                                                                                    aTAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTAc                                                                            <  1:468125/251‑1 (MQ=255)
                                                                                                                                                                                   gggACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCGTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCTATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCaa                                                              <  2:585498/250‑1 (MQ=255)
                                                                                                                                                                                   gggACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGGAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAAt                                                             >  1:739479/1‑251 (MQ=255)
                                                                                                                                                                                     gACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAAACGGGCAAAGGTTGACCAGACAGGAGATTATAACCTTGCCTATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGGGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGGATGGTACTCAAAAAAAGGCAATTg                                                           >  1:248078/1‑251 (MQ=255)
                                                                                                                                                                                        taataaTAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGGAGTAGGAACCGGGCAAAGATTGACCAGAAAGGAGGTTATAACCTTGCATATATTGATCAGGCGGGCAGGGCCAACGATGCCAGTATTTAGCAAGGGGCTTATGGGAATACTGCGATGATTATCCAGAAAGGTTCTGGAAATAAAACAAATAATACAAAGTATGGTACTCAAAAAACGGCAATTgg                                                          >  1:86840/1‑248 (MQ=255)
                                                                                                                                                                                           taataGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGc                                                     >  1:416188/1‑251 (MQ=255)
                                                                                                                                                                                               aGTGCTCAGTTACGGCAGGGAAGCTCAAAACTTTTGGCGGGTGTTGGGCAAGAAGGGAGGAGCAACCGGGCAAAGATTGGACAGGCAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGGAATACTGCGATGATTATCCAGAAAGGGTCTGGGAATAAAGCAAATATTAAACCGGATGGTAATCAAAAAAACGCAATTGTATTGCag                                                   >  1:12165/1‑249 (MQ=255)
                                                                                                                                                                                                     cAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGGTGTTGCGCAAGAAGGGAGGAGCAACCGGGCAAAGATTGGCCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGGGCCAACGATGCCAGTATTTCGCAAGGGGCTTATGGTAATAACGCGATGATTATCCCGAAAGGGTCTGGGAATAAAGCAAATATTACACAGTATGGTACTCAAAAAAAGGCAATTGGAGTGCAGAGACAGTc                                           >  2:265707/1‑251 (MQ=255)
                                                                                                                                                                                                      aGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCg                                          <  1:17399/251‑1 (MQ=255)
                                                                                                                                                                                                               aGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCt                                 <  1:862537/251‑1 (MQ=255)
                                                                                                                                                                                                                 ggAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTAt                               <  2:535418/251‑1 (MQ=255)
                                                                                                                                                                                                                 ggAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTAt                               >  2:288003/1‑251 (MQ=255)
                                                                                                                                                                                                                 ggAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTAt                               >  1:329839/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ggCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGTCAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGACACAACGtt              <  1:871921/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gttgttGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGATATTATAACCTTGCATATATGGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAATCGGCAATTGTAGTGCAGAGACAGTGGCAAATGGATATTCGCTTGACACAACGTTAAtt          <  2:755456/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gttgttGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGc                                                                >  2:515724/1‑197 (MQ=255)
                                                                                                                                                                                                                                         tgttgCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGACACAACGTTAAttt         <  2:2723/250‑1 (MQ=255)
                                                                                                                                                                                                                                             gcgcAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATACCCTTGCATATAGTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGACACAACGTTAATTTCCAtt    <  2:652681/251‑1 (MQ=255)
                                                                                                                                                                                                                                               gcAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGACACAACGTTAATTTCCATTCg  >  2:654237/1‑251 (MQ=255)
                                                                                                                                                                                                                                               gcAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTAGCGTGACACAACGTTAATTTCCATTCg  >  1:704090/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
ACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAGTCTTCATTTAATCAGGCAGCCATAATTGGTCAAGCTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGTGACACAACGTTAATTTCCATTCG  >  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: