Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,271,130 C→T E914E (GAG→GAA)  uvrA ← UvrABC excision nuclease subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,271,1300CT98.6% 242.2 / ‑3.3 70coding (2742/2823 nt)uvrAUvrABC excision nuclease subunit A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (41/28);  minor base . (1/0);  total (42/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

CTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCTCGTCGAGAATATACAGCGTCTGCCCGGTGCCGCGTTTTGACAGTTCACGCGCCAGCTTCACGCGCTGGGCT  >  NC_000913/4270888‑4271370
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
cTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCgccgcc                                                                                                                                                                                                                                          >  2:169472/1‑251 (MQ=255)
       aCCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTg                                                                                                                                                                                                                                   <  1:253991/251‑1 (MQ=255)
            gtATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTgccgcc                                                                                                                                                                                                                              <  1:457967/251‑1 (MQ=255)
                  cGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTgg                                                                                                                                                                                                                        <  1:139032/251‑1 (MQ=255)
                     cGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGACGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGAACTCCGCGCCGGTTTCTGGCGTACCGGAGACGA‑GGAATTCGCCGCCACTACTGCCGCCTTCTGGTcc                                                                                                                                                                                                                     >  1:295741/1‑251 (MQ=255)
                     cGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTcc                                                                                                                                                                                                                     >  1:220240/1‑251 (MQ=255)
                              gCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTCTCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGCTCCCAGGTCGAc                                                                                                                                                                                                            >  2:549573/1‑251 (MQ=255)
                                      cGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGGCCCAGGTCGACAATCCAg                                                                                                                                                                                                     <  2:369526/250‑1 (MQ=255)
                                               attatCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGttt                                                                                                                                                                                           >  2:14670/1‑251 (MQ=255)
                                               attagCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATGTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGTCTTCTGGTCCCAGGTCGACAATCCAGTCAGCgg                                                                                                                                                                                              <  2:489464/248‑1 (MQ=255)
                                                      ggtggATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCa                                                                                                                                                                                     <  1:36605/250‑1 (MQ=255)
                                                                   tATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGc                                                                                                                                                                       >  1:32787/1‑251 (MQ=255)
                                                                         caggctgtgtagttgctgctgGAAGAAAAACGTAAATTGCTGGGTCAACTCTGAAAGGAAAAGGCCTCGCCGAAAGCGGCCGTAGCGATTACAGCATAGGCTTAATGAAGCGTGCCGTGTGTGATGCTTTGCCCTCCGCGACGGTTTCTGGCGTATCTGAGGCGA‑GTATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCTACAATCCAGTCAGCGGTTTTGATCTCGTCGAGATTTTGCTCactc                                                                                                                                                                 <  2:153489/236‑4 (MQ=255)
                                                                                           ccGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCATAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCTTCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGAGTGTGCTCACTCACCACAATGGTGTTGccc                                                                                                                                               <  2:260010/251‑1 (MQ=255)
                                                                                             ggAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTg                                                                                                                                             >  2:198424/1‑251 (MQ=255)
                                                                                              gaagaaAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTg                                                                                                                                             >  2:6128/1‑250 (MQ=255)
                                                                                               aagaaAAACGTAAATTGCTGGTTCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATca                                                                                                                                                                <  1:570592/230‑1 (MQ=255)
                                                                                               aagaaAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGAt                                                                                                                                           <  2:501134/251‑1 (MQ=255)
                                                                                               aagaaAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTAAAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATca                                                                                                                                                                >  2:570592/1‑230 (MQ=255)
                                                                                                agaaAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATc                                                                                                                                          <  1:119109/251‑1 (MQ=255)
                                                                                                      aCGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGCTTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCgcgc                                                                                                                                      >  2:510486/1‑249 (MQ=255)
                                                                                                              tGCTGGTGCAACTCTGCAAGGAAAACGCCGCTAAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACACTCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGc                                                                                                                             >  2:48937/1‑250 (MQ=255)
                                                                                                               gatGGTGCAACTCTGAAAGGAAACTGCCGCTCAGAAAGCGGCCTTACCGATTACAGCATCGGCTTAAAGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCCGTTTAttc                                                                                                                             >  2:414272/3‑247 (MQ=255)
                                                                                                               gCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCa                                                                                                                            <  1:14670/250‑1 (MQ=255)
                                                                                                                cTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAg                                                                                                                                                                                                                  >  1:152827/1‑163 (MQ=255)
                                                                                                                cTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTGCTGGTCCCAg                                                                                                                                                                                                                  <  2:152827/163‑1 (MQ=255)
                                                                                                                   gTGCAACTCTGAAAGGTAAAGACCGCTAAGAAAGCGGCCGTAGCGATTACGGCATTGGGTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACTGTGTTTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACTACTGCCGGCTTCTGGGCCCAGGTCGACACTCCAGTCAGCGGTTTTGATCACGTCGAGATTTTGCTCAATCACCACAATGGTGTTGCTCTGATCGCGCAGTGTATGCAGTACg                                                                                                                       <  1:586557/251‑1 (MQ=255)
                                                                                                                       aaCTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCg                                                                                                                    >  1:348107/1‑250 (MQ=255)
                                                                                                                           ctGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGc                                                                                                                 >  1:267178/1‑249 (MQ=255)
                                                                                                                                    aaaGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAat                                                                                                      >  2:146924/1‑251 (MQ=255)
                                                                                                                                          cGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGCAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGATGTTGCCCTGCTCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGACTATCGGc                                                                                                >  2:316196/1‑251 (MQ=255)
                                                                                                                                            cTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCg                                                                                               >  2:319564/1‑250 (MQ=255)
                                                                                                                                              cAGAACGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGCGCAGTTGCTGAATATCGGCGa                                                                                              >  2:568623/1‑249 (MQ=255)
                                                                                                                                              cAGAAAGCGGCCGTAACGATTCCAGCATCGGTTTAAGGAAGCGTGCGGTGTGTGATGCTTCGCTCGCCGCGACGGTTTCGGGCGTACCGGAGGCGA‑GGATTGCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAAGCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGGTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAg                                                                                            <  2:518066/251‑1 (MQ=255)
                                                                                                                                               aGAAAGCGGCCTTAACGCTTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGt                                                                                           >  2:208663/1‑251 (MQ=255)
                                                                                                                                                    gCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGc                                                                                         >  1:315937/1‑248 (MQ=255)
                                                                                                                                                      ggCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACAACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGAcc                                                                                    >  2:566961/1‑251 (MQ=255)
                                                                                                                                                         cTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCggt                                                                                 >  1:297627/1‑251 (MQ=255)
                                                                                                                                                          ttAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCggtg                                                                                >  2:156467/1‑251 (MQ=255)
                                                                                                                                                           tAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCggtg                                                                                <  2:492625/250‑1 (MQ=255)
                                                                                                                                                             aCGATTACAGCCTCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTc                                                                             >  1:442523/1‑251 (MQ=255)
                                                                                                                                                             aCGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTc                                                                             <  1:156467/251‑1 (MQ=255)
                                                                                                                                                             aCGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTACGCGACGGTTTCTGGAGTAACGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCAGGTCCATGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGCGCAGATGCTGAATCTCGGCGAAGTGCAGACCGGTGGTc                                                                             >  1:525977/1‑251 (MQ=255)
                                                                                                                                                              cGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACGCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCg                                                                            <  2:210509/251‑1 (MQ=255)
                                                                                                                                                                   aCAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGCCGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCtcg                                                                       >  1:171081/1‑251 (MQ=255)
                                                                                                                                                                   aCAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCtc                                                                        <  1:169472/250‑1 (MQ=255)
                                                                                                                                                                    cAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCtcgt                                                                      >  1:32555/1‑251 (MQ=255)
                                                                                                                                                                      gCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCtcgtcg                                                                    <  1:319564/251‑1 (MQ=255)
                                                                                                                                                                       cATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGCGACGA‑GGATTTCGCCGCCACCCCTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGAGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCtcgtcg                                                                    >  1:380549/1‑250 (MQ=255)
                                                                                                                                                                          cGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGTCGAAGTGCAGACCGGTGGTCGGCTCGTCgaga                                                                 <  1:479382/250‑1 (MQ=255)
                                                                                                                                                                                aaGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGGGGTCGGCTCTTCGAGAATATACa                                                          <  2:239213/251‑1 (MQ=255)
                                                                                                                                                                                 aGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCTCGTCGAGAATATAc                                                           >  1:21601/1‑249 (MQ=255)
                                                                                                                                                                                     aGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCTCGTCGAGAATATACAGCGTc                                                     <  2:442523/251‑1 (MQ=255)
                                                                                                                                                                                       cGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCTCGTCGAGAATATACAGCGTCTg                                                   >  2:148160/1‑251 (MQ=255)
                                                                                                                                                                                           ccGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCTCGTCGAGAATATACAGCGTCTGCCCg                                               >  2:471634/1‑251 (MQ=255)
                                                                                                                                                                                              tgtgtgATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCTCGTCGAGAATATACAGCGTCTGCCCgggg                                            >  2:415566/1‑249 (MQ=255)
                                                                                                                                                                                                  tgATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCCCGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCTCGTCGAGAATATACAGCGTCTGCCCGGTGCcgc                                        >  1:567887/1‑251 (MQ=255)
                                                                                                                                                                                                   gATGCTTCGAACTCGGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTACGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGATTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGTTGACTATCGGCGGAGTGCAGACCGGTGGTTGGCTCGTCGAGAATCTACTGCGTCTTCCCGGTGAcgcg                                       >  2:379958/1‑251 (MQ=255)
                                                                                                                                                                                                                ccGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCTCGTCGAGAATATACAGCGTCTGCCCGGTGCCGCGTTTTGACAGtt                            <  1:280793/249‑1 (MQ=255)
                                                                                                                                                                                                                 cgcgACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGCGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCTCGTCGAGAATATACAGCGTCTGCCCGGTGCCGCGTTTTGACAGTTCAc                         >  2:262013/1‑251 (MQ=255)
                                                                                                                                                                                                                    gACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCTCGTCGAGAATATACAGCGTCTGCCCGGTGCCGCGTTTTGACAGTTCAcgcg                      >  1:490544/1‑251 (MQ=255)
                                                                                                                                                                                                                    gACGGTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCTCGTCGAGAATATACAGCGTCTGCCCGGTGCCGCGTTTTGACAGTTCAcgcg                      <  1:198424/251‑1 (MQ=255)
                                                                                                                                                                                                                        gTTTCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCTCGTCGAGAATATACAGCGTCTGCCCGGTGCCGCGTTTTGACAGTTCACGCGCCAg                  <  1:97581/251‑1 (MQ=255)
                                                                                                                                                                                                                         tttCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCTCGTCGAGAATATACAGCGTCTGCCCGGTGCCGCGTTTTGACAGTTCACGCGCCAGc                 >  1:72095/1‑251 (MQ=255)
                                                                                                                                                                                                                         tttCTGGCGTACCGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCTCGTCGAGAATATACAGCGTCTGCCCGGTGCCGCGTTTTGACAGTTCACGCGCCAGc                 >  1:405063/1‑251 (MQ=255)
                                                                                                                                                                                                                             tGGCGTACCGGAGCCGC‑GGATTTAGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACACTCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCAACACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCTCGTCGAGAATATACAGCGTCTGCCCGGTGCCGCGTTTTGAAAGTTCACGCGCCAGCTTc              >  2:42406/1‑250 (MQ=255)
                                                                                                                                                                                                                                   cccGGAGACGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCTCGTCGAGAATATACAGCGTCTGCCCGGTGCCGCGTTTTGACAGTTCACGCGCTAGCTTCACGCGCt       <  2:516142/250‑1 (MQ=255)
                                                                                                                                                                                                                                   aCCGGAGACGATGCAT‑TCGCCACCACCACTGCCGCCTTCTGCTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCCCGTCGGGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGATGATGAATATCGCCGAAGTGCAGACAGGTGGTGGGCTCGTCGAGAATATACCGCGTGTGCCCGGTGCCGCGTTTTAACACTTCACGCGCCAGCTTCAAGCGCt       >  2:127742/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cGGAGACGA‑GGATTTCGCCGCCACCCCTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTGTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCTCGTCGAGAATATACAGCGTCTGCCCGGTGCCGCGTTTTGACAGTTCACGCGCCAGCTTCACGCGCTgg     <  2:348107/251‑1 (MQ=255)
                                                                                                                                                                                                                                        agaCGA‑GGATTTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCTCGTCGAGAATATAAAGCGTCTGCCCGGTGCCGCGTTTTGACAGTTCACGCGCCAGCTTCAAGCGCTGGGCt  >  2:554953/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
CTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGA‑GGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTCGAGCAGTTGCTGAATATCGGCGAAGTGCAGACCGGTGGTCGGCTCGTCGAGAATATACAGCGTCTGCCCGGTGCCGCGTTTTGACAGTTCACGCGCCAGCTTCACGCGCTGGGCT  >  NC_000913/4270888‑4271370

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: ‑