Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,621,542 T→C T220T (ACT→ACC)  deoD → purine nucleoside phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,621,5420TC98.0% 157.1 / ‑4.2 51T220T (ACT→ACC) deoDpurine nucleoside phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (18/32);  total (19/32)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

TAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATT‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCT  >  NC_000913/4621298‑4621748
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                              
gagacgcggCTAAAGCCCTGGGTGTGGATGCTCGCTGGGGTAGCCTGTTCTCCGCTGACCTGTTCTACTCTCCTGACGGCGAAATGGTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTGTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGTCGCTGAGCGTCAGACCACCTTc                                                                                                                                                                                                          <  2:586373/243‑1 (MQ=255)
            aaGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGGGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCTCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATc                                                                                                                                                                                              <  1:293534/251‑1 (MQ=255)
             aGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATc                                                                                                                                                                                              >  1:12557/1‑250 (MQ=255)
                       tttGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGGTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTAGCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCAc                                                                                                                                                                                    <  2:496489/249‑1 (MQ=255)
                       gttGATGCTCGCGTGGGTAACCTGTGCTCCGCTGACCTTTTCTACTCGCCGGACGGCGAAATGTTCGGCGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTCCGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCAc                                                                                                                                                                                    <  2:370306/249‑1 (MQ=255)
                       ggtGATGCGCGCGTGGGTATCCTGTTCTCCGCTGACCTGGTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGTGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGATCAGACCACTGCCGCTGTGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACt                                                                                                                                                                                   <  2:344305/249‑1 (MQ=255)
                          gccgacctcgTGGGTAACCTGTTCTCCGCTGACCTGTTCTACCCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCAGC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCCCCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGa                                                                                                                                                                                <  2:517503/243‑1 (MQ=255)
                          gATGCTCGCGTGGGTAACCTGTGCGCCGCTGACCTGGTCTACTATCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CGCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTTGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCTCACTgg                                                                                                                                                                                 <  2:169865/250‑1 (MQ=255)
                           attcTCGCGGGTGTATCCTGTTCTACTCTGACCTGTTATACTATCCGGACGGCGAAAGTTTCGACGTGCTGGAAAAATACGG‑ATTGCTCGGCTTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGACAGCCCCTACAATCTGCACCGTTTCTAACCACATCCGCACTCCCGAGCAGACCACTGCCGCTGAGCGTCAGACCAACTTAAACGACATGATCAAAATCGCGCTGGaa                                                                                                                                                                               <  2:55135/248‑1 (MQ=255)
                                  cgTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTc                                                                                                                                                                        <  1:273699/251‑1 (MQ=255)
                                  cgTGGGTAACCTGTTCTACGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGa                                                                                                                                                                                >  1:517503/1‑243 (MQ=255)
                                   ggggggatccggTTCTCCGCTGACCTGGTCGACTCTCCGGACGGGGAAATGTGCGACGTGATGGAATAATACGGCAGC‑CTGGGCGGGG‑AAATGGAAGCGGCTGGGATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCGCCTGCAACGACATGATCAAAATCGCACTGGAAGCCGTTct                                                                                                                                                                       <  2:325725/240‑1 (MQ=255)
                                   gTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAAt                                                                                                                                                                              >  1:69877/1‑244 (MQ=255)
                                   gTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCCCCGGACGGCGAAATGGTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAAt                                                                                                                                                                              <  2:69877/244‑1 (MQ=255)
                                          cccgatcttcctcTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAAGACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGGGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCACACTGGAATCCGTTCTGCTGGGc                                                                                                                                                                <  2:590111/239‑1 (MQ=255)
                                             tGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGAt                                                                                                                                                             >  1:555828/1‑251 (MQ=255)
                                                  tcCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCACACTGGAATCCGTTCTGCTGGGc                                                                                                                                                                >  1:590111/1‑243 (MQ=255)
                                                    cGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGt                                                                                                                                                      <  1:599847/251‑1 (MQ=255)
                                                        gACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAAt                                                                                                                                                   <  2:576584/250‑1 (MQ=255)
                                                             gtgcTACTCGCCTGACGGCGAAGTGTGTGACGTGATGGAAAAAGACGGCATC‑CTCGGGGTGG‑TTATGGAAGCGGCTGTTATCTACGGCGTCGCTGCAGTGTGTGGCGCGAAAGCCCCGACCATCTGCACCGGGTCTGACCACATCCGCGCTCACGGGCAGACCACTGCCGCTGAGCTTCAGAGCACCTTCAACGACATGATGGAAAGCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTgtgt                                                                                                                                              <  2:237976/247‑1 (MQ=255)
                                                             ggtCTACTCTCCGGACAGCGAAATGTTCGACGTGGTGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTCTCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGACTGACCACATCCGCACTAACGAGCAGACCACTGCCGCTGAGCGTCAGAACACCTTCAACGACATGATCAAAATCGCACTGGCATCCGTTCCGCCGGGCGATaaatggaaattatgtt                                                                                                                                             >  2:489428/3‑238 (MQ=255)
                                                               tCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGttt                                                                                                                                            <  1:394583/250‑1 (MQ=255)
                                                                       ccGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTgaa                                                                                                                                     >  1:447675/1‑247 (MQ=255)
                                                                        cggacgcagAAATGTTCTACGTGCTGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACTGCCTCGCTGCCGAGTTTGGAGCGCAAGCCCTGACCATCTGCACCGTGTCTGACCACCTCCGCACTCAAGCGCAGACCACTACCGATGAGCGCCAGACCACCTTCCACGACATGATCAAAATCGCAATGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTATGTTTCGCtcaaagg                                                                                                                                  >  2:429925/9‑245 (MQ=255)
                                                                         ggACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTggg                                                                                                                                                                 >  1:601968/1‑219 (MQ=255)
                                                                         ggACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACAGCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCTCACTGGAATCCGTTCTGCTggg                                                                                                                                                                 <  2:601968/219‑1 (MQ=255)
                                                                           aCGGCGAAATGTTCGACGTTATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGATTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATTAAAATCGCACGGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGa                                                                                                                               <  2:45817/251‑1 (MQ=255)
                                                                               cGAAATGTTCGACGTGATGGCAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACGTGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATACAGAGTAATTGTGTTTCGCTGAAAGGCGaaagt                                                                                                                           >  2:597808/1‑247 (MQ=255)
                                                                                        cGACGTGATGGAAAAATACGGCATC‑CTCGGCGGGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGGTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACGGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAAACCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGaacgtcgtgttaag                                                                                                                  <  2:5120/251‑14 (MQ=255)
                                                                                             tGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGCTAAAGAGTAATTGTGTTTCGCTGAAAGGCGaacgtcgtgttaagccgg                                                                                                              >  1:298655/1‑233 (MQ=255)
                                                                                                 ggAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGcag                                                                                                         <  2:587402/251‑3 (MQ=255)
                                                                                                    aaaaaatcgggctcctcGGCGTGGGAAATG‑AAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGga                                                                                                       <  1:605187/237‑1 (MQ=255)
                                                                                                    aaaaTACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGc                                                                                                                                         <  1:61886/216‑1 (MQ=255)
                                                                                                    aaaaTACGGCAGC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTGTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGgag                                                                                                      <  1:342470/251‑1 (MQ=255)
                                                                                                    aaaaTAAGCCATC‑CTCGGCGTGG‑AAATTGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGAGCGATAGCCCTGACCATCTGCACCGTGTCTGACCACAGCAGGAATCACGAGCAGACCACTGCCGCTGGGCGTCAGACCACCTTCACCGACATGATCAAAATCGCACTGGAATCCGATCTGCTGGGCGTTAAAGAGTATTTGTGTTTCGc                                                                                                                                         >  2:61886/1‑216 (MQ=255)
                                                                                                                 cctcGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCt                                                                                          >  1:22900/2‑251 (MQ=255)
                                                                                                                       ggTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCAGCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAg                                                                                    <  2:71038/250‑1 (MQ=255)
                                                                                                                       cGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAAt                                                                                                                                                                                         >  1:417368/1‑150 (MQ=255)
                                                                                                                           g‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATc                                                                                <  1:114074/251‑1 (MQ=255)
                                                                                                                             aaaTGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATCt                                                                               >  1:673515/1‑251 (MQ=255)
                                                                                                                              aaTGGAAGCGGCTGGTATCTCCGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATCTa                                                                              <  2:122353/251‑1 (MQ=255)
                                                                                                                                 ggAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATCTAtt                                                                            <  1:498748/250‑1 (MQ=255)
                                                                                                                                       ggCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCt                                                                                                                                                                                                            >  1:591050/1‑116 (MQ=255)
                                                                                                                                       ggCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCt                                                                                                                                                                                                            <  2:591050/116‑1 (MQ=255)
                                                                                                                                           ggTATCTACGGTGTCGCTTCAGGGTTTGGCGCGAATGCCCGGACCATCTGCTCCTGGGCTGACCACATCCGCCCTCACGCGCAGACCCCTGTCGCTGATCGTCATACCACTTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGATATGTGAACGTCGTGTTAAGCCGGAGCCGGGGACTGCTCCGGCTTTTTAGTAGCTATTcttcaatatc                                                                  <  2:535636/250‑10 (MQ=255)
                                                                                                                                              aTCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATCTATTCAtcaatagcaac                                                               <  1:444011/250‑11 (MQ=255)
                                                                                                                                                cTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCg                                                                                                                                                               >  1:612027/1‑152 (MQ=255)
                                                                                                                                                cTACGGCGTCGCGGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACAGCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACGGGAATCCGTGCTGCTGGGCg                                                                                                                                                               <  2:612027/152‑1 (MQ=255)
                                                                                                                                                  aCGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATCTATTCAtcaatagcaactga                                                            <  1:679443/249‑14 (MQ=255)
                                                                                                                                                    ggCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATCTATTCAtcaatagcaactgattt                                                         >  1:384982/1‑234 (MQ=255)
                                                                                                                                                                                                          cacaTCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGACGATTGCCTTGTGAAGCCGGAGCGGGCGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGACTATTTTCGCTATACTTTTCAGTGAAATGTGTTCATAAATCTATTCAAGTATCt  >  2:331008/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                              
TAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATT‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCT  >  NC_000913/4621298‑4621748

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