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

AAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTG  >  NC_000913/4621246‑4621693
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                          
aaaGTCCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTg                                                                                                                                                                                                        <  1:86812/250‑1 (MQ=255)
   gACCATGACTTGGCCGCTTACGCCGCTTTCGACAAGATGCGTAACGCCGTCGACGCGGCTAAAGCACTGGGGGATGATGCTCGCGTGGGTAACCTGGTCACAGTTGTCCCGTTCTACTCTCCGGCCGGCGAAATGTTCGACGTGATGGAAAAATAGGGCATCCCCGGCGCGGAAATGGAAGCGGCTGGTACCCTCGGCGTCGCTTCAGAGTTTGGCTCGAAAGTCCTGACCATCTGCTCCGTTTCTGTCca                                                                                                                                                                                                    <  2:35702/251‑1 (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                                                                                                                                                     >  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)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                          
AAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTG  >  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: