Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,626,247 A→G L112M (CTA→TTG radA → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,626,2470AG96.4% 163.6 / ‑4.1 56L112M (CTA→TTGradADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (25/29);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

CCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCT  >  NC_000913/4626009‑4626484
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              
ccATCCCCGAGTTGCGTCTTGCTGCGTCTCCAACGGTGGCGCTTACCGTGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACGCTCCGAGATCCGCCGTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTGCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCTGGGAAATCCACGCTGTTGCTGCCAACGCtg                                                                                                                                                                                                                                    <  2:369179/251‑1 (MQ=255)
      ccGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCaaa                                                                                                                                                                                                                              <  1:373779/251‑1 (MQ=255)
                   tgctgcGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGAGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCTGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCcagc                                                                                                                                                                                                                    >  1:414208/1‑248 (MQ=255)
                      tgcGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCcagcagctg                                                                                                                                                                                                               >  2:4909/1‑247 (MQ=255)
                       gcGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCAGCAATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACGGGCCCAGCAGATGa                                                                                                                                                                                                              >  2:141425/1‑250 (MQ=255)
                       gcGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTTTGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTTCCGCGTTTTTCCACCTTATTTAAAGATTTCGACCGCGTACTAGGCTGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGaa                                                                                                                                                                                                             <  1:225868/251‑1 (MQ=255)
                       gcGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGa                                                                                                                                                                                                              <  1:8155/250‑1 (MQ=255)
                                  ggtggCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAATTGCCATTCTGATTGGCGGTAACCCTGGGGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGaaacagctggatg                                                                                                                                                                                                  >  2:192747/1‑241 (MQ=255)
                                     ggggcgTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATAACAGCCTTGAAGAGCGGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCa                                                                                                                                                                                               <  2:46548/248‑1 (MQ=255)
                                                  cGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTc                                                                                                                                                                                  <  1:183484/251‑1 (MQ=255)
                                                  cGTCTCAGCGGATATTCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTACAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCCGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATACACGCTGTTGCTGCACACGCTGTGCAACATGGCCCAGCAGCTGAAAACCCTGTaagacacagcccacgagtc                                                                                                                                                                                  >  2:25736/1‑233 (MQ=255)
                                                        aGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGc                                                                                                                                                                             >  1:256174/1‑250 (MQ=255)
                                                           ggCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCAACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCgaggcaac                                                                                                                                                                          >  2:88115/1‑243 (MQ=255)
                                                                    ggTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCTTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAAt                                                                                                                                                                <  1:9054/251‑1 (MQ=255)
                                                                            cGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCCGATGAAACCGCTGTATGTCACCGGCGGCGAGTCGCTGAAACAGGTGGCAATGCGCGCTc                                                                                                                                                        >  2:43163/1‑251 (MQ=255)
                                                                            cGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATgctccctc                                                                                                                                                        >  2:164641/1‑245 (MQ=255)
                                                                                     aataggccAGAAACGCGCCGAGAGCAGCCGTGAAGAGCTGCCGCGTGTTTCCACCGGATTTAAAGAGTTCGACCGCGGTCTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATGGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGccgtg                                                                                                                                               <  2:240827/245‑4 (MQ=255)
                                                                                     aaaaGTCCAGAAACTCTCCGATATCCGCCTTGAAGAGCTGCCTCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGAGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCACAAGCTGTGCCAACTGGCCCCGCAGATGAAAACGCGGTATGTCCCCGGCGAAGAGTCGCCGCAACAGGTGGCCATGCGCGCACAGCGCCTTg                                                                                                                                               >  2:220887/1‑251 (MQ=255)
                                                                                       aaGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAAACCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAACACGCTGTATGTCACCGGCGCAGAGTAGCTGCACAAGGTGGCAACGCGCGCTCATCGCCTTgac                                                                                                                                             >  2:14916/1‑249 (MQ=255)
                                                                                            cAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGcc                                                                                                                                        >  1:283149/1‑251 (MQ=255)
                                                                                                 aCTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTg                                                                                                                                   <  2:418821/251‑1 (MQ=255)
                                                                                                  ctctCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGAGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATACTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGATGCTGAAAACGCTGTGCAACAGGGCCCAGCAGGTGAAAACGATGTATGTCCCCGGCGAAAAGCCGCTGCCACAGGTGGCAATACACGCTCCTCGCCGTgggcacccgaacg                                                                                                                                   >  2:275642/1‑239 (MQ=255)
                                                                                                  ctctCCGATATCAGCCTTGAAGAGCTGCAGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACACTGGTGCGGGGAAATCCACGCTGTTGCTGCACAAGCTGTGCAAACTGGCCCATCAGATGAAAACGCTGTATGTAACCGGCGGAGAGGCGCTGCAACAGGTGGCAATTCGACCTCATCGCCTTGGCCCGACGACTgg                                                                                                                                  >  2:37313/1‑250 (MQ=255)
                                                                                                   cctcCGATAGCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGAc                                                                                                                                 <  2:250006/250‑1 (MQ=255)
                                                                                                           atCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCa                                                                                                                          >  1:12537/1‑250 (MQ=255)
                                                                                                             cAGCCTTGAAGAGCCGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAata                                                                                                                       <  1:489067/251‑1 (MQ=255)
                                                                                                                        agCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAACGGATGAAAACGCTGTATGTCACCGGCGAAGAATCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGa                                                                                                             >  1:143392/1‑250 (MQ=255)
                                                                                                                                cgTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCGTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCAt                                                                                                    <  2:224043/251‑1 (MQ=255)
                                                                                                                                 gTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGGCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAAGGTCAATATGTTGTCGGAAACCAGCATc                                                                                                   <  1:186907/251‑1 (MQ=255)
                                                                                                                                    tttCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAACCCAGCATCGa                                                                                                 <  2:61006/250‑1 (MQ=255)
                                                                                                                                     ttCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGag                                                                                                <  1:386978/250‑2 (MQ=255)
                                                                                                                                     ttCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGACGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGACTGCCGACTGCCAATCTCAATATTGTGTCGGAAACCAGCAtagggg                                                                                               >  2:270961/1‑246 (MQ=255)
                                                                                                                                      cacactggATTTAAAGAGTTTGACCGCGTACTAGGCGTCAGCGTGGTGTCAGGTAGGGCCAGTCTGAGTGGCGGGAACCCTGGTGCGGGGAAAGCCACGCTGTTGCTGTAAACGCTGGGCAAACTGGCCCAGCAGGTGAAAACGCTGTAGGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCCGCATCGagca                                                                                              <  2:178804/245‑4 (MQ=255)
                                                                                                                                       ccACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGagca                                                                                              <  2:123644/250‑4 (MQ=255)
                                                                                                                                          ccGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATc                                                                                          <  1:188423/251‑1 (MQ=255)
                                                                                                                                           cGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCt                                                                                         <  1:339437/251‑1 (MQ=255)
                                                                                                                                                               gcgggtgatgcggtggcgTGGAGCTAGGAATGGCCATCTTGATGGGCCGTAGCCCTGGTGCGGGGAAATCCCCTCTGTTGCTGCAAACGCGTTGCAAACTTGCGCAGCAGATGAAATCGCTCTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCTCTCATCGCCTTGGCCTGCCGACTTACAATCTCAATATGTTGTCGGAAACCATCATCGAGCAGATCTGCCTGATTGctgcataagag                                                                     <  2:220057/242‑11 (MQ=255)
                                                                                                                                                                   aCTAGGCGGCGGCGTGGTGCCAGGAAGGGCCATTCTTATTGGCGGTAACCCTGGGGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCAGCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAAc                                                                <  2:392836/251‑1 (MQ=255)
                                                                                                                                                                    cTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGGCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATAGACTTTGCCTTCCGGCTGACAATCTCCATATGGTGGCGGAAACCAGCATCGTGCCGATCTTGCTGATTGCCGaagaaggacaaac                                                                >  2:415845/1‑244 (MQ=255)
                                                                                                                                                                     tAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCg                                                              >  1:211470/1‑251 (MQ=255)
                                                                                                                                                                       ggcggcggcGTGGTGCCAGGAAGTGCCATTCTGATTGGCGTTAACCCTGGGGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGaa                                                            <  2:360676/251‑1 (MQ=255)
                                                                                                                                                                       ggcggcggcGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGaa                                                            <  1:164641/251‑1 (MQ=255)
                                                                                                                                                                       ggcggcggcGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCCGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGaa                                                            <  1:470635/251‑1 (MQ=255)
                                                                                                                                                                          ggcggcgTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCt                                                         >  1:269782/1‑251 (MQ=255)
                                                                                                                                                                                       cAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCACACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGGAGAGTCGCTGCAACAGGTTGCAATGCGCGCTCATCGCCTTGGCCTGCCGGCTGACAATCTCCATATGTTTTCGGAAACCAGCATAGAGCAGATCTGCCTTATTTCCGAAGAAAAGCAACCGAATCTACTGGTAATTGAc                                            >  2:870/1‑251 (MQ=255)
                                                                                                                                                                                        aGGAAGTGCCATTCTGATTGGCGGTAAACCTGGTACGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGCACACGCTGTATGTCACCAGCGAAGAGGCGCTGCCACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTTACGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTTCCTGATTTCCGAAGGAGCGCAACCGGAGCTGATGGACATTGACt                                           >  2:39220/1‑251 (MQ=255)
                                                                                                                                                                                          gAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATTGTAATTGACTCg                                         <  2:54653/251‑1 (MQ=255)
                                                                                                                                                                                             gTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGCGCAGATCTGCCTGATTGCCGAGGAAGAGCCACCGAAGATGATGGTAATTGACTCGATc                                      >  2:206336/1‑251 (MQ=255)
                                                                                                                                                                                                          tGGCGGTCACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACACTGTGCAACCTGGCCCAGCAGATGAAACCGCTGTATGTCACCGGCGAAGAGTCCCTGCAACCGGAGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCACTTTGCTGTCGGAAACCCACATCGAGCCGATCTGCCTGTTTGCAGCAGAAAGTCACCCGACGCTGCTGGTAATTGACGCGAGCCAGTTGCTGCat                          >  2:368103/1‑250 (MQ=255)
                                                                                                                                                                                                            gCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATg                       <  1:21878/251‑1 (MQ=255)
                                                                                                                                                                                                              ggTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGc                     <  1:128097/251‑1 (MQ=255)
                                                                                                                                                                                                                   cccTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATg                <  1:221977/251‑1 (MQ=255)
                                                                                                                                                                                                                       ggTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGCAACCCGCATCGAGAAGATCTGCCTGATTGCCGAAGAAGCACAACAGACGCTGATGGGAATTGACTCGATCCAGGAGAGGCATATGGCgggtggaa             >  2:446124/1‑244 (MQ=255)
                                                                                                                                                                                                                       ggTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACa            <  1:4909/251‑1 (MQ=255)
                                                                                                                                                                                                                        gTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGCTGTAccg           >  2:205671/1‑249 (MQ=255)
                                                                                                                                                                                                                                 aaTCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCCTCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGccc  >  1:71731/1‑250 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              
CCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCT  >  NC_000913/4626009‑4626484

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: