Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,397,834 T→C 5.8% V141A (GTG→GCG)  mutL → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,397,8340TC 5.8% 231.2 / 4.1 69V141A (GTG→GCG) mutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base T (30/35);  new base C (2/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

AACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGATTTGCGGCACCGCTTTTCTTGAACAAGCGCTGGC  >  NC_000913/4397552‑4398089
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                               
aaCGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGCCGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCCGTgggtgc                                                                                                                                                                                                                                                           >  1:320739/1‑286 (MQ=255)
 aCGGGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTGTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGgac                                                                                                                                                                                                                                                           <  2:487734/288‑1 (MQ=255)
 aCGCGTTGGGGCAAAACTTAGCCGGATTGGTGATAACGGATGCGGTATCAAAAAAGGTGAGCTGGCGCGGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGgac                                                                                                                                                                                                                                                           <  1:424226/288‑1 (MQ=255)
 aCGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGgacg                                                                                                                                                                                                                                                          >  2:128359/1‑289 (MQ=255)
 aCGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGgacg                                                                                                                                                                                                                                                          >  1:302302/1‑289 (MQ=255)
 aCGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGgac                                                                                                                                                                                                                                                           <  2:337829/288‑1 (MQ=255)
 aCGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGgac                                                                                                                                                                                                                                                           <  2:230272/288‑1 (MQ=255)
         gggCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTgg                                                                                                                                                                                                                                                   <  1:470897/288‑1 (MQ=255)
              aaaCTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGc                                                                                                                                                                                                                                             <  1:298915/289‑1 (MQ=255)
                 cTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTgg                                                                                                                                                                                                                                          <  2:258362/289‑1 (MQ=255)
                   tATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGAt                                                                                                                                                                                                                                        <  2:302302/289‑1 (MQ=255)
                    aTCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGAt                                                                                                                                                                                                                                        >  2:310088/1‑288 (MQ=255)
                      ccactgtggatgtAACGGCGGCGGATCCAAAAAGGAGCGGCGGGGGCTGGGGGTGCGTCGTCATGCCACCAGTAAAGTCGCCTCTCTGGACGGTCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTg                                                                                                                                                                                                                                     <  1:155389/277‑1 (MQ=255)
                      ccGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTg                                                                                                                                                                                                                                     <  1:491558/289‑1 (MQ=255)
                              gTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTacaa                                                                                                                                                                                                                             <  2:320739/289‑1 (MQ=255)
                                         tGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCcggcg                                                                                                                                                                                                                   <  2:402688/288‑1 (MQ=255)
                                                                   gctggcgctggcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGaaa                                                                                                                                                                                        >  2:364876/1‑289 (MQ=255)
                                                                                 gTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACca                                                                                                                                                                          <  2:344406/289‑1 (MQ=255)
                                                                                         accaGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACca                                                                                                                                                                          >  1:184483/1‑281 (MQ=255)
                                                                                         accaGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACca                                                                                                                                                                          <  2:184483/281‑1 (MQ=255)
                                                                                                aaaTCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGatca                                                                                                                                                            >  2:129464/1‑288 (MQ=255)
                                                                                                            tGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCGGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGGACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTg                                                                                                                                                >  1:242176/1‑288 (MQ=255)
                                                                                                                   tctcGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGgac                                                                                                                                                                                                                                                           >  1:138391/1‑174 (MQ=255)
                                                                                                                   tctcGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGgac                                                                                                                                                                                                                                                           <  2:138391/174‑1 (MQ=255)
                                                                                                                    ctcGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGCCGCTCACTTCACGCACCGCAGAACAGCAGGACGCCTGGCAGGCCTATGCCGACGGGCGCGATCTGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGGGCTGTATCTGTTCTCTAACACCGCGGCGCGGCGCAAATTCCTGCGCCCCGAGAACACCCAATGTAACCACCTTGATGAGATCATCCGCCGCATTTCGCTGgcgg                                                                                                                                       >  1:76392/1‑288 (MQ=255)
                                                                                                                              ttATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGt                                                                                                                             <  2:183706/289‑1 (MQ=255)
                                                                                                                                agcAGCATGGGCTTCGCGGGTAGGGAGCTGGCGGGAGTCAGTTCGGTTGCCCGCCTGACCTTTACTTCACGCACCGCAGAACAGGAGGAAGCCTGGCAGGCCTAGGCGGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGGGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGAGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGt                                                                                                                             <  1:310088/285‑1 (MQ=255)
                                                                                                                                    gCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGACAACCGAATTTAACCACATTGATGAGATCAACCGCCGCATTTGGCTGGCCCGGTTCGACGTCACgc                                                                                                                        >  1:183706/1‑287 (MQ=255)
                                                                                                                                      cTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGAt                                                                                                                                                                                                                                        <  2:163838/174‑1 (MQ=255)
                                                                                                                                      cTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGAt                                                                                                                                                                                                                                        >  1:163838/1‑174 (MQ=255)
                                                                                                                                         ggCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCaa                                                                                                                    >  1:297097/1‑287 (MQ=255)
                                                                                                                                         ggCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCaa                                                                                                                    >  1:112149/1‑287 (MQ=255)
                                                                                                                                               cgcgGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCCCACTTCACGCACCGCAGAACCGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGCTCTGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGTGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCCCTGCGCACATTACTGCGCACCGAGAAAACCGACTCTAACCACCTTGCTGAACTCATCCGCCGCCTTGAGCTGGCGCGTTTCGACGTccagttccaccggtca                                                                                                             >  1:413629/1‑273 (MQ=255)
                                                                                                                                                   ggTGAGGCGCTGGCGAGTATCAGTTCGGGTTCCCGCCTGACGCTCACTTCACGCACCGAAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCAt                                                                                                          <  2:123272/288‑1 (MQ=255)
                                                                                                                                                     tGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGAAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATaa                                                                                                        <  2:426729/288‑1 (MQ=255)
                                                                                                                                                               gCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAACACCGAATTTAACCACATTGATGAGATCATCCGCCGCATCGCGCTGGCGCGTTTCGACGTCACGATCACTCTGCCGCACAACGGCAAATTtg                                                                                             >  1:486256/1‑289 (MQ=255)
                                                                                                                                                               gCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGCGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATtg                                                                                             >  1:363541/1‑289 (MQ=255)
                                                                                                                                                                cGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATtg                                                                                             <  2:112149/288‑1 (MQ=255)
                                                                                                                                                                cGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATtg                                                                                             <  2:344501/288‑1 (MQ=255)
                                                                                                                                                                cGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATtg                                                                                             <  2:2911/288‑1 (MQ=255)
                                                                                                                                                                           tCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTa                                                                                   <  1:128359/287‑1 (MQ=255)
                                                                                                                                                                           tCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTa                                                                                   <  2:297097/287‑1 (MQ=255)
                                                                                                                                                                           tCGGTTTCCCGCCTGACGCTACCTTCACGAACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTa                                                                                   <  1:382261/287‑1 (MQ=255)
                                                                                                                                                                              gTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACcgc                                                                               >  2:230251/1‑288 (MQ=255)
                                                                                                                                                                              gTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACcgc                                                                               >  2:382261/1‑288 (MQ=255)
                                                                                                                                                                               tttCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGAt                                                                                                                                                                                                                                        <  2:168836/134‑1 (MQ=255)
                                                                                                                                                                               tttCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGAt                                                                                                                                                                                                                                        >  1:168836/1‑134 (MQ=255)
                                                                                                                                                                               tttCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACcgcg                                                                              <  2:486256/288‑1 (MQ=255)
                                                                                                                                                                               tttCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGCGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACcgcg                                                                              <  2:363541/288‑1 (MQ=255)
                                                                                                                                                                                                            gCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATc                                                                                                                                                                                                                                       >  2:310073/1‑106 (MQ=255)
                                                                                                                                                                                                            gCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATc                                                                                                                                                                                                                                       <  1:310073/106‑1 (MQ=255)
                                                                                                                                                                                                            gCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGCGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCa                                                                            >  2:258296/1‑261 (MQ=255)
                                                                                                                                                                                                            gCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGCGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCa                                                                            <  1:258296/261‑1 (MQ=255)
                                                                                                                                                                                                                       ggAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCCGCGCATCCTGTGGGGACGACGCTGGAGGTGCGGGACCTGTTCTACAACAACCCGGCGCGGCGCAAAGTCCTGCGCACCGAGCACACCGAAGTTAATCACATCGATGAGATCATGCGCCGCATTGCGCTGGCGCGTTTCGGCGGCACGATCAACAGTGCGCACATCGGAACCACTGTGCTCTAGTACCGCGCACTGCCAGCAGGCGGGCAAAAAGAAAGgctgctcggaccc                                     >  1:258144/1‑278 (MQ=255)
                                                                                                                                                                                                                                    gCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGATTTGCGGCAcc                         >  2:279030/1‑288 (MQ=255)
                                                                                                                                                                                                                                         tGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGt                                                                                        <  1:369204/220‑1 (MQ=255)
                                                                                                                                                                                                                                         tGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGt                                                                                        >  2:369204/1‑220 (MQ=255)
                                                                                                                                                                                                                                              aaGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACcgcg                                                                              <  2:336193/225‑1 (MQ=255)
                                                                                                                                                                                                                                              aaGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACcgcg                                                                              >  1:336193/1‑225 (MQ=255)
                                                                                                                                                                                                                                              aaGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAg                                                                           <  1:385889/228‑1 (MQ=255)
                                                                                                                                                                                                                                              aaGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAg                                                                           >  2:385889/1‑228 (MQ=255)
                                                                                                                                                                                                                                              aaGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGATTTGCGGCACCGCTTTTCTTGa              >  2:56729/1‑289 (MQ=255)
                                                                                                                                                                                                                                                   cgcgATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACcgcg                                                                              <  2:11475/220‑1 (MQ=255)
                                                                                                                                                                                                                                                   cgcgATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACcgcg                                                                              >  1:11475/1‑220 (MQ=255)
                                                                                                                                                                                                                                                   cgcgATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGATTTGCGGCACCGCTTTTCTTGCACaa          >  2:150465/1‑288 (MQ=255)
                                                                                                                                                                                                                                                   cgcgATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGATTTGCGGCACCGCTTTTCTTGAACaa          >  2:231855/1‑288 (MQ=255)
                                                                                                                                                                                                                                                     cgATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGATTTGCGGCACCGCTTTTCTTGAACAAgc        >  2:389794/1‑288 (MQ=255)
                                                                                                                                                                                                                                                        tatGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATGGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACcgcg                                                                              <  2:467081/215‑1 (MQ=255)
                                                                                                                                                                                                                                                          tGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGATTTGCGGCACCGCTTTTCTTGAACAAGCGCTGGc  <  1:364876/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                               
AACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGATTTGCGGCACCGCTTTTCTTGAACAAGCGCTGGC  >  NC_000913/4397552‑4398089

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: