Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,619,121 (CG)4→3 coding (893‑894/1323 nt) deoA → thymidine phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,619,1140G.98.2% 245.3 / ‑2.5 55coding (886/1323 nt)deoAthymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (27/27);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01
*NC_0009134,619,1150C.98.2% 245.8 / ‑3.2 55coding (887/1323 nt)deoAthymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (27/27);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAAT‑GGGCGGCGGACGCCGTCAGGCATCTGACACCA  >  NC_000913/4618864‑4619357
                                                                                                                                                                                                                                                          ||                                                                                                                                                                                                                                                     
gACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAg                                                                                                                                                                                                                                                        <  2:234217/251‑1 (MQ=255)
               cgAAGCCCGTGCCGACGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCGCCGCGCTGCTCACCGACACGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAGCCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCAGAAgcgcgcttaattgcag                                                                                                                                                                                                                                         <  2:49516/250‑11 (MQ=255)
                          ccGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGa                                                                                                                                                                                                                            <  1:134154/251‑1 (MQ=255)
                           cGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGAc                                                                                                                                                                                                                           <  1:5591/251‑1 (MQ=255)
                               gCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCCGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGCGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACg                                                                                                                                                                                                                       >  2:407683/1‑251 (MQ=255)
                                 gATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTCCTGGCCTACCGTGCAGGTAACGCGGTTGAAGCTCGTGAAGCGGTGCAGTTCCTGACGGGCGCATATCGCAACCCGCGCCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAACGATGACGCCGAAGCGCGCGC‑‑AATTGCAGGCGGAGCTGGACAACggc                                                                                                                                                                                                                     >  1:47606/1‑250 (MQ=255)
                                            cctacacgacgctcgcccgatcgaccaccGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATAGCGTAACCCGCGTCTGTTTGATGTCACGAGGGCGCTGTGCGGGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGGGCTGGACAACGGTAAAGCTGCAGa                                                                                                                                                                                                          <  2:421921/228‑1 (MQ=255)
                                                         tGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGt                                                                                                                                                                                            >  2:14783/1‑251 (MQ=255)
                                                                cgcACCACCGCGCTGCTCACCGACATGAATCAGGTACGGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAg                                                                                                                                                                                     <  2:216237/251‑1 (MQ=255)
                                                                 gcccgcccGCGCTGCTTACTGACATGAGTCAGGGACTGGCCTCGAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATTTCGTAACCCGCGTCTGTGTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCGCCGGCAAACTGGGGAGAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGATAACGGTAAAGCGGCATAAGTCT‑TTGGTCGTATGGTAGc                                                                                                                                                                                    <  2:292656/245‑1 (MQ=255)
                                                                   accaccGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCCGTTCCTGCCGGGTGAATATCGTACCCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGGAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTTGCCAACGGTAAAACGGCAGAAgctcgggcgagcacccgtctg                                                                                                                                                                                    >  1:421921/1‑230 (MQ=255)
                                                                                tCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGGAAAGCGGCAGAAGTCT‑TTGGGCGTAGGGGAGCGGGCA‑AAAAAGGccc                                                                                                                                                                    >  2:44567/1‑251 (MQ=255)
                                                                                      aCATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAAAGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CGCAAAAAGGCccgaccg                                                                                                                                                               >  1:100743/1‑250 (MQ=255)
                                                                                       cATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGAAAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCcgaccgac                                                                                                                                                             >  1:81281/1‑251 (MQ=255)
                                                                                         tGAATCAGTTCCTGGCCTCCAGTGCAGGGAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTACCTGACGGGTGAATATCGTAACCTGCGGCTGTTTTATGTCACGATGGCGCTGTGCGTGGAGATGCTGAACGCCGGCAAACTGGCGAAAGATGTCGCCGAA‑‑GCGCGCGAAATGGCAGGCGGTGCTGGACAACGGTAAGGCGTCATAAGTCTGTTGG‑CGTATGGCATCGG‑TACAAAAAGGGCCGTCCGACt                                                                                                                                                            <  2:303637/250‑1 (MQ=255)
                                                                                         tGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGACGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGGTTGATGTCACGATTGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGGTGGACAACGGGAAAGCGGGGGAAGACT‑TTGGGCGGTTGGTAGCGG‑GACAAAAAGGCCcgaacgacgt                                                                                                                                                           >  2:331217/1‑249 (MQ=255)
                                                                                         tGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAAAAGGCT‑TTGGGTGGATGGTAGCGG‑CACAAAAAGGCCCGACCTACtt                                                                                                                                                           >  1:25694/1‑251 (MQ=255)
                                                                                            aTCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGGAAAGCGGGAGAAGGCT‑TTGGTCGGATGGTAGAGG‑CACAAAAAGGCCCGGACGACTTCGt                                                                                                                                                        >  1:84511/1‑251 (MQ=255)
                                                                                                   aCTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGAGGTGGACCACGGTAAAGCGGCCAACGTCT‑TTGGTCGTCTTGGAGGGG‑CACAAAAAGGGCCGACCGGCCTCGTTgagga                                                                                                                                                  >  2:262642/1‑248 (MQ=255)
                                                                                                    cTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACt                                                                                                                                                >  1:78188/1‑251 (MQ=255)
                                                                                                      ggCCTCCAGTGCAGTTAACGCTGTTGAAGTTCGTGAAGCGGTGCAGTTCCTTACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGTAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAACAGGTCCGACCGACTTCGTTGAGAAc                                                                                                                                                 <  1:404513/248‑1 (MQ=255)
                                                                                                                       aCGCGGGTGAAGTTCGGGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGGTGAGATGCTGATCGCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAATCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTGGAGAACTACGCGAAGTATCTGCCGAc                                                                                                                             <  2:18499/251‑1 (MQ=255)
                                                                                                                           ggTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGATGATCTCCGGCAAACTGGCGGAAGATGACGCCGAA‑‑GCGCGCGAAATTGAAGGCGGGGCGGGACAAAGGGAAAGGGGCAGGAGACT‑TTGGTCGGATGGTGGAGG‑CCCAAAAAAGCACGACCGGCTTCGGTGGGAACTAcgggaggagtgtgccgaccacg                                                                                                                         >  1:43843/1‑231 (MQ=255)
                                                                                                                               gAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATTGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATg                                                                                                                      <  1:205754/250‑1 (MQ=255)
                                                                                                                                aaGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGAAAAAGGGAAAGCGGCAGAAGTCT‑TTGGTCGTATTGGAGCGG‑CAAAAAAAGGGCCGGCCGGCTTCGTTGAGAACTACGCGAAGGATTTGCCGACAACGATGcg                                                                                                                    >  2:278964/1‑250 (MQ=255)
                                                                                                                                 aGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGGTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGGCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTg                                                                                                                   <  2:395941/251‑1 (MQ=255)
                                                                                                                                            cGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGa                                                                                                                              <  1:411342/229‑1 (MQ=255)
                                                                                                                                            cGGTGAAGTTCCTGCCGGGTGAATATCGTCACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGGCAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTGTGGGAGCGG‑CACAAAAAGGACCGACCGACTTCGTTTCGAACTACGCGAAGAATCTGCCGa                                                                                                                              >  2:411342/1‑229 (MQ=255)
                                                                                                                                             ggTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGGCAACGGTAAAGAGGCAGAAGTCT‑TTGGGCGTATGGTAGCGG‑CACCAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGACGACAGCGAGGCTGACGAAAGCAGTc                                                                                                       >  1:434813/1‑251 (MQ=255)
                                                                                                                                             ggTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTc                                                                                                       <  1:209850/251‑1 (MQ=255)
                                                                                                                                             ggTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTc                                                                                                       <  2:81281/251‑1 (MQ=255)
                                                                                                                                                  aGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGccc                                                                                                                                                                    <  1:212995/185‑1 (MQ=255)
                                                                                                                                                      ccTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACCTCGTTGAGAACTCCGCGAAGTATCTGCCGACAGCGAAGCCGACGAAAACAGTCTATGCTGAt                                                                                              >  2:315810/1‑251 (MQ=255)
                                                                                                                                                               tGAATATCGTAAACCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGGTGGACAACGGGAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACTAAGCGAAGTATCTGCCGACAGCGATGCTGGCGAAAGCAGTCTATTCTGATACCGAAggg                                                                                     >  1:96303/1‑250 (MQ=255)
                                                                                                                                                                gAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGtt                                                                                    <  2:202643/251‑1 (MQ=255)
                                                                                                                                                                  atatCGTAACCCGCGTCGGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGAGGACGCCGAA‑‑GCGCGCGAAAGTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTTAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCGATGCTGATACCTAAGGtttt                                                                                  <  2:227687/251‑1 (MQ=255)
                                                                                                                                                                      cGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGGGATGCTGATCTCAGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGAAGCGG‑CACAAAAAGGCCCGACCGATTTCGTTGATAACTACGCGAAGTAGCTGACGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGtta                                                                              <  2:213133/251‑3 (MQ=255)
                                                                                                                                                                                gTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGATGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGGAAAGCGGCAGGAGTCT‑TTGGTCGTATGGGAGCGG‑CACAAAAAGGCCCGACCGACATCGTTGAGAACTAAGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGGCTATGCTGAAACCGAAGGTTTTGTCAGGGCAATggg                                                                    >  2:89550/1‑250 (MQ=255)
                                                                                                                                                                                  cTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGGAAAGCGGCAGAAGGCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCCGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGAt                                                                   >  2:353504/1‑250 (MQ=255)
                                                                                                                                                                                   tGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGGCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATAc                                                                 >  1:390282/1‑251 (MQ=255)
                                                                                                                                                                                    gTTTGATGGCACGATGGCGCTGTGCGGGGAGATGCTGATCTCCGGCAAACGGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCTGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑GTGGTCGTATGGCAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAGCTACGCGAATGATCTGCCGACAGCGAGGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGAt                                                                   <  1:250344/248‑1 (MQ=255)
                                                                                                                                                                                      ttGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTAGGCTGATACCGAAGGTTTTGTCAGTGAAATGGATAccc                                                               <  2:414500/250‑1 (MQ=255)
                                                                                                                                                                                       tGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCcg                                                              <  2:139178/250‑1 (MQ=255)
                                                                                                                                                                                               cGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CCCAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTggg                                                     >  1:207836/1‑251 (MQ=255)
                                                                                                                                                                                                   ggCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCCGAAGGCT‑TTGGTCGGATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGGATCTGCCGACCGCGATGCTGAAGAAAGAAGGCTATGCTGATACCGAAGGTTTTGTCAGTGAAATTGATACCCGCGCGCTGGGGATg                                                 >  1:91603/1‑251 (MQ=255)
                                                                                                                                                                                                     cgcTGTGCGTGGAGATGCGGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑GTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGGCAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGTATACCCGCGCGCTGGGGATGGc                                               <  1:376054/251‑1 (MQ=255)
                                                                                                                                                                                                             gTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTg                                       <  1:278964/251‑1 (MQ=255)
                                                                                                                                                                                                                        gATCTCCGGCAAACTGGCGAAAGATGACGCCAAA‑‑GCGCGCGACATTGCAGGCGGTGCTGGCCAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTCTAGAAGCGG‑CCCAACAAGGCCCGGACGACTTCGTTGGGAACTAAGCGCAGTATCCGGCGACAGGGATGCTGCCGAAAGCCGTCTATGCTGGTACCGAAGGTTTTGTCGGTGAAATGGCTACCCGCGCGCTGGGGATGGCAGTGGTTGCAAG‑Gggcggg                           >  1:385158/1‑250 (MQ=255)
                                                                                                                                                                                                                             ccGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAAT‑GGgcggcg                          <  2:43843/247‑1 (MQ=255)
                                                                                                                                                                                                                             ccGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAAT‑GGGCGGCGGAc                       <  2:239663/250‑1 (MQ=255)
                                                                                                                                                                                                                               ggCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGAGCTGGACAACGGTCAAGCGGCCGAAGTCT‑TTGGTCGTATGGTAGCGG‑CCCAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGGATCTGCCGACAGCGCTGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCCGTGAAATGGCTACCCGCGGGCCGGGGGTGGCCGTGGTTGCACT‑GGGGGGCGGACGcc                    >  2:30313/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aaaCTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGGTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGGGAACTACGCGAAGTATCTGCCGGCAGCGATGCTGGCGAAAGCAGTCTATGc                                                                                                  >  1:129994/1‑171 (MQ=255)
                                                                                                                                                                                                                                  aaaCTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGc                                                                                                  <  2:129994/171‑1 (MQ=255)
                                                                                                                                                                                                                                    aCTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAAT‑GGGCGGCGGACGCCGTCAg               <  2:78188/251‑1 (MQ=255)
                                                                                                                                                                                                                                       ggCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGGAAAGCGGCAGAAGGCT‑TTGGGCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACCTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGGCTATGCTGGTAACGAAGGGTTTGTCAGTGAAATGGAAACCCGCGCGCTGGGGGTGGCAGTGGTTGCAAT‑GGGGGGCGGagggcggcaggcc            >  2:2230/1‑239 (MQ=255)
                                                                                                                                                                                                                                            aaGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGGAGCGG‑CACAAAAAAGCCCGACCGACTTCGTTGAGAACTACGCGAAGGATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCCCGCGGTGGG‑ATGGCAGTGGGTGCAATTGggcgggggcgcccgccgggacctgg        >  1:218526/1‑230 (MQ=255)
                                                                                                                                                                                                                                                 gACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAACGCGGCAGAAGTCG‑TTGGTCGTATGGTAGCTG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGGCAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCTCGCTGGGGATGGCAGTGGTTGCAAG‑GGGCGGCGGACGCCGTCAGGCATCTGACACCa  <  2:73005/251‑1 (MQ=255)
                                                                                                                                                                                                                                                          ||                                                                                                                                                                                                                                                     
GACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAGCGG‑CACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAAT‑GGGCGGCGGACGCCGTCAGGCATCTGACACCA  >  NC_000913/4618864‑4619357

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑