Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,2730TC53.6% 9.7 / 105.0 69I75I (ATT→ATC) stfRputative prophage side tail fiber protein StfR
Reads supporting (aligned to +/- strand):  ref base T (20/12);  new base C (26/11);  total (46/23)
Fisher's exact test for biased strand distribution p-value = 6.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAGCCCGTGAGGCGGCAACCCATGCGGCTGATGCTGCG  >  NC_000913/1429025‑1429495
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                              
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                                                                                                                                                                                                          >  2:351369/1‑193 (MQ=255)
                                                                              aaCTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAAc                                                                                                                                                                                                          <  1:351369/193‑1 (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)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                              
GGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAGCCCGTGAGGCGGCAACCCATGCGGCTGATGCTGCG  >  NC_000913/1429025‑1429495

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