Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913731,2500TC25.5% 94.0 / 24.7 51R556R (CGT→CGCrhsCrhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base T (22/16);  new base C (7/6);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACG  >  NC_000913/731002‑731494
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
ttAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGttt                                                                                                                                                                                                                                                    >  2:145729/1‑251 (MQ=40)
 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)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
TTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACG  >  NC_000913/731002‑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: