Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,621,491 A→G V203V (GTA→GTG)  deoD → purine nucleoside phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,621,4910AG96.7% 193.3 / NA 60V203V (GTA→GTG) deoDpurine nucleoside phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (20/38);  minor base T (0/2);  total (20/40)
Fisher's exact test for biased strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

AAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATT‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTG  >  NC_000913/4621246‑4621693
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                          
aaaGTCCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATC‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTg                                                                                                                                                                                                        <  1:86812/250‑1 (MQ=255)
   gACCATGACTTGGCCGCTTACGCCGCTTTCGACAAGATGCGTAACGCCGTCGACGCGGCTAAAGCACTGGGGGATGATGCTCGCGTGGGTAACCTGGTCACAGTTGTCCCGTTCTACTCTCCGGCCGGCGAAATGTTCGACGTGATGGAAAAATAGGGCATC‑CCCGGCGCGG‑AAATGGAAGCGGCTGGTACCCTCGGCGTCGCTTCAGAGTTTGGCTCGAAAGTCCTGACCATCTGCTCCGTTTCTGTCca                                                                                                                                                                                                    <  2:35702/251‑1 (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                                                                                                                                                     >  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)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                          
AAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATT‑CTCGGCGTGG‑AAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTG  >  NC_000913/4621246‑4621693

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