Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913731,2530TC22.4% 104.1 / 22.1 51F557F (TTT→TTCrhsCrhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base T (23/15);  new base C (5/6);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACG  >  NC_000913/731003‑731494
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 
tAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGCCGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGCTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGtttt                                                                                                                                                                                                                                                   >  2:23616/1‑251 (MQ=255)
      cAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCa                                                                                                                                                                                                                                              <  1:251919/250‑1 (MQ=255)
          cACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCACCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAAAGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCAGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTg                                                                                                                                                                                                                                          >  1:203406/1‑250 (MQ=255)
                             gAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCATCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAggg                                                                                                                                                                                                                      >  2:220174/1‑250 (MQ=255)
                                 tACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGACACTGCCACTGACGGCGATATCACCCGCTACCGATATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCCCCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAggg                                                                                                                                                                                                                  >  1:23536/1‑251 (MQ=255)
                                  aCGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGc                                                                                                                                                                                                                 <  1:74733/251‑1 (MQ=255)
                                                    aTGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCAt                                                                                                                                                                                               <  1:121491/251‑1 (MQ=255)
                                                           cGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAg                                                                                                                                                                                        <  1:317627/251‑1 (MQ=40)
                                                            gTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGc                                                                                                                                                                                       >  1:799/1‑251 (MQ=40)
                                                                                  tGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTg                                                                                                                                                                  >  2:110742/1‑250 (MQ=40)
                                                                                     cGGCGATATTACCCGCTACCGTTATGATACTCAACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCCTGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGaaac                                                                                                                                                              >  1:245483/1‑250 (MQ=255)
                                                                                     cGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGaa                                                                                                                                                                >  1:33493/1‑249 (MQ=40)
                                                                                                      aCCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGaa                                                                                                                                             <  1:67308/251‑1 (MQ=40)
                                                                                                                 aaTCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATg                                                                                                                                   <  2:78978/250‑1 (MQ=40)
                                                                                                                     cacacaGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATAc                                                                                                                              >  2:282414/1‑251 (MQ=40)
                                                                                                                           gTGACTTTCCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGTCTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAAcgcc                                                                                                                        <  2:239756/251‑1 (MQ=255)
                                                                                                                            ggACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAAcgccg                                                                                                                       <  2:260172/250‑1 (MQ=255)
                                                                                                                              aCTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACgccgcc                                                                                                                     >  1:12892/1‑251 (MQ=255)
                                                                                                                                  aCCCTGCACACCGGAAGCTGCCACCGGCAGCCGGAAAACCATGACGTGGAGACGTTACGGTCAGTTGCTGCGCTTCACTGACTGTTCCGGATATGTAACCCGCTATGACCATGACCGTTTTGGTCCGGTGACGGCGGTGCACCGCGAAGAAGGGCTGAGTCAGTACCGCGCATACGAAAGCCGTGGACCGATAATTGCCGTGAAAGAAACGCAGGGCCACGAAAAGCGGTATGCATACAACCCCGCCgggg                                                                                                                 >  2:211803/1‑249 (MQ=14)
                                                                                                                                          cAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGtgacctgacc                                                                                                         >  1:230244/1‑251 (MQ=255)
                                                                                                                                                        cccGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGcccc                                                                                           <  1:253590/250‑1 (MQ=255)
                                                                                                                                                        aCCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGcccc                                                                                           >  1:266442/1‑251 (MQ=255)
                                                                                                                                                        aCCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGcccc                                                                                           >  1:118229/1‑251 (MQ=255)
                                                                                                                                                                   gAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGgcagc                                                                                 >  1:90463/1‑250 (MQ=255)
                                                                                                                                                                                       gTTACGGTCAGTTGCTGAGCTTCACTGCCTGTTCCGGTTATGTACCCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGAt                                                            <  2:164876/251‑1 (MQ=255)
                                                                                                                                                                                                aGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGa                                                   <  1:145729/251‑1 (MQ=255)
                                                                                                                                                                                                     cTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGGTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGc                                              >  2:232224/1‑251 (MQ=255)
                                                                                                                                                                                                     cTGAGCTTCACTGACTGTTCCGGCTATGTAACCCGCTATGACCATGACCGTTTTGCTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGCCAGCCGTGGACAGGTAATTGCCGTGAAAGACACGACGAGCCATGAAACGCGGGATGAATACAACGACGCCGGAGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACCAAGTAAGCTGCGTGGGGAAAAGc                                              >  2:245149/1‑251 (MQ=14)
                                                                                                                                                                                                         gCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACCCAGTACGATGCGTGGGGCAAAGCCATc                                          >  2:198524/1‑251 (MQ=255)
                                                                                                                                                                                                              accgactgTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGCCAGATGCCGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCCTGAAACGCGGTATGAATACAACGCCGCCGGCGACCGGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACTATGCGTGGGGAAAAGCCATCTGTAc                                     >  2:60171/4‑251 (MQ=11)
                                                                                                                                                                                                                 gactgTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCAc                                  <  1:318974/251‑1 (MQ=9)
                                                                                                                                                                                                                 gactgTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCAc                                  >  1:316907/1‑251 (MQ=9)
                                                                                                                                                                                                                 gactgTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCAc                                  <  1:143978/251‑1 (MQ=255)
                                                                                                                                                                                                                    tgTTCAGGTTATGTAACCCGCTATGACCATGACCGCTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACCGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCCTTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGACAAGCCATCTGTACCACGcc                               >  2:40825/1‑250 (MQ=255)
                                                                                                                                                                                                                     gTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCa                               <  1:84305/250‑1 (MQ=9)
                                                                                                                                                                                                                     gTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCCCCGCGAGGCAGGGATGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGACACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGc                                >  1:132489/1‑249 (MQ=255)
                                                                                                                                                                                                                         cGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCg                          <  1:87751/251‑1 (MQ=9)
                                                                                                                                                                                                                             tATGTAACCCGCTATGACCATGACCGTTTGGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTTATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGGAGGGCGGTc                       <  2:158282/250‑1 (MQ=255)
                                                                                                                                                                                                                              aTGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTg                     <  1:171425/251‑1 (MQ=255)
                                                                                                                                                                                                                                gTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGATGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGAc                   <  1:30438/251‑1 (MQ=255)
                                                                                                                                                                                                                                 tAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGAGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAAAGGGACACAGTACGATGCGGGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGAcg                  >  2:15022/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aaCCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGAcgc                 >  1:256450/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cccGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACgcgc               <  1:220174/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cgttATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAg             <  1:312072/248‑1 (MQ=255)
                                                                                                                                                                                                                                       gttATGACCATGACCGCGTCGGGCAGAGGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGGGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGt            <  2:122743/249‑1 (MQ=11)
                                                                                                                                                                                                                                        cTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGt            >  1:95837/1‑250 (MQ=255)
                                                                                                                                                                                                                                           tgaccatgaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATgg        >  1:283568/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gaccatgaccGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGa       <  1:155144/251‑1 (MQ=255)
                                                                                                                                                                                                                                             accatgaccGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGTCTGAAGCGCAGTATGGa       >  2:97816/1‑250 (MQ=255)
                                                                                                                                                                                                                                              ccatgaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGCAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGaa      >  2:83326/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 tgaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGCCGCGCAGTATGGAATACg  >  2:26903/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 tgaccGCTTAGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACg  >  1:60089/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 
TAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACG  >  NC_000913/731003‑731494

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