Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 65,811 G→T intergenic (‑31/+44) polB ← / ← araD DNA polymerase II/L‑ribulose‑5‑phosphate 4‑epimerase AraD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091365,8110GT96.7% 184.1 / ‑0.4 60intergenic (‑31/+44)polB/araDDNA polymerase II/L‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (27/31);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.73e-01

TAACGCCAGCGGTGTCAGGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑AAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATC‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAAC  >  NC_000913/65571‑66055
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
tAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTG‑AAAACGTAACCCGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCCGTGTCGGGGTAAAA‑TAAAACCTGCCTGCGCCACGCGGTAATTTCCA‑CAAAAAACCAGGCTTTATTATAGcccgg                                                                                                                                                                                                                                               >  1:326304/1‑248 (MQ=255)
tAACGCCAGCGGTGTCAGGCGAACGCCTTGTTAACCCTGCAAAATATGCGGAGCGCGGGGAACCTGACCGGCGGGAATAAACGCCACGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTg                                                                                                                                                                                                                                                <  1:483051/250‑1 (MQ=255)
      cAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGACCGGCGGGAATAAACGCCACGGACTCTGGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGttt                                                                                                                                                                                                                                            <  1:655918/248‑1 (MQ=255)
                    gatcgCCTTGTGCACCCTGAAAAATTGGCTGGGCGCGGGGAACCGGATCGGCGGGAATAAACGCCACGGCCGCGTGTGGGG‑AAAGGTTAACCTGCAACGGCCCGTTGGCCGGCGCCAGCAAGAAGGAGATTTCTGGCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAATCCGGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTgg                                                                                                                                                                                                                           <  2:688177/247‑1 (MQ=255)
                    gaacgCCTTGTTCCCCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTgg                                                                                                                                                                                                                           <  1:601165/247‑1 (MQ=255)
                    gaacgCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTgg                                                                                                                                                                                                                           <  1:20548/247‑1 (MQ=255)
                              gtCACCCTGCAAAATATGCGGAGCGCGGGTAGCCTGATCGGCGGAAATAAACGCCACGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCTTTGTCCGTCGCCAGCCAGAAGGAGACTGCTGTCCCT‑GGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAg                                                                                                                                                                                                                 <  1:109550/250‑1 (MQ=255)
                                 ccccggcaaaatatgtggggggcgGGGAACCTGAGCGGGGGGAATAAACGCCACGGACTCTGGCGTGG‑CAAGCGTAACCGGAAACGGCCCGTTGTCCGTCGCCAGCCAGAGGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCa                                                                                                                                                                                                              <  2:677838/230‑1 (MQ=255)
                                   ccggaaaaaTATGCGGGGCGCGGGGAACGTGTCCGGGGGAAATAACCGCCCTGGCCTCTTCGGGGG‑CAAGCGTAACTTGCACCGGCCGGTGGTCCGCTGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGTTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCAtt                                                                                                                                                                                                            <  2:75618/246‑1 (MQ=255)
                                   agccgtgagctcttccgatcgggggAGCCTGATCGGCGGGAATAAACGCCACGGACGCTTG‑GGTGGCAAGCGAAGCATGCAACGGGCGGTTGTCCGTCGCCAGCAAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGCCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCAtt                                                                                                                                                                                                            <  1:74293/230‑1 (MQ=255)
                                              tGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGGTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCCGTGTCGGGTTTAGA‑TAAAAACTGCCTGCGCCCCCCTG‑GAATTTCATCAAAAAACCAAGCTTGATTTTATCCTGGGTTTGTTTGTTTTTCTGTGGGTTTTTACAGttcttattggccgta                                                                                                                                                                                                  >  1:521394/1‑236 (MQ=255)
                                               gCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTaa                                                                                                                                                                                                 <  1:64940/250‑1 (MQ=255)
                                                cTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGGG‑CAAACGTAACCTGCAACGGCCCGGTGGCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCC‑TGCGGGGTGGCCCGCCCGTGTCGGGTTAAGAGAAAAACCT‑CCTTCGCCACCCAG‑TAATTTCATCAAAAAAACAGGCTTTATTATAACCTGGGTTCGTTTGATTGGGTGGGGTTTTATACAggtattaattccccttat                                                                                                                                                                                                >  1:501655/1‑233 (MQ=255)
                                                                    gggcggGAATAAACGCCACGGACTCTGGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTAc                                                                                                                                                                           <  1:277487/250‑1 (MQ=255)
                                                                     ggcggAAATAAACGCCACGGACTCTGGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTGGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACg                                                                                                                                                                          <  1:323794/251‑1 (MQ=255)
                                                                              aaaCGCCACGGACTCTTGCGGTG‑CACGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGCGTCCCGCCAGTGACGGGGTAAGA‑TAAAACCTGCCTGCGCCCCGCTG‑TAATTCCATCAACAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGGGGTTTTATACATTCCTTACTGCCCGTAATAACACTTCGCGCCATGCTTACgacacaaggg                                                                                                                                                                 >  2:368701/1‑242 (MQ=255)
                                                                              aaaCGCCACGGACTCTGGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGGCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTta                                                                                                                                                                                                                      <  1:705372/198‑1 (MQ=255)
                                                                                ccGCCACGTACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCTT‑TGCGGGATGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGtt                                                                                                                                                              <  1:121297/250‑1 (MQ=255)
                                                                                  gCCACGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGtttt                                                                                                                                                            >  1:134433/1‑250 (MQ=255)
                                                                                     aCGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCtt                                                                                                                                                                                         >  1:631855/1‑220 (MQ=255)
                                                                                     aCGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTAttc                                                                                                                                                         >  1:473885/1‑249 (MQ=255)
                                                                                                      caaacGTAACCTGCAACGGCCCGTTGTCCGCCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑a                                                                                                                                       <  1:587924/247‑1 (MQ=255)
                                                                                                      cAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGGTTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCtt                                                                                                                                                                                         >  1:368622/1‑204 (MQ=255)
                                                                                                             aaCCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCgg                                                                                                                                <  1:89367/251‑1 (MQ=255)
                                                                                                                  gCAACGGCCCTTTGTCCTGCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTTTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTaa                                                                                                                             <  1:303213/249‑1 (MQ=255)
                                                                                                                    tcccgatcttggcccGTCCCCAGCCAGAAGGAGACTTC‑GGCCCTGTGCGGGTTGTCCCGCCATTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTg                                                                                                                          <  2:444271/237‑1 (MQ=255)
                                                                                                                         cccGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCCCGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATACCCTGGTTTCGTGTGATTGGCTGCGGTTTTATACAGCCATTACTGCCCCTAATATGCCTTTGCTCCCTGCTTCAGC‑AGATAGACTTTATT‑CAGCAGCGATCGCTGCCAT‑TCCGGCAACTGCGgcgc                                                                                                                    >  1:42895/1‑251 (MQ=255)
                                                                                                                             ttGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTg                                                                                                                          >  1:444271/1‑241 (MQ=255)
                                                                                                                               gTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTg                                                                                                              >  1:376790/1‑251 (MQ=255)
                                                                                                                                cccGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGa                                                                                                             <  1:533309/250‑1 (MQ=255)
                                                                                                                                                aGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGGTAAGA‑TAAAACCTGCCTGCGCCCCGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCCAGTTTCGTTTTATTTGCTGTGGTTTTATAAAGGCATTACTGACCGTAATATGCCTTCGCGCCATTCTTTCGCCAG‑TAGTTTTTTTTTCA‑CAGCTGTTGCTGC‑ATATCCGGGTACTTCGGCGCTAACTGACGGCAGaaattccccc                                                                                             >  2:458397/1‑243 (MQ=255)
                                                                                                                                                aGGAGACTTCTCTCCCT‑TGCGCGGTGTCCCCACAGCGTCGGGTTAAGC‑TAAAACCTGCCCTCGCCACGCAG‑CAATTCCCTCACCAAACCAGGCTTGATTATAGCCTGGTTTCGCTTGAATGGCTGTGGTTTTATAAAATAATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑ATATAGAGTTTATT‑CAGAAGCGTTTGGTGC‑ATCTCCGGTAACTGCGGCGCCAACAGAAGGCAGAatttcaccc                                                                                             >  2:653594/1‑244 (MQ=255)
                                                                                                                                                           gTCCCT‑TGGGGGGTGACCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGGGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGa                                                                                    <  1:400670/249‑1 (MQ=255)
                                                                                                                                                           gTCCCT‑GGGGGGGTGTCCCGCCAGTGGCGGGTTAAGA‑TAAAACCGGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGAcc                                                                                  <  2:114550/251‑1 (MQ=255)
                                                                                                                                                                    gggggTGCCCCGACAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAg                                                                          <  2:276447/250‑1 (MQ=255)
                                                                                                                                                                      gggTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTTTT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGc                                                                         >  1:476957/1‑250 (MQ=255)
                                                                                                                                                                         tgtCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACg                                                                      >  1:35960/1‑250 (MQ=255)
                                                                                                                                                                            cccGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑AAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATC‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATgg                                                                  <  1:39847/251‑1 (MQ=255)
                                                                                                                                                                            cccGCAAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTTATTATAGCCTGGTTGCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTGCCAGCACGATgg                                                                  <  1:572596/251‑1 (MQ=255)
                                                                                                                                                                             ccGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAatat                                                                                                                                                                                              <  1:309577/129‑1 (MQ=255)
                                                                                                                                                                             ccGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAatat                                                                                                                                                                                              >  2:309577/1‑129 (MQ=255)
                                                                                                                                                                                cacgtggagggTTACGA‑TAACACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCGTGATTATAGACTGTGTTCTTTTGATTGGCTGTGGTTTTATACAGGCATTACTGCCCGTAATATGCCGTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCAGATAAGCGTCCTCTTCCAGCGCGATGGCGtt                                                              <  2:119869/243‑1 (MQ=255)
                                                                                                                                                                                    tgtCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑AAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATC‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGc                                                          >  1:214703/1‑251 (MQ=255)
                                                                                                                                                                                          gTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGc                                                     >  1:707532/1‑250 (MQ=255)
                                                                                                                                                                                           ttAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCAt                                                   >  2:143621/1‑251 (MQ=255)
                                                                                                                                                                                             aaGA‑TAAACCCTTCCGGCGCCACGCTG‑TAATTCCATAAAAAAACCAGGCGTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCCCCGCATcg                                                 <  2:642351/251‑2 (MQ=255)
                                                                                                                                                                                               gA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCCAGCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTc                                               <  1:466218/251‑1 (MQ=255)
                                                                                                                                                                                                      aCCTGCCTGCGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACCGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCAtt                                         >  1:265468/1‑251 (MQ=255)
                                                                                                                                                                                                        ctgcctgcGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGCATATCCCCATATAAGCGACCTCTTCCAGCACACTGGCGTTATGCACCGCATCTTCGGCTtttt                                       >  2:147062/1‑251 (MQ=255)
                                                                                                                                                                                                         tgcctgcGCCACGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATCTAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCAttttt                                      >  2:291300/1‑251 (MQ=255)
                                                                                                                                                                                                                    cGCTG‑TAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCaa                           >  1:342363/1‑250 (MQ=255)
                                                                                                                                                                                                                      gctgtaattCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCaaat                         <  1:680083/243‑2 (MQ=255)
                                                                                                                                                                                                                         gtaattCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGACAGAGGTTTTATACATTCATTACTGCCCGTAATATGCCTTCGCTCCATGTTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCAGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATACGCGACCTCTTCCACCACGAAGGAGTTATGCACCTCATCTTCGGCCTTTTTGcccaagcaaatggg                       >  2:274872/6‑240 (MQ=255)
                                                                                                                                                                                                                                   cAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGGTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGaa             <  1:570850/250‑1 (MQ=255)
                                                                                                                                                                                                                                    aaaaaaCCAGGCTTGATTATAGCCTGGTTTCGTTTGTTTGGCTGAGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGA‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGACTATCCCCCTATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCTCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTgagactg           >  2:362852/1‑245 (MQ=255)
                                                                                                                                                                                                                                       aaaCCAGGCTTGATTATATCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCAGTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGAc        <  2:379925/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAgg    >  1:37232/1‑249 (MQ=255)
                                                                                                                                                                                                                                             ggCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCa                                           <  1:645006/211‑1 (MQ=255)
                                                                                                                                                                                                                                             ggCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCa                                           >  2:645006/1‑211 (MQ=255)
                                                                                                                                                                                                                                             ggCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGATTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAggac  >  2:307966/1‑248 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
TAACGCCAGCGGTGTCAGGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑AAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATC‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAAC  >  NC_000913/65571‑66055

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