Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,397,833 T→C 5.9% 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,8330TC 5.9% 228.5 / 4.2 68V141A (GTG→GCG) mutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base T (29/35);  new base C (2/2);  total (31/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.45e-01

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

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