Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,621,458 A→G E192E (GAA→GAG)  deoD → purine nucleoside phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,621,4580AG98.5% 203.2 / NA 68E192E (GAA→GAG) deoDpurine nucleoside phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (27/40);  minor base T (0/1);  total (27/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATT‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTG  >  NC_000913/4621211‑4621693
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                           
ccTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGttt                                                                                                                                                                                                                                          >  1:43502/1‑251 (MQ=255)
      ccGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑CAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGcgcg                                                                                                                                                                                                                                    >  1:300355/1‑251 (MQ=255)
             cAAAGTTAACCGTATCCATTTTAAAGATCATGACGTTGTCGCTATCGTTGACTTTGACAGGGTGAGTAACGCTGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCTGCTGATCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATTGAAGCTGCTGGTATCTACGGCGTAGATGCAGAGTTTGGCGCGAAAGccc                                                                                                                                                                                                                             <  2:472841/251‑1 (MQ=255)
                      ccgtaTCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGCCGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGCCGTGATGGACAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCTTCGCTGCAGAGTTTGGCGCGAAAGCCCTGCCCATCt                                                                                                                                                                                                                    >  2:462255/5‑251 (MQ=255)
                                tttAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTCACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAACGGCTGGTATTTACGGCGGTGCTGCAGAGTTTGGCCCGGAAGCCCTGACCctctcaccctggtc                                                                                                                                                                                                          >  2:529508/1‑241 (MQ=255)
                                tttAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGtgtc                                                                                                                                                                                                          >  1:370306/1‑248 (MQ=255)
                                   aaaGTCCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTg                                                                                                                                                                                                        <  1:86812/250‑1 (MQ=255)
                                   aaaGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGAAACGCGGTAGACGCGGATAAAGCACTGGGTGTTGATGCTCGCGTGGGAAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAACATACGGAATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTAAGGCGCCGCCGCAGAGTTTGGCGCGAAAGCCCTGTCCATCCGCACCgggtctgc                                                                                                                                                                                                       >  1:471707/1‑244 (MQ=255)
                                      gACCATGACTTGGCCGCTTACGCCGCTTTCGACAAGATGCGTAACGCCGTCGACGCGGCTAAAGCACTGGGGGATGATGCTCGCGTGGGTAACCTGGTCACAGTTGTCCCGTTCTACTCTCCGGCCGGCGAAATGTTCGACGTGATGGAAAAATAGGGCATC‑CCCGGCGCGG‑AAATGGAAGCGGCTGGTACCCTCGGCGTCGCTTCAGAGTTTGGCTCGAAAGTCCTGACCATCTGCTCCGTTTCTGTCca                                                                                                                                                                                                    <  2:35702/251‑1 (MQ=255)
                                         cATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATTCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGGGG‑AAATGGAAGCGGCTGGTATCTAAGGCGTCACTGCAGAGTTTGGGGCGAAAACCCTGAACATCTGCACCgagtcaggccacac                                                                                                                                                                                                 >  2:599847/1‑238 (MQ=255)
                                                  gCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTc                                                                                                                                                                                         <  1:127473/250‑1 (MQ=255)
                                                              gtctttgACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCGACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCa                                                                                                                                                                             <  2:50365/244‑1 (MQ=255)
                                                                         ggtgcgtaacgcgggagcctcggCTAGAGCATTGGGTTTTGATTCCCGCGTGTGTAACCTGTTCCCCGCTGACCGGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATTCGGCATC‑CTCTGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCTCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGGCCACTTCCGCACTCACGAGCAGACCACTGCCGCTGaa                                                                                                                                                                  <  2:471707/228‑2 (MQ=255)
                                                                              gtaccgctgTAGACGCGGCTAAAGCACTGGTTGTTTATGCCCGCGTGGGTAACCTGTTCTCCGCTGACCTTTTTTCCTCCCCGTCCGGCGTAATTGTCGCCGTGATGGAAAAATCCGGCATC‑CCCGGCGTGG‑AAATGGAAGCGTCTGGTATCTCCGGTGTCGCTTCTGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCTCTTCCGCCCTCACGAGCAGACCACTTCCGCTGAGCGTCAg                                                                                                                                                            <  2:422671/243‑1 (MQ=255)
                                                                              gTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGCGCGTc                                                                                                                                                              >  1:229620/1‑249 (MQ=255)
                                                                                aacgtggtagactcggCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTTACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CGCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGa                                                                                                                                                           <  2:323818/235‑1 (MQ=255)
                                                                                aaCGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGAc                                                                                                                                                          >  1:587402/1‑251 (MQ=255)
                                                                                 aCGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGAcc                                                                                                                                                         <  1:462255/251‑2 (MQ=255)
                                                                                       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                                                                                                                                                     <  2:591050/116‑1 (MQ=255)
                                                                                                                                                                                                                              ggCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCt                                                                                                                                                     >  1:591050/1‑116 (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)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                           
CCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATT‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTG  >  NC_000913/4621211‑4621693

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