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

GACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCA  >  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‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGt                                                                                                                                                                                            >  2:14783/1‑251 (MQ=255)
                                                                cgcACCACCGCGCTGCTCACCGACATGAATCAGGTACGGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAg                                                                                                                                                                                     <  2:216237/251‑1 (MQ=255)
                                                                 gcccgcccGCGCTGCTTACTGACATGAGTCAGGGACTGGCCTCGAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATTTCGTAACCCGCGTCTGTGTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCGCCGGCAAACTGGGGAGAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGATAACGGTAAAGCGGCATAAGTCTTTGGTCGTATGGTAGc                                                                                                                                                                                    <  2:292656/245‑1 (MQ=255)
                                                                   accaccGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCCGTTCCTGCCGGGTGAATATCGTACCCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGGAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTTGCCAACGGTAAAACGGCAGAAgctcgggcgagcacccgtctg                                                                                                                                                                                    >  1:421921/1‑230 (MQ=255)
                                                                                tCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGGAAAGCGGCAGAAGTCTTTGGGCGTAGGGGAGCGGGCAAAAAAGGccc                                                                                                                                                                    >  2:44567/1‑251 (MQ=255)
                                                                                      aCATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAAAGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCGCAAAAAGGCccgaccg                                                                                                                                                               >  1:100743/1‑250 (MQ=255)
                                                                                       cATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGAAAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCcgaccgac                                                                                                                                                             >  1:81281/1‑251 (MQ=255)
                                                                                         tGAATCAGTTCCTGGCCTCCAGTGCAGGGAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTACCTGACGGGTGAATATCGTAACCTGCGGCTGTTTTATGTCACGATGGCGCTGTGCGTGGAGATGCTGAACGCCGGCAAACTGGCGAAAGATGTCGCCGAA‑‑GCGCGCGAAATGGCAGGCGGTGCTGGACAACGGTAAGGCGTCATAAGTCTGTTGGCGTATGGCATCGGTACAAAAAGGGCCGTCCGACt                                                                                                                                                            <  2:303637/250‑1 (MQ=255)
                                                                                         tGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGACGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGGTTGATGTCACGATTGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGGTGGACAACGGGAAAGCGGGGGAAGACTTTGGGCGGTTGGTAGCGGGACAAAAAGGCCcgaacgacgt                                                                                                                                                           >  2:331217/1‑249 (MQ=255)
                                                                                         tGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAAAAGGCTTTGGGTGGATGGTAGCGGCACAAAAAGGCCCGACCTACtt                                                                                                                                                           >  1:25694/1‑251 (MQ=255)
                                                                                            aTCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGGAAAGCGGGAGAAGGCTTTGGTCGGATGGTAGAGGCACAAAAAGGCCCGGACGACTTCGt                                                                                                                                                        >  1:84511/1‑251 (MQ=255)
                                                                                                   aCTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGAGGTGGACCACGGTAAAGCGGCCAACGTCTTTGGTCGTCTTGGAGGGGCACAAAAAGGGCCGACCGGCCTCGTTgagga                                                                                                                                                  >  2:262642/1‑248 (MQ=255)
                                                                                                    cTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACt                                                                                                                                                >  1:78188/1‑251 (MQ=255)
                                                                                                      ggCCTCCAGTGCAGTTAACGCTGTTGAAGTTCGTGAAGCGGTGCAGTTCCTTACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGTAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAACAGGTCCGACCGACTTCGTTGAGAAc                                                                                                                                                 <  1:404513/248‑1 (MQ=255)
                                                                                                                       aCGCGGGTGAAGTTCGGGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGGTGAGATGCTGATCGCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAATCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTGGAGAACTACGCGAAGTATCTGCCGAc                                                                                                                             <  2:18499/251‑1 (MQ=255)
                                                                                                                           ggTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGATGATCTCCGGCAAACTGGCGGAAGATGACGCCGAA‑‑GCGCGCGAAATTGAAGGCGGGGCGGGACAAAGGGAAAGGGGCAGGAGACTTTGGTCGGATGGTGGAGGCCCAAAAAAGCACGACCGGCTTCGGTGGGAACTAcgggaggagtgtgccgaccacg                                                                                                                         >  1:43843/1‑231 (MQ=255)
                                                                                                                               gAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATTGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATg                                                                                                                      <  1:205754/250‑1 (MQ=255)
                                                                                                                                aaGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGAAAAAGGGAAAGCGGCAGAAGTCTTTGGTCGTATTGGAGCGGCAAAAAAAGGGCCGGCCGGCTTCGTTGAGAACTACGCGAAGGATTTGCCGACAACGATGcg                                                                                                                    >  2:278964/1‑250 (MQ=255)
                                                                                                                                 aGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGGTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGGCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTg                                                                                                                   <  2:395941/251‑1 (MQ=255)
                                                                                                                                            cGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGa                                                                                                                              <  1:411342/229‑1 (MQ=255)
                                                                                                                                            cGGTGAAGTTCCTGCCGGGTGAATATCGTCACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGGCAACGGTAAAGCGGCAGAAGTCTTTGGTCGTGTGGGAGCGGCACAAAAAGGACCGACCGACTTCGTTTCGAACTACGCGAAGAATCTGCCGa                                                                                                                              >  2:411342/1‑229 (MQ=255)
                                                                                                                                             ggTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGGCAACGGTAAAGAGGCAGAAGTCTTTGGGCGTATGGTAGCGGCACCAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGACGACAGCGAGGCTGACGAAAGCAGTc                                                                                                       >  1:434813/1‑251 (MQ=255)
                                                                                                                                             ggTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTc                                                                                                       <  1:209850/251‑1 (MQ=255)
                                                                                                                                             ggTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTc                                                                                                       <  2:81281/251‑1 (MQ=255)
                                                                                                                                                  aGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGccc                                                                                                                                                                    <  1:212995/185‑1 (MQ=255)
                                                                                                                                                      ccTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACCTCGTTGAGAACTCCGCGAAGTATCTGCCGACAGCGAAGCCGACGAAAACAGTCTATGCTGAt                                                                                              >  2:315810/1‑251 (MQ=255)
                                                                                                                                                               tGAATATCGTAAACCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGGTGGACAACGGGAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTAAGCGAAGTATCTGCCGACAGCGATGCTGGCGAAAGCAGTCTATTCTGATACCGAAggg                                                                                     >  1:96303/1‑250 (MQ=255)
                                                                                                                                                                gAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGtt                                                                                    <  2:202643/251‑1 (MQ=255)
                                                                                                                                                                  atatCGTAACCCGCGTCGGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGAGGACGCCGAA‑‑GCGCGCGAAAGTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTTAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCGATGCTGATACCTAAGGtttt                                                                                  <  2:227687/251‑1 (MQ=255)
                                                                                                                                                                      cGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGGGATGCTGATCTCAGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGAAGCGGCACAAAAAGGCCCGACCGATTTCGTTGATAACTACGCGAAGTAGCTGACGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGtta                                                                              <  2:213133/251‑3 (MQ=255)
                                                                                                                                                                                gTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGATGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGGAAAGCGGCAGGAGTCTTTGGTCGTATGGGAGCGGCACAAAAAGGCCCGACCGACATCGTTGAGAACTAAGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGGCTATGCTGAAACCGAAGGTTTTGTCAGGGCAATggg                                                                    >  2:89550/1‑250 (MQ=255)
                                                                                                                                                                                  cTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGGAAAGCGGCAGAAGGCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCCGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGAt                                                                   >  2:353504/1‑250 (MQ=255)
                                                                                                                                                                                   tGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGGCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATAc                                                                 >  1:390282/1‑251 (MQ=255)
                                                                                                                                                                                    gTTTGATGGCACGATGGCGCTGTGCGGGGAGATGCTGATCTCCGGCAAACGGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCTGTGCTGGACAACGGTAAAGCGGCAGAAGTCTGTGGTCGTATGGCAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAGCTACGCGAATGATCTGCCGACAGCGAGGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGAt                                                                   <  1:250344/248‑1 (MQ=255)
                                                                                                                                                                                      ttGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTAGGCTGATACCGAAGGTTTTGTCAGTGAAATGGATAccc                                                               <  2:414500/250‑1 (MQ=255)
                                                                                                                                                                                       tGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCcg                                                              <  2:139178/250‑1 (MQ=255)
                                                                                                                                                                                               cGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCCCAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTggg                                                     >  1:207836/1‑251 (MQ=255)
                                                                                                                                                                                                   ggCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCCGAAGGCTTTGGTCGGATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGGATCTGCCGACCGCGATGCTGAAGAAAGAAGGCTATGCTGATACCGAAGGTTTTGTCAGTGAAATTGATACCCGCGCGCTGGGGATg                                                 >  1:91603/1‑251 (MQ=255)
                                                                                                                                                                                                     cgcTGTGCGTGGAGATGCGGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTGTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGGCAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGTATACCCGCGCGCTGGGGATGGc                                               <  1:376054/251‑1 (MQ=255)
                                                                                                                                                                                                             gTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTg                                       <  1:278964/251‑1 (MQ=255)
                                                                                                                                                                                                                        gATCTCCGGCAAACTGGCGAAAGATGACGCCAAA‑‑GCGCGCGACATTGCAGGCGGTGCTGGCCAACGGTAAAGCGGCAGAAGTCTTTGGTCGTCTAGAAGCGGCCCAACAAGGCCCGGACGACTTCGTTGGGAACTAAGCGCAGTATCCGGCGACAGGGATGCTGCCGAAAGCCGTCTATGCTGGTACCGAAGGTTTTGTCGGTGAAATGGCTACCCGCGCGCTGGGGATGGCAGTGGTTGCAAGGggcggg                           >  1:385158/1‑250 (MQ=255)
                                                                                                                                                                                                                             ccGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGgcggcg                          <  2:43843/247‑1 (MQ=255)
                                                                                                                                                                                                                             ccGGCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGAc                       <  2:239663/250‑1 (MQ=255)
                                                                                                                                                                                                                               ggCAAACTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGAGCTGGACAACGGTCAAGCGGCCGAAGTCTTTGGTCGTATGGTAGCGGCCCAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGGATCTGCCGACAGCGCTGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCCGTGAAATGGCTACCCGCGGGCCGGGGGTGGCCGTGGTTGCACTGGGGGGCGGACGcc                    >  2:30313/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aaaCTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGGTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGGGAACTACGCGAAGTATCTGCCGGCAGCGATGCTGGCGAAAGCAGTCTATGc                                                                                                  >  1:129994/1‑171 (MQ=255)
                                                                                                                                                                                                                                  aaaCTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGc                                                                                                  <  2:129994/171‑1 (MQ=255)
                                                                                                                                                                                                                                    aCTGGCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAg               <  2:78188/251‑1 (MQ=255)
                                                                                                                                                                                                                                       ggCGAAAGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGGAAAGCGGCAGAAGGCTTTGGGCGTATGGTAGCGGCACAAAAAGGCCCGACCGACCTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGGCTATGCTGGTAACGAAGGGTTTGTCAGTGAAATGGAAACCCGCGCGCTGGGGGTGGCAGTGGTTGCAATGGGGGGCGGagggcggcaggcc            >  2:2230/1‑239 (MQ=255)
                                                                                                                                                                                                                                            aaGATGACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGGAGCGGCACAAAAAAGCCCGACCGACTTCGTTGAGAACTACGCGAAGGATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCCCGCGGTGGGATGGCAGTGGGTGCAATTGggcgggggcgcccgccgggacctgg        >  1:218526/1‑230 (MQ=255)
                                                                                                                                                                                                                                                 gACGCCGAA‑‑GCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAACGCGGCAGAAGTCGTTGGTCGTATGGTAGCTGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGGCAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCTCGCTGGGGATGGCAGTGGTTGCAAGGGGCGGCGGACGCCGTCAGGCATCTGACACCa  <  2:73005/251‑1 (MQ=255)
                                                                                                                                                                                                                                                          ||                                                                                                                                                                                                                                                     
GACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCA  >  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: