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→GAGdeoDpurine 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

CCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTG  >  NC_000913/4621211‑4621693
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                           
ccTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGttt                                                                                                                                                                                                                                          >  1:43502/1‑251 (MQ=255)
      ccGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGCAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGcgcg                                                                                                                                                                                                                                    >  1:300355/1‑251 (MQ=255)
             cAAAGTTAACCGTATCCATTTTAAAGATCATGACGTTGTCGCTATCGTTGACTTTGACAGGGTGAGTAACGCTGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCTGCTGATCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATTGAAGCTGCTGGTATCTACGGCGTAGATGCAGAGTTTGGCGCGAAAGccc                                                                                                                                                                                                                             <  2:472841/251‑1 (MQ=255)
                      ccgtaTCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGCCGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGCCGTGATGGACAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCTTCGCTGCAGAGTTTGGCGCGAAAGCCCTGCCCATCt                                                                                                                                                                                                                    >  2:462255/5‑251 (MQ=255)
                                tttAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTCACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAACGGCTGGTATTTACGGCGGTGCTGCAGAGTTTGGCCCGGAAGCCCTGACCctctcaccctggtc                                                                                                                                                                                                          >  2:529508/1‑241 (MQ=255)
                                tttAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGtgtc                                                                                                                                                                                                          >  1:370306/1‑248 (MQ=255)
                                   aaaGTCCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTg                                                                                                                                                                                                        <  1:86812/250‑1 (MQ=255)
                                   aaaGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGAAACGCGGTAGACGCGGATAAAGCACTGGGTGTTGATGCTCGCGTGGGAAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAACATACGGAATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTAAGGCGCCGCCGCAGAGTTTGGCGCGAAAGCCCTGTCCATCCGCACCgggtctgc                                                                                                                                                                                                       >  1:471707/1‑244 (MQ=255)
                                      gACCATGACTTGGCCGCTTACGCCGCTTTCGACAAGATGCGTAACGCCGTCGACGCGGCTAAAGCACTGGGGGATGATGCTCGCGTGGGTAACCTGGTCACAGTTGTCCCGTTCTACTCTCCGGCCGGCGAAATGTTCGACGTGATGGAAAAATAGGGCATCCCCGGCGCGGAAATGGAAGCGGCTGGTACCCTCGGCGTCGCTTCAGAGTTTGGCTCGAAAGTCCTGACCATCTGCTCCGTTTCTGTCca                                                                                                                                                                                                    <  2:35702/251‑1 (MQ=255)
                                         cATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATTCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGGGGAAATGGAAGCGGCTGGTATCTAAGGCGTCACTGCAGAGTTTGGGGCGAAAACCCTGAACATCTGCACCgagtcaggccacac                                                                                                                                                                                                 >  2:599847/1‑238 (MQ=255)
                                                  gCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTc                                                                                                                                                                                         <  1:127473/250‑1 (MQ=255)
                                                              gtctttgACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCGACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCa                                                                                                                                                                             <  2:50365/244‑1 (MQ=255)
                                                                         ggtgcgtaacgcgggagcctcggCTAGAGCATTGGGTTTTGATTCCCGCGTGTGTAACCTGTTCCCCGCTGACCGGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATTCGGCATCCTCTGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCTCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGGCCACTTCCGCACTCACGAGCAGACCACTGCCGCTGaa                                                                                                                                                                  <  2:471707/228‑2 (MQ=255)
                                                                              gtaccgctgTAGACGCGGCTAAAGCACTGGTTGTTTATGCCCGCGTGGGTAACCTGTTCTCCGCTGACCTTTTTTCCTCCCCGTCCGGCGTAATTGTCGCCGTGATGGAAAAATCCGGCATCCCCGGCGTGGAAATGGAAGCGTCTGGTATCTCCGGTGTCGCTTCTGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCTCTTCCGCCCTCACGAGCAGACCACTTCCGCTGAGCGTCAg                                                                                                                                                            <  2:422671/243‑1 (MQ=255)
                                                                              gTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGCGCGTc                                                                                                                                                              >  1:229620/1‑249 (MQ=255)
                                                                                aacgtggtagactcggCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTTACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCGCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGa                                                                                                                                                           <  2:323818/235‑1 (MQ=255)
                                                                                aaCGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGAc                                                                                                                                                          >  1:587402/1‑251 (MQ=255)
                                                                                 aCGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGAcc                                                                                                                                                         <  1:462255/251‑2 (MQ=255)
                                                                                       gagacgcggCTAAAGCCCTGGGTGTGGATGCTCGCTGGGGTAGCCTGTTCTCCGCTGACCTGTTCTACTCTCCTGACGGCGAAATGGTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTGTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGTCGCTGAGCGTCAGACCACCTTc                                                                                                                                                   <  2:586373/243‑1 (MQ=255)
                                                                                                   aaGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGGGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCTCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATc                                                                                                                                       <  1:293534/251‑1 (MQ=255)
                                                                                                    aGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATc                                                                                                                                       >  1:12557/1‑250 (MQ=255)
                                                                                                              tttGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGGTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTAGCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCAc                                                                                                                             <  2:496489/249‑1 (MQ=255)
                                                                                                              gttGATGCTCGCGTGGGTAACCTGTGCTCCGCTGACCTTTTCTACTCGCCGGACGGCGAAATGTTCGGCGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTCCGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCAc                                                                                                                             <  2:370306/249‑1 (MQ=255)
                                                                                                              ggtGATGCGCGCGTGGGTATCCTGTTCTCCGCTGACCTGGTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGTGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGATCAGACCACTGCCGCTGTGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACt                                                                                                                            <  2:344305/249‑1 (MQ=255)
                                                                                                                 gccgacctcgTGGGTAACCTGTTCTCCGCTGACCTGTTCTACCCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCAGCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCCCCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGa                                                                                                                         <  2:517503/243‑1 (MQ=255)
                                                                                                                 gATGCTCGCGTGGGTAACCTGTGCGCCGCTGACCTGGTCTACTATCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCGCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTTGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCTCACTgg                                                                                                                          <  2:169865/250‑1 (MQ=255)
                                                                                                                  attcTCGCGGGTGTATCCTGTTCTACTCTGACCTGTTATACTATCCGGACGGCGAAAGTTTCGACGTGCTGGAAAAATACGGATTGCTCGGCTTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGACAGCCCCTACAATCTGCACCGTTTCTAACCACATCCGCACTCCCGAGCAGACCACTGCCGCTGAGCGTCAGACCAACTTAAACGACATGATCAAAATCGCGCTGGaa                                                                                                                        <  2:55135/248‑1 (MQ=255)
                                                                                                                         cgTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTc                                                                                                                 <  1:273699/251‑1 (MQ=255)
                                                                                                                         cgTGGGTAACCTGTTCTACGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGa                                                                                                                         >  1:517503/1‑243 (MQ=255)
                                                                                                                          ggggggatccggTTCTCCGCTGACCTGGTCGACTCTCCGGACGGGGAAATGTGCGACGTGATGGAATAATACGGCAGCCTGGGCGGGGAAATGGAAGCGGCTGGGATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCGCCTGCAACGACATGATCAAAATCGCACTGGAAGCCGTTct                                                                                                                <  2:325725/240‑1 (MQ=255)
                                                                                                                          gTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAAt                                                                                                                       >  1:69877/1‑244 (MQ=255)
                                                                                                                          gTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCCCCGGACGGCGAAATGGTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAAt                                                                                                                       <  2:69877/244‑1 (MQ=255)
                                                                                                                                 cccgatcttcctcTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAAGACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGGGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCACACTGGAATCCGTTCTGCTGGGc                                                                                                         <  2:590111/239‑1 (MQ=255)
                                                                                                                                    tGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGAt                                                                                                      >  1:555828/1‑251 (MQ=255)
                                                                                                                                         tcCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCACACTGGAATCCGTTCTGCTGGGc                                                                                                         >  1:590111/1‑243 (MQ=255)
                                                                                                                                           cGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGt                                                                                               <  1:599847/251‑1 (MQ=255)
                                                                                                                                               gACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAAt                                                                                            <  2:576584/250‑1 (MQ=255)
                                                                                                                                                    gtgcTACTCGCCTGACGGCGAAGTGTGTGACGTGATGGAAAAAGACGGCATCCTCGGGGTGGTTATGGAAGCGGCTGTTATCTACGGCGTCGCTGCAGTGTGTGGCGCGAAAGCCCCGACCATCTGCACCGGGTCTGACCACATCCGCGCTCACGGGCAGACCACTGCCGCTGAGCTTCAGAGCACCTTCAACGACATGATGGAAAGCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTgtgt                                                                                       <  2:237976/247‑1 (MQ=255)
                                                                                                                                                    ggtCTACTCTCCGGACAGCGAAATGTTCGACGTGGTGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTCTCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGACTGACCACATCCGCACTAACGAGCAGACCACTGCCGCTGAGCGTCAGAACACCTTCAACGACATGATCAAAATCGCACTGGCATCCGTTCCGCCGGGCGATaaatggaaattatgtt                                                                                      >  2:489428/3‑238 (MQ=255)
                                                                                                                                                      tCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGttt                                                                                     <  1:394583/250‑1 (MQ=255)
                                                                                                                                                              ccGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTgaa                                                                              >  1:447675/1‑247 (MQ=255)
                                                                                                                                                               cggacgcagAAATGTTCTACGTGCTGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACTGCCTCGCTGCCGAGTTTGGAGCGCAAGCCCTGACCATCTGCACCGTGTCTGACCACCTCCGCACTCAAGCGCAGACCACTACCGATGAGCGCCAGACCACCTTCCACGACATGATCAAAATCGCAATGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTATGTTTCGCtcaaagg                                                                           >  2:429925/9‑245 (MQ=255)
                                                                                                                                                                ggACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTggg                                                                                                          >  1:601968/1‑219 (MQ=255)
                                                                                                                                                                ggACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACAGCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCTCACTGGAATCCGTTCTGCTggg                                                                                                          <  2:601968/219‑1 (MQ=255)
                                                                                                                                                                  aCGGCGAAATGTTCGACGTTATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGATTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATTAAAATCGCACGGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGa                                                                        <  2:45817/251‑1 (MQ=255)
                                                                                                                                                                      cGAAATGTTCGACGTGATGGCAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACGTGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATACAGAGTAATTGTGTTTCGCTGAAAGGCGaaagt                                                                    >  2:597808/1‑247 (MQ=255)
                                                                                                                                                                               cGACGTGATGGAAAAATACGGCATCCTCGGCGGGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGGTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACGGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAAACCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGaacgtcgtgttaag                                                           <  2:5120/251‑14 (MQ=255)
                                                                                                                                                                                    tGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGCTAAAGAGTAATTGTGTTTCGCTGAAAGGCGaacgtcgtgttaagccgg                                                       >  1:298655/1‑233 (MQ=255)
                                                                                                                                                                                        ggAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGcag                                                  <  2:587402/251‑3 (MQ=255)
                                                                                                                                                                                           aaaaaatcgggctcctcGGCGTGGGAAATGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGga                                                <  1:605187/237‑1 (MQ=255)
                                                                                                                                                                                           aaaaTACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGc                                                                                  <  1:61886/216‑1 (MQ=255)
                                                                                                                                                                                           aaaaTACGGCAGCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTGTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGgag                                               <  1:342470/251‑1 (MQ=255)
                                                                                                                                                                                           aaaaTAAGCCATCCTCGGCGTGGAAATTGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGAGCGATAGCCCTGACCATCTGCACCGTGTCTGACCACAGCAGGAATCACGAGCAGACCACTGCCGCTGGGCGTCAGACCACCTTCACCGACATGATCAAAATCGCACTGGAATCCGATCTGCTGGGCGTTAAAGAGTATTTGTGTTTCGc                                                                                  >  2:61886/1‑216 (MQ=255)
                                                                                                                                                                                                        cctcGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCt                                   >  1:22900/2‑251 (MQ=255)
                                                                                                                                                                                                              ggTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCAGCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAg                             <  2:71038/250‑1 (MQ=255)
                                                                                                                                                                                                              cGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAAt                                                                                                                                  >  1:417368/1‑150 (MQ=255)
                                                                                                                                                                                                                  gAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATc                         <  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)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                           
CCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTG  >  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: