Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,626,391 A→G E160E (GAA→GAG radA → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,626,3910AG96.6% 174.9 / ‑2.4 60E160E (GAA→GAGradADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (24/33);  total (24/36)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.76e-01

CCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCT  >  NC_000913/4626147‑4626636
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
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)
                                                                                           cacgctgttgcggccaaCGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCCGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCCACGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTCATTGACTCGATCCAGGGGATGCATATGCCGGGTGGAGAGACATCGCCttgcc                                                                                                                                                      >  2:57779/16‑247 (MQ=255)
                                                                                              gctgttgcTGCAATCGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCg                                                                                                                                                   <  1:506932/244‑1 (MQ=255)
                                                                                                gtgtgcTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAATACGCTGTATGTCACCGGTGAAGCGTCGCTGCAACAGGTGGCAATGCGCTCTCATCGCCTTGGCCTGCCGACTGACAAGCTTAATAAGTTGTCGGAAACCAGCATCGAGCAGATCTGCCGGATTGCCGAAGAAGAGCAACCGAAGCTGAGGGTAATTGACTCGATCCAGGGGATGCATATGGCGGATGTACAGTCATCGccggtcagcgtg                                                                                                                                                 <  2:498435/246‑11 (MQ=255)
                                                                                                   tgcTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGgcg                                                                                                                                              >  1:124382/3‑251 (MQ=255)
                                                                                                   tgcTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGgc                                                                                                                                               <  1:104495/248‑1 (MQ=255)
                                                                                                               tgtgCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTAATCGCCTGGCAGCGTGGCGCAGGTGCGTGa                                                                                                                                   >  1:219252/1‑250 (MQ=255)
                                                                                                                     aaCTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGACACCAGCATCGAGCAGATCTGCCTGCTTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTAGATCCAGGTGATGCATATGGCGGATGTACAATCATCGCCGGGCAGCGTGGCGCAgtggcggggaacgggg                                                                                                                            >  2:482654/1‑236 (MQ=255)
                                                                                                                           cccAGCAGATGAAAACGCTGTATGTCACCCGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATGTGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGGTTGCCGAAGAAGAACAACCGACGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTTGCCGCGTGGCGCAGGGGAGGGAAACTGCGACTTAt                                                                                                                      >  1:129399/1‑251 (MQ=255)
                                                                                                                                cagATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGac                                                                                                                 <  1:393101/251‑1 (MQ=255)
                                                                                                                                 agATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAATCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAAcggcgggttttctgacc                                                                                                                >  2:404952/1‑240 (MQ=255)
                                                                                                                                  gATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCCACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGTCGGATGTACAGTCATCGCCTGGCAGCGTTGCGCGGGTGCGTGTAACGGCGGCttttctgacgc                                                                                                               >  2:73599/1‑242 (MQ=255)
                                                                                                                                     gAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCt                                                                                                            <  1:447006/251‑1 (MQ=255)
                                                                                                                                     gAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCAGGCAACGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCt                                                                                                            >  1:162638/1‑251 (MQ=255)
                                                                                                                                           gCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCACCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCa                                                                                                      <  1:368103/251‑1 (MQ=255)
                                                                                                                                              gTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCaa                                                                                                     <  2:269782/249‑1 (MQ=255)
                                                                                                                                                        ggCGAAGAGTCGCTGCACCAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGTATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGgtgt                                                                                         <  1:43163/251‑1 (MQ=255)
                                                                                                                                                          cGAAGAGTCGCTTCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGGCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTTACGCGCTTCGCCAAAACGCGCGGTGTgg                                                                                       <  2:164512/251‑1 (MQ=255)
                                                                                                                                                                  ggCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAAGCTCAATATGGTGTCGGAAACCAGCATCGAGCAGAGCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTc                                                                               <  2:256174/250‑1 (MQ=255)
                                                                                                                                                                  cGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTc                                                                               <  1:205671/251‑1 (MQ=255)
                                                                                                                                                                  cGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTc                                                                               <  1:39220/251‑1 (MQ=255)
                                                                                                                                                                      ggaaCAGGGGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCGAAATATGTTGGCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCAtgg                                                                           <  2:390621/249‑1 (MQ=255)
                                                                                                                                                                           aGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTg                                                                                                                                    <  2:319256/189‑1 (MQ=255)
                                                                                                                                                                           aGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTg                                                                                                                                    >  1:319256/1‑189 (MQ=255)
                                                                                                                                                                            ggTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCAtggttggac                                                                     >  1:450199/1‑246 (MQ=255)
                                                                                                                                                                                               cGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATggct                                                   <  1:88115/250‑3 (MQ=255)
                                                                                                                                                                                               cGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTAAAGTCATCGCCGGGCAGCGTGGCGCAGGTACGTGAAACGGCGGCTTATCTGACGCGCTTCGCCACAACGCGCGGGGGGGCGATTGTCAtgggtggacacgagatcgggaaagcac                                                   >  1:433846/1‑226 (MQ=255)
                                                                                                                                                                                               cGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTc                                                  <  2:143392/251‑1 (MQ=255)
                                                                                                                                                                                               cGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCAGCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATggct                                                   <  1:230845/250‑3 (MQ=255)
                                                                                                                                                                                               cGCCTTGGCCTGCCGACTGACAATCTCAATAGGTTGGCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGAATCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAATCGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACg                                                                  <  2:433846/235‑1 (MQ=255)
                                                                                                                                                                                                                 gACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGt                                <  1:66870/251‑1 (MQ=255)
                                                                                                                                                                                                                 gACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGaaa                                  <  1:141425/249‑1 (MQ=255)
                                                                                                                                                                                                                  aCAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGGCACGTAACCAAAACTGGCTCGCTGGCGGGCCCCAAAggg                               >  1:81246/1‑249 (MQ=255)
                                                                                                                                                                                                                         cAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGa                          <  1:411433/249‑1 (MQ=255)
                                                                                                                                                                                                                         cAATATGTAGTCGGAAACCAGCATCGAGCAGAGCTGCCTGATTGCCGAAGAAGAGCATCCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTTCAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGTCGATTGTCATGGTTGGACTCGTAACCAAAGATGGCTCGCTGGCGGGTCCTAAAGTGCTg                            <  1:56120/247‑1 (MQ=255)
                                                                                                                                                                                                                         cAATATGGTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTg                            <  1:418521/247‑1 (MQ=255)
                                                                                                                                                                                                                             atGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGCTCCAGGTGATGCATATGGCGGATGTACAGTCCTCGCCTTGCAGCGTGGCGCAGGTGCGTGAAACTGCGGCTTATATGACGCGCTTCGCCCAAACCCGCGGTGTGGCGCTTTTCATTGTTTGGCCCCTAACCCAAGATGGCTCCCTGGCGGGCCCTCAAGTGTTGGTACACtt                    >  2:38037/1‑250 (MQ=255)
                                                                                                                                                                                                                                 tgtCGGAAACCGGCATCGATCAGATCTGCCTCGTTACCGAAGAAGAGCAACCGAAGCTGATGGTAAGTGACTCTGTCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAATGTGCTGGAACACTGTAt                 <  2:217997/250‑1 (MQ=255)
                                                                                                                                                                                                                                               aTCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGAAGCGCTTCGCCCAAACGCGCGGTGTGGCGGTTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGTCGGGCCCTAAAGGGCTTGAAAACTTTATAGGCTGTTCGGTGCt  >  1:204899/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
CCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCT  >  NC_000913/4626147‑4626636

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