Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913731,2600GA24.0% 96.1 / 24.8 50V560M (GTG→ATG) rhsCrhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base G (23/15);  new base A (6/6);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 7.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACG  >  NC_000913/731013‑731494
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          
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:316907/1‑251 (MQ=9)
                                                                                                                                                                                                       gactgTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCAc                                  <  1:318974/251‑1 (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)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          
CACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACG  >  NC_000913/731013‑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: