Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,2610TC52.9% 4.5 / 97.9 70H71H (CAT→CAC) stfRputative prophage side tail fiber protein StfR
Reads supporting (aligned to +/- strand):  ref base T (21/12);  new base C (26/11);  total (47/23)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.29e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAGCCCGTGAGGCGGCAACCCATGCGGCTGATGCTGCG  >  NC_000913/1429019‑1429495
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                          
ggTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGAc                                                                                                                                                                                                                                    >  1:86022/1‑251 (MQ=255)
      ggCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTGCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACATAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGGGGTGGGGAACAGGCTGGCCTCAGAAAATCGGGATGAAGCCGGGCTTTATAGCATGGACGTTGAGTACGGTCATTACAGCGTTACTCTGTTGTTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTAc                                                                                                                                                                                                                              <  1:212371/251‑1 (MQ=255)
           gCTTTTTTGTGGGGTAAATATGGCAGTAAAGATTTCAGGTTAACGGAAAGACGGCACAGGAAAACCGGTACAGAACGGCAAAATCCAGCTGAAAGCAAAACGTAACAGCCCCACGGTTGTGGTGAACACGCTGGCCTCATAAAATCCGGATGAAGCCGGGCGTTACATCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCgtgt                                                                                                                                                                                                                         <  1:24294/251‑1 (MQ=255)
                   gtgGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGAt                                                                                                                                                                                                                 <  1:8352/251‑1 (MQ=255)
                     ggggTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAACCCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTACCAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATtc                                                                                                                                                                                                               <  2:63714/251‑1 (MQ=255)
                        gTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATtcac                                                                                                                                                                                                             >  1:190936/1‑248 (MQ=255)
                            atatGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGCACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAAccg                                                                                                                                                                                                        >  2:177041/1‑250 (MQ=255)
                              atGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAAccggg                                                                                                                                                                                                      <  1:119242/251‑4 (MQ=255)
                                  cAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACg                                                                                                                                                                                                  >  2:31399/1‑251 (MQ=255)
                                    gTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCt                                                                                                                                                                                                <  1:19093/251‑1 (MQ=255)
                                         gATTTCAGGTGTACTGAAAGACGGCACAGGAACACCGGTACAGTACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGGGGAAGGATTACCGCCGTCACATGCCGGGCCCATTACCGTGTATGAAGTTTCTCAACCCGGTACGCTGAATg                                                                                                                                                                                           >  2:259220/1‑251 (MQ=255)
                                                  agtccTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCGGCAGGTTGACGGTTTTCGACCATCGCACGCCGGGACCATCACCGTGTATGCAGATTCACAACAGGGGACGCTGAATGCTTTtctct                                                                                                                                                                                  >  2:261102/5‑250 (MQ=255)
                                                     ccTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTtctctgt                                                                                                                                                                                <  1:268630/249‑4 (MQ=255)
                                                       tGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGcc                                                                                                                                                                             >  2:281113/1‑251 (MQ=255)
                                                               ggCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACgg                                                                                                                                                                      >  1:336756/1‑250 (MQ=255)
                                                                gCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACgg                                                                                                                                                                      >  2:236010/1‑249 (MQ=255)
                                                                   caGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGGGgat                                                                                                                                                                 >  2:192830/1‑251 (MQ=255)
                                                                       aaaaCCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGatgatg                                                                                                                                                             >  2:8547/1‑251 (MQ=255)
                                                                           ccGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCg                                                                                                                                                         >  1:168910/1‑251 (MQ=255)
                                                                           acGGTACAGAACTGCACCATTCAGCTGAAAGCAAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGAAGGTTGACGGTTTTCCAACATCGCACGCCGGGACCATCACCGTGTATGAAGTTTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCg                                                                                                                                                         >  1:168911/2‑251 (MQ=255)
                                                                            cGGTACAGAACTGCACAATCAAGCTGAAAGCAAAACGAAACCGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGACGCAGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATTACGGAGGATGATGCCCGt                                                                                                                                                        >  2:181632/1‑251 (MQ=255)
                                                                               tACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCg                                                                                                                                                     >  1:252208/1‑251 (MQ=255)
                                                                                 cAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCccggccgg                                                                                                                                                    >  1:8018/1‑250 (MQ=255)
                                                                                   gAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAgg                                                                                                                                                 >  1:321511/1‑251 (MQ=255)
                                                                                   gAACTGCACAATCCAGCTGAAAGCAAAACGTAAAAGCACCACGGGGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAgg                                                                                                                                                 <  2:320419/251‑1 (MQ=255)
                                                                                    aaCTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAAc                                                                                                                                                                                                          <  1:351369/193‑1 (MQ=255)
                                                                                    aaCTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAAc                                                                                                                                                                                                          >  2:351369/1‑193 (MQ=255)
                                                                                    aaCTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAggt                                                                                                                                                <  2:86532/251‑2 (MQ=255)
                                                                                        gCACAACCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGAGTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTg                                                                                                                                            <  1:157979/251‑1 (MQ=255)
                                                                                          accattcAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAggtgctgcg                                                                                                                                          >  1:115521/7‑244 (MQ=255)
                                                                                             attcAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACACTCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGACTGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAggtgctgcgtcg                                                                                                                                       >  1:129026/4‑241 (MQ=255)
                                                                                             attcAGCTGAAAGCCAGACGTAACAGCACCACGGTAGAGGTGAACACGGTGGGATCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGAGGAGTACGGTCAGTACAGTGTCATCCTGCAGGATGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAggtgctgcgtcg                                                                                                                                       >  2:252336/4‑241 (MQ=255)
                                                                                             aTCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATTCCGGAGGATGATGCCCGTCCGGAGGCACTGCGCCg                                                                                                                                       >  2:192412/1‑251 (MQ=255)
                                                                                                 aGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGtttt                                                                                                                                   >  1:164849/1‑251 (MQ=255)
                                                                                                    tGAAAGCCAGACGTAACAGCACAACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCAACGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGa                                                                                                                                 >  1:129556/1‑250 (MQ=255)
                                                                                                     gAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCg                                                                                                                                                                                                       <  2:296333/179‑1 (MQ=255)
                                                                                                     gAAAGCAAAACGTAACAGCACCACGGTGGTGGCGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCg                                                                                                                                                                                                       >  1:296333/1‑179 (MQ=255)
                                                                                                                 tAACAGCACCACGGTGGGGGTGAACACGGTGGTCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGCCGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGCTGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGa                                                                                                                     >  2:91613/1‑249 (MQ=255)
                                                                                                                         ccacGGGGGGGGTGAACCGGCTGGCCTCAGAAAATCCGGAGGAAGCGTGGCGTTACAGCAGGGACGTTGATTACGGTCAGTACAGCGTTATTCGGTGGGGGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGgcg                                                                                                           <  2:318751/251‑1 (MQ=255)
                                                                                                                            cGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGTCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGCTGGTGGAAGAGGTGGCGCGt                                                                                                        >  2:132350/1‑251 (MQ=255)
                                                                                                                             ggtggtggtGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGACGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGt                                                                                                        >  2:76729/1‑250 (MQ=255)
                                                                                                                                  tggtgAACACGCTGGCATAAGACAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGGATGAAGATTCTCAACCCGGCACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGAGGAAGAGGTGGCGCGTAAagcg                                                                                                  >  2:215578/1‑251 (MQ=255)
                                                                                                                                   ggtgAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGt                                                                                                 <  1:281113/251‑1 (MQ=255)
                                                                                                                                      gactacGCTGGCGTGAGAAAATCCGGAGGAGGCGGGGCGTTACAGCATGGAGGTGGAGTACGCTCAGTACAGCGTTATGCTGTTGGTGGAAGGATTCCCGCCGGCACATGCCGGGACCATTACCGTGTATGAAGATTCGCAACCCGTTACGCTGAATGATTTTCTCGGTGCCATGACGGAGTATGAGGCCCGGCCGGAGGGACTGCGCCGGTTTGAACTGATGGTGGAAGAGGGGGCGAGTAACGCGTcc                                                                                               <  2:20417/246‑1 (MQ=255)
                                                                                                                                       aaCACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGTATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCcgt                                                                                             <  1:233291/251‑2 (MQ=255)
                                                                                                                                       aaCACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGGCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCcgt                                                                                             <  2:190936/251‑2 (MQ=255)
                                                                                                                                           cggtgggcacAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGAACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGCGGATGATGCCCGGCCAGCGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCgtggtg                                                                                         >  2:105976/10‑251 (MQ=255)
                                                                                                                                             gtgggcTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGACATCCTGCAGGATGACGGTTTTCCACCATCGCACGACGGGACCATCACCGTGTATGAAGATTCACAAACGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGAAGGTGCTGCGTCGTCTTGACCTGTTGGTGGAAGAGGTGGCGCGTAACGCGTCcgtgggggc                                                                                       >  2:190062/6‑244 (MQ=255)
                                                                                                                                             cTGGCATCAGCAAATTCGGGTGAAGCCGGGCCTTACAGCATGGACGCTGAGCACGGTCAGTACCGCGTCATTCAGTTCGTGGAAGGATTCCCGCCGTCACATGCCGGGAACATTACCGTGTATGAAGATTCTCACCCCAGTACGCTGAATGATTTTCTCGGCGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGAAACGCGTCCGCGGTGGc                                                                                       >  2:55603/1‑251 (MQ=255)
                                                                                                                                              tGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGc                                                                                       >  2:8412/1‑250 (MQ=255)
                                                                                                                                                       aaaaTCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCAAATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGCATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGAGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACg                                                                             >  1:148090/1‑251 (MQ=255)
                                                                                                                                                             cGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGcc                                                                       >  1:190553/1‑251 (MQ=255)
                                                                                                                                                             cGGATGAAGACGGGCGTTACAGCATGGATGTGGAGAACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCAACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGACTGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAgac                                                                       >  2:176546/1‑249 (MQ=255)
                                                                                                                                                              ggaggAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAgacg                                                                      <  2:276128/247‑4 (MQ=255)
                                                                                                                                                              ggATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAgacg                                                                      >  1:106096/1‑248 (MQ=255)
                                                                                                                                                                       cGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACTGTTTTACACCATCGCACGCCGGGACCATCACCTTTTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTACGTGGTGGCACAGAGTACGGCATACGCgaagaaat                                                             <  1:287257/251‑3 (MQ=255)
                                                                                                                                                                          gCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAg                                                          <  2:245726/251‑1 (MQ=255)
                                                                                                                                                                              tACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGccgg                                                      >  1:203506/1‑249 (MQ=255)
                                                                                                                                                                                    aTGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTCTGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGCCGGCGATGc                                                >  2:31913/1‑251 (MQ=255)
                                                                                                                                                                                     tggatgtggAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGCAATCAGCCGGCGATGcc                                               >  2:318646/9‑251 (MQ=255)
                                                                                                                                                                                       gatgtggAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGCAGAAATCAGCCGGCGATGcc                                               >  1:326441/7‑249 (MQ=255)
                                                                                                                                                                                              aGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGCCGGCGATGCCAgtgcatc                                       >  2:335518/1‑244 (MQ=255)
                                                                                                                                                                                               gTACGGTCAGTCCAGAGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGCCCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGAGGACGAGGTGGCGCGTACCGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGCAGTCAGCCAGTGATGCCAGCACATc                                       >  2:216535/1‑249 (MQ=255)
                                                                                                                                                                                                   ggTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGGGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCAGGGATGCCAGCACATCAGCCCg                                 >  2:198789/1‑251 (MQ=255)
                                                                                                                                                                                                      cAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAGCCCGTGa                              <  1:291443/251‑1 (MQ=255)
                                                                                                                                                                                                          aCAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGCCGGCGGATGATGCCCGTCCGGCGGCACTGCGCCGTTTTGAACTGATGGTGGCAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAAAACACGGCAGCCGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAACCCGTGAggcg                          >  1:344114/1‑251 (MQ=255)
                                                                                                                                                                                                            aGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCCCTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGCAGAAGTCAGCCAGTGATGCCAGCACATCAACCCGTGAggcggc                        >  1:113556/1‑251 (MQ=255)
                                                                                                                                                                                                                  aTTCTGTTGGTGGAAGTATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAGCCCGTGAGGCGGCAAccc                   <  1:31399/250‑1 (MQ=255)
                                                                                                                                                                                                                     atGTTGGTGGAAGGATTCCCGCCATCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCAGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGACGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGCACAAGGCAGCCGCGAAGAAGTCAGCCACTGATTCCAGCAAATCAGCCCGTGAGACGGCCACCAATGc               >  2:185470/2‑251 (MQ=255)
                                                                                                                                                                                                                                  gtttCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAGCCCGTGAGGCGGCAACCCATGCGGCTGATGCTGCg  <  1:198789/249‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                          
GGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAGCCCGTGAGGCGGCAACCCATGCGGCTGATGCTGCG  >  NC_000913/1429019‑1429495

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