Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913494,5500GC28.9% 81.6 / 11.6 45R49P (CGG→CCG) recRrecombination mediator protein RecR
Reads supporting (aligned to +/- strand):  ref base G (5/27);  new base C (13/0);  total (18/27)
Fisher's exact test for biased strand distribution p-value = 1.17e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGG  >  NC_000913/494309‑494789
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
cgacgCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTc                                                                                                                                                                                                                                           <  2:338008/251‑1 (MQ=255)
           ggtgtATTGAAGAAACGCAGAAAGAATAAAGGGCCTCTGTAGCCGCCGGAATGCAGCTGCCGGCTGGCTTTAAGATGACGTTTTGACGCAAACCAGCCTGATGTTAACACAGCTTATGGAAGCACTGGGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCATCGTATTGAGTTCACGCTGCTGCATCGCGATCTGAGTGGCGGGATGCTTCGGGCGCATGCGCTCACCCGGGCGATGGCGGAAATC‑GGcc                                                                                                                                                                                                                               <  2:277510/248‑1 (MQ=255)
                aTTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGAGGCAACCCAGCCCGCTGTTAGCACAGCTTAGGGAAGCACTGCGCTGGCTGCCGGGCGTGGGCCTGAAGACGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGCTGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTg                                                                                                                                                                                                                          <  2:109251/250‑1 (MQ=255)
                                  gaaaAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑c                                                                                                                                                                                                      <  1:228155/251‑1 (MQ=255)
                                     aaaaTGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CAcc                                                                                                                                                                                                   <  1:430655/251‑1 (MQ=255)
                                      aaaTGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGATTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGGAGCGTCTGGGGGCGGCGCCCCCCCCGGGGGTGTCGGAAAAC‑GGGC‑ACTGCGcccagtgccccccctttcccg                                                                                                                                                                                                   >  1:387349/1‑232 (MQ=255)
                                           gCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAg                                                                                                                                                                                             <  1:25157/251‑1 (MQ=255)
                                                 gTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGGGGCGTCTGGGGCAGGCGCTCACCCCGGCGGTGTCGGGAATC‑TGCC‑ACCGCGCCGATTTCCGCCCTTT‑CACCGACAAGGAAGTc                                                                                                                                                                                       >  1:87166/1‑251 (MQ=255)
                                                  tATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGGTGCGTCTTGGGCAGGGGCTCCCCCCGGGGATGTCGGAAAAC‑GGCC‑CCCGCGCCGCTTTCCGCCCTTT‑CCCCGAACAGGaaagct                                                                                                                                                                                      >  1:453055/1‑247 (MQ=255)
                                                   aTCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCt                                                                                                                                                                                      <  1:238042/250‑1 (MQ=255)
                                                     cctccGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGGGTCTGGCGCAGGCGCGCAACCCGGCGGTGTCGGAAATC‑GGCCCACT‑CGCCGGTTTCCGCACTTT‑CACCGAACAGGAAgacttt                                                                                                                                                                                    >  1:274503/1‑245 (MQ=255)
                                                       tccGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGGTGCGTCTGGGGCAGGCGCCCCCCCCGGCGCTGTCGGAAAAC‑GGCC‑CCCGCGCCGCTTGCCGCCCTTT‑CCCCCAACACGAAGGCTGTaa                                                                                                                                                                                  >  2:177611/1‑250 (MQ=255)
                                                        ccGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCTGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGGCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCTGTaa                                                                                                                                                                                  <  2:43318/249‑1 (MQ=255)
                                                         cGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCTGTAACAt                                                                                                                                                                               <  2:315064/251‑1 (MQ=255)
                                                                 aGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGGCTGGGGCAGGCGCTCCCCCCGGCGGTGGCGGAAAACAGGCC‑‑CTGCGCCCGGTGCCGCACGTT‑CACCGGACAGGGAGTCTGTAAAATCTTTTccg                                                                                                                                                                       >  1:14963/1‑249 (MQ=255)
                                                                             gCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCgtcgtc                                                                                                                                                           <  1:315499/251‑1 (MQ=255)
                                                                                          gTTCTGCTGCAAACCAGCCCGCTGTCAACACCGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGGTGCGGCTGGCGCAGGCGCTCCCCCCGGGGATGTCGGGAATC‑GGCC‑CCCGCGCCGATTGCCGCCCTTT‑CACCGGACAGGAAGGCTGTAACATCTGGTTGGATCCGCgtcgtccggcaaaagggc                                                                                                                                               >  2:62479/1‑238 (MQ=255)
                                                                                               gATGCAATCCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGTATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑GCTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAAt                                                                                                                                           <  1:264787/249‑1 (MQ=255)
                                                                                                           cccGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCTTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCGCACCCGTGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTTCCGCACTTT‑CACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCTCTTCGTCAGGAAAACGGTCAAATCTGCGTGGTGgaga                                                                                                                             <  1:424626/251‑1 (MQ=255)
                                                                                                             cGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGAGGCGTCTGGCGCCGGCGGCCACCCCGGCGCGGTCGGGAATC‑GGCC‑ACCGCGCCGCGTGCCGCCCTTTTCACC‑AAAAGGAAAGCTGTAACAACTGTTCGAATCCGCTGCGTCCGGGAAAAGGTCAAATcctcgtgggggggagg                                                                                                                           >  2:357861/1‑236 (MQ=255)
                                                                                                                   tAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGGCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTCCGCCGAATGCCGCCCTTT‑CACCCCACAGGAAGTCTGTAACAACTGTTCGAATCCGCGTCGTCAGGAAAACGGTAAAAACTGCGTGGTGGAGAGTCcggcg                                                                                                                     >  1:472268/1‑251 (MQ=255)
                                                                                                                      cacaGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGAc                                                                                                                  <  1:177611/251‑1 (MQ=255)
                                                                                                                            ttATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGAGGTGTCTGGCGCAGGGGCTCCCCCCGGCGATGTGGGGAAAC‑GGCC‑ACCGCGCCGCTTGCCGCCCCTT‑CACCCAACAGGAAGCCTGTTACATCTGTTCGAATCCCCGGCGTCAGGGAAACGGGGAAATCTTCGGGGTGGAGAGGGcggcgggcagctcc                                                                                                            >  2:237460/1‑243 (MQ=255)
                                                                                                                                        tGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGAAGGGTCGGGGGCAGGGGGCCACCCGGGCGATGTCGGAAATA‑GGCC‑CCAGCGGCGACTGCCGCAACTT‑CACCGAAACGGAAGTCTGTAAAATTTGGTGGAATCCGCGGCGGCCGGAAAAAGGTCAAATCTGCGTGGGGGGGAGGCCGGCGGGCAATTTCGCCATTGAga                                                                                                  >  1:224159/1‑248 (MQ=255)
                                                                                                                                         gcgcTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAg                                                                                                <  1:473834/250‑1 (MQ=255)
                                                                                                                                                  gCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCTGTCACGCTGTTTCGTCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCGTTAACATCGGTGCGTATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGGGGTGGAGGGTGCGGCGTACATCTACGCCATTGAGCGGACGGGGc                                                                                         <  2:322202/248‑1 (MQ=255)
                                                                                                                                                   ccgcgCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTGCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCGCACCCGGGCGATGTCGGAAATC‑GGCC‑ACGGCGCCGATTGCCGCACTGT‑CACCGAACAGGAAGTCTGTAACATCTTTTCGAATCCGCGTCGTCAGTAAAACGGTCAAATCTGCGTTGTGGAGAGTCCGGCGGTCATCTACGCCATTGAGCAGACGGTGCAGtt                                                                                     <  2:197806/247‑1 (MQ=255)
                                                                                                                                                       gCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCa                                                                                 <  2:14952/251‑1 (MQ=255)
                                                                                                                                                                  aaGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGACTGGCGCCGGCACTCACCCCGGCGATGTTGGAAAAC‑GGCC‑CCCGCGCCGGTTGCCGCACATT‑CCCCGAAAAAGAAGTCTGGAAAATCTTTTCGGATCCGCGTCGCCAGGAAAACCGTCAAATCTGCGAGGTGGAAAGACCGGGGGGCATCTACCCCATTGGGCAGACGGGGcgggcttccaggcgttattt                                                                      >  1:157744/1‑232 (MQ=255)
                                                                                                                                                                    gTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTtg                                                                    <  1:26355/251‑1 (MQ=255)
                                                                                                                                                                    gTCGGCGCAGCGTATGGCGGTCACGCTGCGTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGAGGGCGGAAATC‑GGCC‑ATTGCGCCGATTGCCGCACGTT‑CACCGAACAGGAAGGCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTGAAATCTGCGTGGTGGAGAGTCCTGCGGACGTGTACGCCATTGAGCAGACTGGGCAGTTTGCAGGTCGATATTTtg                                                                    <  2:20650/251‑1 (MQ=255)
                                                                                                                                                                           cAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCCGGCGCTCCCCCCGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCCATTGCCGCACTTT‑CACCGCACAGGGAGTCTGTAACATCTGTTCGAATCCCCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGGGAGTCCGGCGGGCATCTACGCCCTTTAGCAGACGGGGGCGTTTTCAGGTCGTTTTTTTTTGtttgt                                                             >  1:40712/1‑248 (MQ=255)
                                                                                                                                                                                            gctgctTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGGGCAGGCGCTCACCCCGGCGATGGCGGGAATC‑GGGC‑ACCGCGCCGATTGCCGCACCTT‑CTCCCAACAGGGAGTCTGTAACAACTGTTCGAATTCCCGCCGTCCGGAAAACGGGCCAATCTGCGTGGTGGGGAGGTCGGGGGGCATCTACGCCCTTGAGCAGGCGGGGGAGGTTTCTGGGCGTTTTTTTTTGTTGGTgggggctcttgcaccgc                                            >  2:363197/1‑238 (MQ=255)
                                                                                                                                                                                                ctTCAGCGCGATCGTAGCGGCGGGATGCGGCTGGCGCAGGCGCTCCCCCCGGCGATGTCGGGAATC‑GGCC‑ACTGCGCCGATTGCCGCCCTTT‑CCACGAACAGGGAGTCTGTAACAACTGTTCGAATCCGCGGCGGCCGGGAAAAGGGCAAAACTTCGTGGTGGGGAGTCCGGCGGACAACTACGCCATTGAGCAGACGGGGGAGGTTTCCGGGCGTTTTTTTGTGTTGATGGGGGATCTGTCACCCCTggg                                        >  1:196852/1‑250 (MQ=255)
                                                                                                                                                                                                    aGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTGATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGc                                    <  1:174546/251‑1 (MQ=255)
                                                                                                                                                                                                       gcgATCGTAGCGGCGGTATTCGTCTGGCGCAGTCGCGCACCCGGGCGATGTCGGAAATC‑GGCC‑ACGGCGCCGATTGCCGCACTTT‑CACCTAACAGTAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGATGCCTGCGGACAGCGACGCCATTGAGCATACGGGGCAGGTGTCAGGTCTTGATGCTGGGTTTATTGGGCATCTGTCACCGCTGGAAGGCATc                                 <  2:438124/251‑1 (MQ=255)
                                                                                                                                                                                                       gcgAGCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATc                                 <  1:424036/251‑1 (MQ=255)
                                                                                                                                                                                                             gTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGGc                                                                                         <  1:89719/189‑1 (MQ=255)
                                                                                                                                                                                                                  ggcggGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGgatg                       >  1:251416/1‑250 (MQ=255)
                                                                                                                                                                                                                     gggATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGatat                   >  2:113740/1‑251 (MQ=255)
                                                                                                                                                                                                                               tGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGgag                                                                                                                              >  2:18838/1‑134 (MQ=255)
                                                                                                                                                                                                                               tGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGgag                                                                                                                              <  1:18838/134‑1 (MQ=255)
                                                                                                                                                                                                                                ggCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTg           <  1:28671/248‑1 (MQ=255)
                                                                                                                                                                                                                                ggCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGa          <  1:386890/249‑1 (MQ=255)
                                                                                                                                                                                                                                      ggCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTgg  <  1:153702/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
CGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATC‑GGCC‑ACTGCGCCGATTGCCGCACTTT‑CACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGG  >  NC_000913/494309‑494789

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