Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 250,472 A→G intergenic (+2338/‑426) rayT → / → dinB REP‑associated tyrosine transposase/DNA polymerase IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913250,4720AG86.2% 144.6 / 14.8 65intergenic (+2338/‑426)rayT/dinBREP‑associated tyrosine transposase/DNA polymerase IV
Reads supporting (aligned to +/- strand):  ref base A (4/5);  new base G (29/27);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

GTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAGCGCAATGGCGGACCAGATGCTGCTGGATTCCAAAAATCCGCAAAGCGCGGGCAACCGGCGCATTGAGATTATGGTGCTGACC  >  NC_000913/250226‑250713
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
gTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGActcg                                                                                                                                                                                                                                               <  2:166489/251‑1 (MQ=255)
     gtgAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGctcgctcg                                                                                                                                                                                                                                           >  2:196166/1‑250 (MQ=255)
         aTGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGataat                                                                                                                                                                                                                                      >  1:33095/1‑251 (MQ=255)
           gagaAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGAtaataa                                                                                                                                                                                                                                    <  1:510220/251‑1 (MQ=255)
                 aCGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAattat                                                                                                                                                                                                                              >  1:328184/1‑251 (MQ=255)
                      gctgcaGAAGAGGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCACCACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAattattatt                                                                                                                                                                                                                          >  2:41329/1‑250 (MQ=255)
                         gctgAAGAAGAAATCCGCCCCCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTACTATTAc                                                                                                                                                                                                                        <  1:495808/249‑1 (MQ=255)
                              agaagaaATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCAt                                                                                                                                                                                                                 <  2:223358/251‑1 (MQ=255)
                                   aaaTCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACcg                                                                                                                                                                                                             >  1:63038/1‑250 (MQ=255)
                                    aaTCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACcgat                                                                                                                                                                                                           <  1:52752/251‑1 (MQ=255)
                                     aTCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACcgatg                                                                                                                                                                                                          >  2:502282/1‑251 (MQ=255)
                                         gCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACcgatgcgat                                                                                                                                                                                                      <  1:467034/251‑1 (MQ=255)
                                          ccACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACcgatgcgctg                                                                                                                                                                                                     >  2:436205/1‑248 (MQ=255)
                                           cACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCACTATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTCTTATTACCGGGCCTACCGATGCGATgg                                                                                                                                                                                                    >  2:327994/1‑251 (MQ=255)
                                             ccGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAGTTATTATTACCGGGCCTACCGATGCGATGGcc                                                                                                                                                                                                  >  2:173714/1‑251 (MQ=255)
                                                    tGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACaaa                                                                                                                                                                                            <  1:586073/250‑1 (MQ=255)
                                                     ggTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACaaaaa                                                                                                                                                                                          <  1:398251/251‑1 (MQ=255)
                                                     ggTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACaaaa                                                                                                                                                                                           >  1:98760/1‑250 (MQ=255)
                                                             ggccaCCAGTATCAACACCATTGCCATAGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCACAGGGATTACGCGTGCTGATCAAAGACGACCAGAACCGCAATATGGTTGCACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGGTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATTGCCTACAAAAACAATATCt                                                                                                                                                                                  <  2:75739/247‑1 (MQ=255)
                                                             ggCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATGATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAatat                                                                                                                                                                                    <  1:232523/249‑1 (MQ=255)
                                                                aaCCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTaca                                                                                                                                                                               <  2:177975/251‑1 (MQ=255)
                                                                 aCCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCATGGAGTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTGAAAACGCTGCTGGTGGAGCTCTCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTacaa                                                                                                                                                                              <  2:382435/251‑1 (MQ=255)
                                                                 aCCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTaca                                                                                                                                                                               >  2:200313/1‑250 (MQ=255)
                                                                    aGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCACTATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTacaaca                                                                                                                                                                            >  2:357100/1‑250 (MQ=255)
                                                                            cacCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAAcc                                                                                                                                                                   >  1:135113/1‑251 (MQ=255)
                                                                                   gCCCGTGATCCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCggg                                                                                                                                                            >  1:169764/1‑251 (MQ=255)
                                                                                                  atGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGGTGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTc                                                                                                                                             <  2:187913/251‑1 (MQ=255)
                                                                                                     gAAGCGAATCTGGAGCTGGAGATTGATCCCCAGGGATTACGCGTGCTGATTACAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATACAATTATTATTACCGGGCATCCCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGGGCCCGCGCACTTTCgac                                                                                                                                          >  2:305537/1‑249 (MQ=255)
                                                                                                         cGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCtcgt                                                                                                                                      >  1:563793/1‑251 (MQ=255)
                                                                                                                    atggagatTGTTTCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGTCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGa                                                                                                                           <  1:356741/251‑1 (MQ=255)
                                                                                                                        agatTGTGCCCCAGGGATTACGCGGGCTGAGTAAAGACGGCCAGAACCGCAATATGGTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCGGGTGGAGCTCGCGCGAGTGTTCGGCTCGCTCGGTAATAAAATTATTAGTACCGGGCATACCGATGCGATGGGCTACAAAAACAAGATCTACAGCAACTGGAACCTTGCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAagag                                                                                                                       <  2:212911/251‑1 (MQ=255)
                                                                                                                           ttGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGcc                                                                                                                    >  2:550975/1‑251 (MQ=255)
                                                                                                                                             gcgTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTGCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGTTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATaaa                                                                                                  <  1:327994/251‑1 (MQ=255)
                                                                                                                                               gTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGt                                                                                                <  2:346831/251‑1 (MQ=255)
                                                                                                                                                 gCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGa                                                                                              <  2:264214/251‑1 (MQ=255)
                                                                                                                                                       tAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCGGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTGCGGGTGACCGCGCGCGTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAg                                                                                         <  2:619714/250‑1 (MQ=255)
                                                                                                                                                          aGACGACCAGAACCGCAATATGTTTGACCGCGGCCGCGCGAAGATTCTACCGTTCTTTAAAAAGCTGCTGGTGGAGCTCGCGACAGTGTTCGACTCGCTCGATAATAAAATTATTATTACCGAGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCTCGCTTTCGGCTAGTCGGGTGCTGGAAGAGGACGGAATGCCGGGCGATAAAGGGATGCAGGTa                                                                                      >  2:416725/1‑250 (MQ=255)
                                                                                                                                                          aGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATg                                                                                                              <  1:550430/226‑1 (MQ=255)
                                                                                                                                                          aGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGGTCTGGAAGAGGCCGGAATg                                                                                                              >  2:550430/1‑226 (MQ=255)
                                                                                                                                                            acgacCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAgc                                                                                   <  1:583665/251‑1 (MQ=255)
                                                                                                                                                                     aCCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAGCGCAATGGCg                                                                          >  2:137979/1‑251 (MQ=255)
                                                                                                                                                                                           cAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGCAGGGGCCGGAATGCCGGAAGATAAAGGGATGAAGGGAAGCGCAATGGCGGACCAGATGAGGCTGGATTcc                                                     >  2:15948/1‑250 (MQ=255)
                                                                                                                                                                                            aGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAGCGCAATGGCGGACCAGATGCTGCTGGATTCCa                                                    >  1:355106/1‑250 (MQ=255)
                                                                                                                                                                                              cgcgAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAGCGCAATGGCGGACCAGATGCTGCTGGATTCCaaa                                                  >  1:261314/1‑250 (MQ=255)
                                                                                                                                                                                               gcgAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCATTGGTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAGCGCAATGGCGGACCAGATGCTGCTGGATTCCaaaaa                                                <  2:355106/251‑1 (MQ=255)
                                                                                                                                                                                                 gAAGATTCTGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATACTCAAATTCTTATTACCGGGCATCCCGATGCGATGGCCTACAAAAACCATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTc                                                                                                                                             >  2:375973/1‑156 (MQ=255)
                                                                                                                                                                                                 gAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAGCGCAATGGCGGACCAGATGCTGCTGGATTCCAAAAATc                                              >  2:234353/1‑251 (MQ=255)
                                                                                                                                                                                                 gAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAAGAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTc                                                                                                                                             <  1:375973/156‑1 (MQ=255)
                                                                                                                                                                                                  aaGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGACGGAAGATAAAGTGAAGCAGGTAAGCGCAATGGCGGGCCAGATGCTGCTGGATTCCAAAAAttc                                             >  2:475418/1‑249 (MQ=255)
                                                                                                                                                                                                        aTGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAGCGCAATGGCGGACCAGATGCTGCTGGATTCCAAAAATCCGCaaa                                        <  2:507233/250‑1 (MQ=255)
                                                                                                                                                                                                              ttCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAGCGCCATGGCGGACCCGAAGCTGCTGGATTCCAAAAATCCGCAAAGCGCggg                                 >  2:89987/1‑251 (MQ=255)
                                                                                                                                                                                                                  ttAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAGCGCAATGGCGGACCAGATGCTGCTGGATTCCAAAAATCCGCAAAGCGCGGGc                                <  2:98760/248‑1 (MQ=255)
                                                                                                                                                                                                                   tAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATCATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGATGACCGCGCGCTTTCAGCTCGTCGGGTGTTGGAAGAGGCCGGAATGCCGGAAGATAAAGGGATGCAGGTAAGCGCAATGGCGGACCAGATGCTGCTGGATTCCAAAAATACGCAAAGCGCGGGCAAcc                            >  1:74192/1‑251 (MQ=255)
                                                                                                                                                                                                                      aaCGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGGAGATAAAGGGATGCAGGTAAgc                                                                                   >  2:501218/1‑193 (MQ=255)
                                                                                                                                                                                                                      aaCGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAgc                                                                                   <  1:501218/193‑1 (MQ=255)
                                                                                                                                                                                                                      aaCGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAGCGCAATGGCGGACCAGATGCTGCTGGATTCCAAAAATCCGCAAAGCGCGGGCAACCGgc                         <  2:162877/251‑1 (MQ=255)
                                                                                                                                                                                                                       aCGCTGCTGGTGGAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAGCGCAATGGCGGACCAGATGCTGCTGGATTCCAAAAATCCGCAAAGCGCGGGCAACCGgcg                        >  1:358681/1‑251 (MQ=255)
                                                                                                                                                                                                                              tggtggAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTGCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTaa                                                                                     <  2:523574/183‑1 (MQ=255)
                                                                                                                                                                                                                                gtggAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGGAGATAAAGGGATGCAGGGAAGCGCAATGGCGGAcccga                                                                   >  2:205728/1‑196 (MQ=255)
                                                                                                                                                                                                                                gtggAGCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAGCGCAATGGCGGACCAGa                                                                   <  1:205728/199‑1 (MQ=255)
                                                                                                                                                                                                                                gtggAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAGCGCAATGGCGGACCAGATGCTGCTGGATTCCAAAAATCCGCAAAGCGCGGGCAACCGGCGCATTGAGAt               <  2:261314/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gCTCGCGCCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAGCGCAATGGCGGACCAGATGCTGCTGGATTCCAAAAATCCGCAAAGCGCGGGCAACCGGCGCATTGAGATTATgg          <  1:436205/251‑1 (MQ=255)
                                                                                                                                                                                                                                           gcCAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATAACGATGCGCTGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGCAGATAAAGTGATGCAGGTAAGCGCAATGGCGGACCAGCTGCTGCTGGATTCCCAAAATCCGCAAAGCGCGGGCAACCGGCGCATTGGGATTATGGTGCTgc    >  2:82368/1‑250 (MQ=255)
                                                                                                                                                                                                                                           gcCAGTGGTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAGCGCAATGGCGGACCAGATGCTGCTGGATTCCAAAAATCCGCAAAGCGCGGGCTACCGGCGCAGTGAGATTATGGTGCTGa    <  2:421901/251‑1 (MQ=255)
                                                                                                                                                                                                                                             cAGTGTTCGGCTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAACACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAGCGCAATGGCGGACCAGATGCTGCTGGATTCCAAAAATCCGCAAAGCGCGGGCAACCGGCGCATTGGGATTATGGTGCTGAcc  >  1:33191/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
GTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAGCGCAATGGCGGACCAGATGCTGCTGGATTCCAAAAATCCGCAAAGCGCGGGCAACCGGCGCATTGAGATTATGGTGCTGACC  >  NC_000913/250226‑250713

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