Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,621,278 C→T F132F (TTC→TTT deoD → purine nucleoside phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,621,2780CT98.1% 159.6 / ‑4.4 52F132F (TTC→TTTdeoDpurine nucleoside phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (33/18);  total (33/19)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

AAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCA  >  NC_000913/4621043‑4621522
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    
aaGGCCGCAAAATTTCCGTAATTGGTCACGGTCTGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCCGTTCTGCCTCACGTAACACTGCGCGACGTCGTTCTCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCTCTATCGCTGACTTTGACATGGTCCGTAAc                                                                                                                                                                                                                                       >  2:393051/1‑251 (MQ=255)
  ggCCGCAAACTTTCCGTAATGGGTCACGGTATGCGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGCTCACCGATTTCGGCGTGCAGAAAATTATCCGCGTGGGTTCCTGTGGCGCTGTTCTGCCGCACGTACAACTGCGCGACGTCGTTATCGGTCTGGGTGCCTGCACCGATTCCACAGTTAACCGTATCCGTTTTAAAGGCAATGACTTTGCCGCTATCACTGGCTTTGACCTGGTGGGTAACgg                                                                                                                                                                                                                                     >  2:216769/1‑250 (MQ=255)
     cGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCGTCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGACGCTATCGCTGACTTTGACATGGTGCGTAACGcggt                                                                                                                                                                                                                                  >  1:148278/1‑248 (MQ=255)
         aaaatcccatAATTGGTCACTGTATTGGTATCACGTCCAGCTCCATCTACACCAAAGAACTGATCACCGATTGCGGCGTGAAGAAAATTATCCGCGTGGGTGCCTGTGGCGTATTTCTGCCGCACGTAAAACTTCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCGATCGCTGACTTTGACAGGGTGCGTACCGCGGTAGac                                                                                                                                                                                                                              <  2:320015/242‑2 (MQ=255)
                cGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAACTTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGATTTCGCTGACTTTGACATGGTGCGTAACGCGGTAGacgaggct                                                                                                                                                                                                                        >  2:171812/1‑243 (MQ=255)
                           aCGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTTAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATTGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGt                                                                                                                                                                                                            <  1:315684/251‑1 (MQ=255)
                                      aTCCCGTCCTGCTCCATCTACACCAAAGAACTGATCCCCGAGTGCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTTCCGCACGTAAAACTGCGCGACGTCGTTAGCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTTTCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTc                                                                                                                                                                                                  <  2:294848/250‑1 (MQ=255)
                                             ccTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGt                                                                                                                                                                                          >  1:326775/1‑251 (MQ=255)
                                               tGCTCCATCTACAACAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTa                                                                                                                                                                                         >  1:367154/1‑250 (MQ=255)
                                                gCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTGGAGGCTCGCGTGGGTAAc                                                                                                                                                                                       <  1:142056/251‑1 (MQ=255)
                                                  tCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGAGTGGGTAACCt                                                                                                                                                                                     >  1:340995/1‑251 (MQ=255)
                                                   ccATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGGTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTg                                                                                                                                                                                    <  2:469876/251‑1 (MQ=255)
                                                        tACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTtctc                                                                                                                                                                               >  1:362736/1‑251 (MQ=255)
                                                              aaaGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTg                                                                                                                                                                          >  1:106409/1‑250 (MQ=255)
                                                               aaGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGAc                                                                                                                                                                        >  1:397611/1‑251 (MQ=255)
                                                                aGAACTGATCCCCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGAcc                                                                                                                                                                       <  1:380078/251‑1 (MQ=255)
                                                                           ccGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGGTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGGGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTActct                                                                                                                                                            <  1:487621/251‑1 (MQ=255)
                                                                            cGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGCCCTGTTCTACtctc                                                                                                                                                           >  1:226388/1‑251 (MQ=255)
                                                                                            aaaaTTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGtt                                                                                                                                           <  1:210113/251‑1 (MQ=255)
                                                                                                 tATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGAc                                                                                                                                       <  2:106409/250‑1 (MQ=255)
                                                                                                     cgcgTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAACGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAATCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGAt                                                                                                                                  >  1:182984/1‑251 (MQ=255)
                                                                                                     cgcgTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGCCCTGTTCTACTCTCCGGACGGCGCACTGTTCgccgcga                                                                                                                                   >  2:100621/1‑244 (MQ=255)
                                                                                                               ccTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATAc                                                                                                                        >  1:340948/1‑251 (MQ=255)
                                                                                                               ccTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGGTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTGTTAAAGACCATTACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGTTGCTCGCGTGGGTAACCTGTTCTCCGCGGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATAc                                                                                                                        <  1:501456/251‑1 (MQ=255)
                                                                                                                  gtgGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATgg                                                                                                                                >  1:308409/1‑240 (MQ=255)
                                                                                                                  gtgGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATgg                                                                                                                                <  2:308409/240‑1 (MQ=255)
                                                                                                                      cgcAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCAtcc                                                                                                                 <  1:277667/251‑3 (MQ=255)
                                                                                                                          gTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTCGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGCTGGAAGAATACTGCATCCTCgg                                                                                                             >  2:116660/1‑251 (MQ=255)
                                                                                                                                         aaaCTGCTCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCAGCCGTTGTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCGAAAGCGCTGGGTATTGATGCGCGCGGGGGTAACCTGTGCTCCGCTGACCTGTTCGACTCGCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGGGGAAATGGAAGc                                                                                              <  2:49904/251‑1 (MQ=255)
                                                                                                                                            cTGCGCGACGTCGTTATCGGTATGGGAGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCACTGACCTGTTCTACTCTCCGGACGGCGAACTGTTCGCCGTGGTGGGAAAATAGCGCATCCTCGGCGGGGAAATGGATGCGGc                                                                                           >  2:198614/1‑251 (MQ=255)
                                                                                                                                                     gtcgtTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGCTGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCt                                                                                   >  1:83385/1‑250 (MQ=255)
                                                                                                                                                          tATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGGGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTTTACTCTCCGGCCGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGCACTGGAAGCGGCTggggtttatggcg                                                                             >  2:356145/1‑240 (MQ=255)
                                                                                                                                                             cGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGCCCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGGATCTAAGGCGTcc                                                                          >  2:71460/1‑250 (MQ=255)
                                                                                                                                                                  tGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGc                                                                      >  1:258474/1‑250 (MQ=255)
                                                                                                                                                                      ggCCTTCACCGATTCCAAAGTTAACCGTATCCGGTTTAAAGACCATGACGTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGGTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCGCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCGGCAGagt                                                                 <  2:207056/250‑3 (MQ=255)
                                                                                                                                                                                 aTTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGaaaa                                                                                                                            <  1:497244/181‑1 (MQ=255)
                                                                                                                                                                                 aTTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGaaaa                                                                                                                            >  2:497244/1‑181 (MQ=255)
                                                                                                                                                                                              ccgtttCCGTTTTAAAGACCATGACTTTGCCGTTATCGCTGACTTTGACATGGTGCGAGCCGCTGTTGACGCGGCTAAAGCCCTGGGGTTTTATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTTGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCagct                                         <  2:360155/246‑4 (MQ=255)
                                                                                                                                                                                                     cGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGt                                  >  1:354334/1‑251 (MQ=255)
                                                                                                                                                                                                              gACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGAccc                         >  1:395023/1‑250 (MQ=255)
                                                                                                                                                                                                               aCCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTACCGCGGTAGACGCTGCTAAAGCACTGGGGGTTGATGCTCGCGTGTGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCAGCCTCGGCGTGGAAATGGAAGCGGCTGGTGGCTACCGCCTCGCTGCAGAGTTTTGCGCGAACGCCCTTACCATCTgacacgtctctgaaccc                        >  1:116767/1‑235 (MQ=255)
                                                                                                                                                                                                               aCCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAACGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTCCGGGGTCGCTGCCGAGTTTGGCGCGAAAGCCCCGACCCTCTGCCCCGTGTCTgcaccc                        >  2:247808/1‑246 (MQ=255)
                                                                                                                                                                                                               aCCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGGCGCTGCGGAGTTTGGCGCGAAAGCCCTGGCCATCTGCACCGTGCCTGACcac                        >  2:133097/1‑251 (MQ=255)
                                                                                                                                                                                                                       tttGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCAc                >  1:243764/1‑251 (MQ=255)
                                                                                                                                                                                                                             gCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGGTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTAAGGCGTCGCTGAAGAGTTTGGCGCGAAAGCCCTTAACATCTGGACCggggcttaccccctcctccctcactg          >  2:231353/1‑226 (MQ=255)
                                                                                                                                                                                                                               tATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACgcgc        >  1:77924/1‑248 (MQ=255)
                                                                                                                                                                                                                                aTCGCTGACTTTGACATGGGGCGGAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCGGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGGTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGGGTCTGACCACATCCGCACTCGCGAGc        <  2:471840/250‑1 (MQ=255)
                                                                                                                                                                                                                                  cGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGtt                                                                                                                                           <  2:375907/117‑1 (MQ=255)
                                                                                                                                                                                                                                  cGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGtt                                                                                                                                           >  1:375907/1‑117 (MQ=255)
                                                                                                                                                                                                                                  cGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAg      <  2:326775/250‑1 (MQ=255)
                                                                                                                                                                                                                                     tGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGAccc  >  1:43323/1‑250 (MQ=255)
                                                                                                                                                                                                                                      gACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGGAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGGCGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGTATCTACGGCGTCGCTGCAg                                                                    >  1:123652/1‑183 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    
AAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCA  >  NC_000913/4621043‑4621522

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: