Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 248,203 G→A intergenic (+69/‑2695) 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_000913248,2030GA96.0% 213.0 / NA 76intergenic (+69/‑2695)rayT/dinBREP‑associated tyrosine transposase/DNA polymerase IV
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (41/32);  minor base C (3/0);  total (44/32)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGTTGATTATATTTATATAAATCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCAT  >  NC_000913/247970‑248449
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      
catGTTGATTATATTTATATAAATCCAGTAAAGCAGGTTTGGTTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACTGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGTATAAGGCGttt                                                                                                                                                                                                                                        <  1:373322/250‑1 (MQ=255)
   gTTGCTTATATTTATATAAATCCAGTAAAGCATGGTTGGGTAAAGCCAGTGAGGGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATAAATTGGGCGGGGGACGTAACGGATTTTAGAGCCGGAGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCAAGCGCGTCTTATCATTCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTAttc                                                                                                                                                                                                                                    >  1:283161/1‑249 (MQ=255)
   gTTGATTATATTTATATAAATCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGTTAAGCCGGATAAGGCGTTTATGc                                                                                                                                                                                                                                    >  1:50714/1‑251 (MQ=255)
    ttGATAATATTTATATAACTCCGGTACAGCATGCTTGGGTATAGCACGTGTGTGATTGGCCATTCTAAACGTTCCCTCGCGATATAGCGCGAGGGTTGTATCCCATCGATTGGGAGGGGGACGTAACGGATTTTAGTGCAGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCATGATGCGACGCTAGCGAGTCTTCTCATGCCTACCAACTTGTGCCGGGTCGGTAAGCCGGATAAGGCGTTTATGcc                                                                                                                                                                                                                                   >  2:455341/1‑251 (MQ=255)
            atttatatAAATCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATAGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCgg                                                                                                                                                                                                                           >  1:297118/1‑251 (MQ=255)
            aattatatAAATCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTAGCCATTCTCAACGTTCCATCGTGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCgg                                                                                                                                                                                                                           >  2:52411/3‑251 (MQ=255)
                  ataAATCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGtt                                                                                                                                                                                                                      <  2:445976/250‑1 (MQ=255)
                    aaaTCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATg                                                                                                                                                                                                                   <  1:153591/251‑1 (MQ=255)
                      aTCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATgcg                                                                                                                                                                                                                 <  2:156193/251‑1 (MQ=255)
                       tCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGtt                                                                                                                                                                                                                                         >  1:34335/1‑226 (MQ=255)
                       tCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGtt                                                                                                                                                                                                                                         <  2:34335/226‑1 (MQ=255)
                        ccAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGTTCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATgcgc                                                                                                                                                                                                                >  2:106950/1‑250 (MQ=255)
                             acagcagggaggagTAAAGAAAGTGAGGGATTGGCCATTCTAAACGTTCCATCGCTATGTCGAGCGAGGGTTGTAGCCCAACGGTTGGGCGGGGGACGTAACGGAATTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATg                                                                                                                                                                                                           <  2:128374/237‑1 (MQ=37)
                             aaaGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGc                                                                                                                                                                                                          >  1:337977/1‑251 (MQ=255)
                              aaGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGcc                                                                                                                                                                                                         >  1:477670/1‑251 (MQ=255)
                                gCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGCTTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACACACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTg                                                                                                                                                                                                       >  2:69241/1‑251 (MQ=255)
                                    ggTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGc                                                                                                                                                                                                   >  1:430792/1‑251 (MQ=255)
                                     gTTGGGTAAAGCAAGTGCGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGTTTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCCGTTATGCGCAGATGCCTGATGCg                                                                                                                                                                                                  >  1:418976/1‑251 (MQ=255)
                                           tAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTg                                                                                                                                                                                            >  1:174970/1‑251 (MQ=255)
                                              aGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGgcg                                                                                                                                                                                         >  1:306243/1‑251 (MQ=255)
                                                     gagtgaTTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGCTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCtt                                                                                                                                                                                   >  1:208322/1‑250 (MQ=255)
                                                       gtgaTTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATc                                                                                                                                                                                <  1:388856/251‑1 (MQ=255)
                                                        tgaTTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATc                                                                                                                                                                                >  2:16884/1‑250 (MQ=255)
                                                         gaTTGGCCATTCTCAACGTTCCATCGAGATGTCGCGCGAGGGTTATATACCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGTCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCAt                                                                                                                                                                              <  2:299743/251‑1 (MQ=255)
                                                           ttGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGc                                                                                                                                                                            >  2:186541/1‑251 (MQ=255)
                                                                    ctcAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACt                                                                                                                                                                   >  1:418008/1‑251 (MQ=255)
                                                                      cAACGTTCCATCGCGATGTCGGGTGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTGAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCGGATGCGTCGCGAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATTGGTAAGCTGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACt                                                                                                                                                                   <  1:74612/249‑1 (MQ=38)
                                                                       aaCGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTtgt                                                                                                                                                                >  1:4425/1‑251 (MQ=255)
                                                                          gTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGcc                                                                                                                                                             >  2:248893/1‑251 (MQ=255)
                                                                                  gcgATGTCTCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGGTTATGCCGCATCCGTCATGTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTAcgagcccgcgaatatttgcg                                                                                                                                                     <  1:455341/251‑20 (MQ=37)
                                                                                  gcgATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGt                                                                                                                                                     >  1:185400/1‑251 (MQ=255)
                                                                                     aTGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTACGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAgg                                                                                                                                                  >  2:103364/1‑251 (MQ=255)
                                                                                              gAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGACCGGTAGGACGGATa                                                                                                                                           >  2:374180/1‑249 (MQ=255)
                                                                                                gggTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGc                                                                                                                                       <  1:346202/251‑1 (MQ=255)
                                                                                                      tatCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGATAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGAGTTCAc                                                                                                                                 >  1:225514/1‑251 (MQ=255)
                                                                                                      tatCCCATCGATTGGGCGGGGGACGTAACGGATCTTAGATCCGGGGCGCGAATCATTTCATAATTGAGCGCAGCTGCCTGATGCGAAGCTAGCGCGTCATATCATGCCTCCAAACTTGTGCCGCATCGGTACGCCGGATAAGGAGTTTATGCCGCATCCGGCACTTATGAGCAGATGCCTGCTGCGTCGCTGGCGCGTCTTATCATCCCTACAAACTTCTGCCTGAACGTTCGGACGCATAAGGCGTTCAc                                                                                                                                 >  2:455302/1‑251 (MQ=255)
                                                                                                            aTCGATTGGGCGGGGGTCGTAACGGATTTTAGGGGCGGGGGGCTTATCAGTTTATAATTGTGCTAAGATGCGTGATGCGACGCTTGCGCGTCTTATCTTGCATGCAAAGTTGTGCCGGAGCGGTAAGTCGGATAAGGCGTTTATGCCGCATCCGGCAGTGATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATTATGCCTACAAACTGGTGCCTGAGCGGTAGGACGGATAAGGCGTGCACGCCGCa                                                                                                                           <  2:456012/251‑1 (MQ=255)
                                                                                                              cGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATc                                                                                                                         >  1:155783/1‑251 (MQ=255)
                                                                                                              cGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGAGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCCCGCCGCATc                                                                                                                         >  2:416778/1‑251 (MQ=255)
                                                                                                                            aCGTAATGGATTTTAGTTCCGGTTCGCCTATCATTTCCTAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCTGTAAGCCGTATAAGGCGTTTATGCCGCATCCTGCAGTTATGCGCAGATGTCTGATGCGACGCGGGCGCGTCTTATCATGCCTACAAACTTTTGCCTGAACGGTATGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCa                                                                                                            <  2:283880/250‑1 (MQ=25)
                                                                                                                             cGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCAcc                                                                                                          <  1:248893/251‑1 (MQ=255)
                                                                                                                                 aCGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGAAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGACTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTTGTGCACCAATg                                                                                                      >  1:197118/1‑251 (MQ=255)
                                                                                                                                     atgttAGTGCCGGGGAGCGTATTATTTCATAATTGTGGGAGGATGCCGGATGCGACGCTTGCGCGTCTTATCAGGCCTACAAACTTGTGCCGGATCGGTAATCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTgg                                                                                                  <  2:380232/248‑1 (MQ=37)
                                                                                                                                          agtggcgtggAGCGTATCATTTCATAATTGGGCGAAGAGGCTTGATGCGACGCTAGCGCGGCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGAGAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCTCCAATGCTGGATACt                                                                                             <  2:442903/243‑1 (MQ=37)
                                                                                                                                              ccGGGGAGCGTATCATTTCATAATTGTGCGCAGAGGCCTGAGGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACa                                                                                         <  1:416778/251‑1 (MQ=255)
                                                                                                                                               cGGGGCGCGAATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGGCACCAATGCTGGATACTTACAt                                                                                        >  2:29774/1‑251 (MQ=255)
                                                                                                                                                    aGCGTATCATTTCATACTTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACACACTTGTGCCGGCTCGGTAAGCCGGATCAGGCGTTTATGCCGCATCCGGCAGCTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTAAAAACTTCAGCCTGAACGGTCGGACAGATAAGGCGTTCACGCCGCATCCGGCAGTTGGGCACCACAGCTGGATACTTACATCAACg                                                                                   >  2:225374/1‑251 (MQ=255)
                                                                                                                                                         aTCATTTCATCACTGTGCGCAGATGCCTGATGCGACGGTAGAGTGTCTTATCATGCCTACAAACTTGTGCAGGCTCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCCGATGCCTCATGAGACGCTGGCGCGTCTTATCATGCCTCCAAACTTGTGCCAGAACGGTAGGACGGATAATGCGTTCACGCCGCATTCTGCAGTTCAGCCCCACTGCTGGATTCTTAAATCACCGccccc                                                                              >  2:479319/1‑251 (MQ=255)
                                                                                                                                                         aTCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGccccc                                                                              >  1:167820/1‑251 (MQ=255)
                                                                                                                                                         aTCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGcccc                                                                               <  1:227395/250‑1 (MQ=255)
                                                                                                                                                         aTCATTTCATAATTGTGCGCAGATGCATGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACCAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACCCCCCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGccccc                                                                              >  2:432060/1‑251 (MQ=255)
                                                                                                                                                           cATTTCATAATCGTGCTCAGATACCTGATGCGACTCTAGCGCGTCTTATCTTGCCTACAAACCTGGGCCGTATCTGTAAGCGTGATAAGGCTTTTTTGCCGCATCGGACAGTTATGCGCAGATGCCTGATGCTACGCTGGCGCGGCGTATCATGCCTACAAACTTGTGCCTGATCGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCGCCAATGCTTTATACTTACATCAACGcccc                                                                               <  2:295636/248‑1 (MQ=255)
                                                                                                                                                                 aTAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAGCGCCCCCATAATCtt                                                                      <  2:383844/251‑1 (MQ=255)
                                                                                                                                                                      tgtgCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCCCGTCTTATCATGCCTACCAACTTGTGCCTGAACGGTAGGACGGATAAGTCGTTCACGCCCCATCCGGCAGTTGGGCACCAATGCTGGGTCCTTACATCAAAGCCCCCATAATCTTCacc                                                                  >  2:302823/1‑248 (MQ=38)
                                                                                                                                                                       gtgCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGGGTCTTATCATGCCTACAAACTTGTGCTTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTc                                                                <  1:225374/251‑1 (MQ=255)
                                                                                                                                                                        tgCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTcc                                                               >  1:423733/1‑251 (MQ=255)
                                                                                                                                                                            cAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGAAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCt                                                           >  1:82724/1‑251 (MQ=255)
                                                                                                                                                                                 gCCTGATGCGACGCGAGCGCGTCTTATCATGCCTACAAACTTGTTCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATcc                                                      <  1:277963/251‑1 (MQ=255)
                                                                                                                                                                                     gATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGTTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCGGGCGCGTCTTATCATGCCTACAAACGTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAAGCTACAACTCCATCTTATCCgg                                                    <  2:208322/249‑1 (MQ=37)
                                                                                                                                                                                     gATGCGACGCTAGCGCGTCTTATCATGCCGACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCgg                                                    <  2:126269/249‑1 (MQ=255)
                                                                                                                                                                                          gACGCTAGCGCGTCTTATCATGCCTACACACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTCACATCAACGCCCCCATCATCTTCAACTCCAGCTCATCCGGCACTTCg                                             >  1:309632/1‑251 (MQ=255)
                                                                                                                                                                                          gACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCATACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCg                                             >  1:86762/1‑251 (MQ=255)
                                                                                                                                                                                                agctcgTCTTATCATGCCTACAAACTTGGGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCGGGCGCGTCTTAGCATGCCTACAAACTTGTGCCGGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACGTACAGCATCGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTtat                                         <  1:244923/245‑1 (MQ=38)
                                                                                                                                                                                                   gcgTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGAc                                    <  2:174970/251‑1 (MQ=255)
                                                                                                                                                                                                          aTCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACAt                             >  1:116494/1‑251 (MQ=255)
                                                                                                                                                                                                           tCATGCCTACAAACTTGTGCCGGATCGGTAAGCAGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATg                            >  1:4557/1‑251 (MQ=255)
                                                                                                                                                                                                            cATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGc                           <  1:211811/251‑1 (MQ=255)
                                                                                                                                                                                                            cATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGGGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATg                            <  1:484149/250‑1 (MQ=255)
                                                                                                                                                                                                                   aCAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCg                    <  1:29774/251‑1 (MQ=255)
                                                                                                                                                                                                                       ccTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGTAGTTATGCGCAGATGCCTGATGCGCCGCCGGCGCGTCTTATCATGCCTCCAAACTTGTGCCTTCACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCACTGCTGGATACTTACATCCCCCCCCCCATAATCTTCAACTCCAGCTCATCCGCCACTTCGTTATACGCCAGCACATACAGCCCCGgccc                >  2:246603/2‑249 (MQ=37)
                                                                                                                                                                                                                          tgtgCCGGATCGGTACGACGGATAAGGCGTTTCTGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACCAACTTGTGCCTGCACGGTAGGACGGCTAAGGCGTTCCAGCCGCATCAGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGAACGCACATGCAGCCCCGgcgcca              >  1:409582/1‑248 (MQ=38)
                                                                                                                                                                                                                                  aTCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGCTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGc                                                                                                     >  2:82462/1‑155 (MQ=255)
                                                                                                                                                                                                                                  aTCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGc                                                                                                     <  1:82462/155‑1 (MQ=255)
                                                                                                                                                                                                                                   gcGGTAAGCCGGATAAGGCTTTTATGCCGCATCCGGCAGTTAGGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCGGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCTCAAACAACCTTGc    <  1:170031/250‑1 (MQ=255)
                                                                                                                                                                                                                                   gcGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGc    <  2:273894/250‑1 (MQ=255)
                                                                                                                                                                                                                                     ggTAAGCCGGATAAGGCGATTATGCCGCCTCCGGCAGTTATGCGAAGATGCCTGGTGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCCGTTGTGAACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCCCTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTTCAt  >  2:188900/1‑251 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      
CATGTTGATTATATTTATATAAATCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCAT  >  NC_000913/247970‑248449

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: