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→GAG) radADNA 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

CCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCT  >  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)
                         aCTAGGCGGCGGCGTGGTGCCAGGAAGGGCCATTCTTATTGGCGGTAACCCTGGGGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCAGCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAAc                                                                                                                                                                                                                        <  2:392836/251‑1 (MQ=255)
                          cTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGGCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATAGACTTTGCCTTCCGGCTGACAATCTCCATATGGTGGCGGAAACCAGCATCGTGCCGATCTTGCTGATTGCCGaagaaggacaaac                                                                                                                                                                                                                        >  2:415845/1‑244 (MQ=255)
                           tAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCg                                                                                                                                                                                                                      >  1:211470/1‑251 (MQ=255)
                             ggcggcggcGTGGTGCCAGGAAGTGCCATTCTGATTGGCGTTAACCCTGGGGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGaa                                                                                                                                                                                                                    <  2:360676/251‑1 (MQ=255)
                             ggcggcggcGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGaa                                                                                                                                                                                                                    <  1:164641/251‑1 (MQ=255)
                             ggcggcggcGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCCGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGaa                                                                                                                                                                                                                    <  1:470635/251‑1 (MQ=255)
                                ggcggcgTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCt                                                                                                                                                                                                                 >  1:269782/1‑251 (MQ=255)
                                             cAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCACACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGGAGAGTCGCTGCAACAGGTTGCAATGCGCGCTCATCGCCTTGGCCTGCCGGCTGACAATCTCCATATGTTTTCGGAAACCAGCATAGAGCAGATCTGCCTTATTTCCGAAGAAAAG‑CAACCGAATCTACTGGTAATTGAc                                                                                                                                                                                                    >  2:870/1‑251 (MQ=255)
                                              aGGAAGTGCCATTCTGATTGGCGGTAAACCTGGTACGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGCACACGCTGTATGTCACCAGCGAAGAGGCGCTGCCACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTTACGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTTCCTGATTTCCGAAGGAGCG‑CAACCGGAGCTGATGGACATTGACt                                                                                                                                                                                                   >  2:39220/1‑251 (MQ=255)
                                                gAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATTGTAATTGACTCg                                                                                                                                                                                                 <  2:54653/251‑1 (MQ=255)
                                                   gTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGCGCAGATCTGCCTGATTGCCGAGGAAGAG‑CCACCGAAGATGATGGTAATTGACTCGATc                                                                                                                                                                                              >  2:206336/1‑251 (MQ=255)
                                                                tGGCGGTCACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACACTGTGCAACCTGGCCCAGCAGATGAAACCGCTGTATGTCACCGGCGAAGAGTCCCTGCAACCGGAGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCACTTTGCTGTCGGAAACCCACATCGAGCCGATCTGCCTGTTTGCAGCAGAA‑AGTCACCCGACGCTGCTGGTAATTGACGCGAGCCAGTTGCTGCat                                                                                                                                                                                  >  2:368103/1‑250 (MQ=255)
                                                                  gCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATg                                                                                                                                                                               <  1:21878/251‑1 (MQ=255)
                                                                    ggTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGc                                                                                                                                                                             <  1:128097/251‑1 (MQ=255)
                                                                         cccTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATg                                                                                                                                                                        <  1:221977/251‑1 (MQ=255)
                                                                             ggTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGCAACCCGCATCGAGAAGATCTGCCTGATTGCCGAAGAAGCA‑CAACAGACGCTGATGGGAATTGACTCGATCCAGGAGAGGCATATGGCgggtggaa                                                                                                                                                                     >  2:446124/1‑244 (MQ=255)
                                                                             ggTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACa                                                                                                                                                                    <  1:4909/251‑1 (MQ=255)
                                                                              gTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGCTGTAccg                                                                                                                                                                   >  2:205671/1‑249 (MQ=255)
                                                                                       aaTCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCCTCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGccc                                                                                                                                                          >  1:71731/1‑250 (MQ=255)
                                                                                           cacgctgttgcggccaaCGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCCGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCCACGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTCATTGACTCGATCCAGGGGATGCATATGCCGGGTGGAGAGACATCGCCttgcc                                                                                                                                                      >  2:57779/16‑247 (MQ=255)
                                                                                              gctgttgcTGCAATCGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCg                                                                                                                                                   <  1:506932/244‑1 (MQ=255)
                                                                                                gtgtgcTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAATACGCTGTATGTCACCGGTGAAGCGTCGCTGCAACAGGTGGCAATGCGCTCTCATCGCCTTGGCCTGCCGACTGACAAGCTTAATAAGTTGTCGGAAACCAGCATCGAGCAGATCTGCCGGATTGCCGAAGAAGAG‑CAACCGAAGCTGAGGGTAATTGACTCGATCCAGGGGATGCATATGGCGGATGTACAGTCATCGccggtcagcgtg                                                                                                                                                 <  2:498435/246‑11 (MQ=255)
                                                                                                   tgcTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGgcg                                                                                                                                              >  1:124382/3‑251 (MQ=255)
                                                                                                   tgcTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGgc                                                                                                                                               <  1:104495/248‑1 (MQ=255)
                                                                                                               tgtgCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTAATCGCCTGGCAGCGTGGCGCAGGTGCGTGa                                                                                                                                   >  1:219252/1‑250 (MQ=255)
                                                                                                                     aaCTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGACACCAGCATCGAGCAGATCTGCCTGCTTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTAGATCCAGGTGATGCATATGGCGGATGTACAATCATCGCCGGGCAGCGTGGCGCAgtggcggggaacgggg                                                                                                                            >  2:482654/1‑236 (MQ=255)
                                                                                                                           cccAGCAGATGAAAACGCTGTATGTCACCCGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATGTGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGGTTGCCGAAGAAGAA‑CAACCGACGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTTGCCGCGTGGCGCAGGGGAGGGAAACTGCGACTTAt                                                                                                                      >  1:129399/1‑251 (MQ=255)
                                                                                                                                cagATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGac                                                                                                                 <  1:393101/251‑1 (MQ=255)
                                                                                                                                 agATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAATCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAAcggcgggttttctgacc                                                                                                                >  2:404952/1‑240 (MQ=255)
                                                                                                                                  gATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCCACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGTCGGATGTACAGTCATCGCCTGGCAGCGTTGCGCGGGTGCGTGTAACGGCGGCttttctgacgc                                                                                                               >  2:73599/1‑242 (MQ=255)
                                                                                                                                     gAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCt                                                                                                            <  1:447006/251‑1 (MQ=255)
                                                                                                                                     gAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCAGGCAACGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCt                                                                                                            >  1:162638/1‑251 (MQ=255)
                                                                                                                                           gCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCACCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCa                                                                                                      <  1:368103/251‑1 (MQ=255)
                                                                                                                                              gTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCaa                                                                                                     <  2:269782/249‑1 (MQ=255)
                                                                                                                                                        ggCGAAGAGTCGCTGCACCAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGTATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGgtgt                                                                                         <  1:43163/251‑1 (MQ=255)
                                                                                                                                                          cGAAGAGTCGCTTCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGGCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTTACGCGCTTCGCCAAAACGCGCGGTGTgg                                                                                       <  2:164512/251‑1 (MQ=255)
                                                                                                                                                                  ggCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAAGCTCAATATGGTGTCGGAAACCAGCATCGAGCAGAGCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTc                                                                               <  2:256174/250‑1 (MQ=255)
                                                                                                                                                                  cGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTc                                                                               <  1:205671/251‑1 (MQ=255)
                                                                                                                                                                  cGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTc                                                                               <  1:39220/251‑1 (MQ=255)
                                                                                                                                                                      ggaaCAGGGGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCGAAATATGTTGGCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCAtgg                                                                           <  2:390621/249‑1 (MQ=255)
                                                                                                                                                                           aGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTg                                                                                                                                    <  2:319256/189‑1 (MQ=255)
                                                                                                                                                                           aGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTg                                                                                                                                    >  1:319256/1‑189 (MQ=255)
                                                                                                                                                                            ggTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCAtggttggac                                                                     >  1:450199/1‑246 (MQ=255)
                                                                                                                                                                                               cGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATggct                                                   <  1:88115/250‑3 (MQ=255)
                                                                                                                                                                                               cGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTAAAGTCATCGCCGGGCAGCGTGGCGCAGGTACGTGAAACGGCGGCTTATCTGACGCGCTTCGCCACAACGCGCGGGGGGGCGATTGTCAtgggtggacacgagatcgggaaagcac                                                   >  1:433846/1‑226 (MQ=255)
                                                                                                                                                                                               cGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTc                                                  <  2:143392/251‑1 (MQ=255)
                                                                                                                                                                                               cGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCAGCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATggct                                                   <  1:230845/250‑3 (MQ=255)
                                                                                                                                                                                               cGCCTTGGCCTGCCGACTGACAATCTCAATAGGTTGGCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGAATCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAATCGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACg                                                                  <  2:433846/235‑1 (MQ=255)
                                                                                                                                                                                                                 gACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGt                                <  1:66870/251‑1 (MQ=255)
                                                                                                                                                                                                                 gACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGaaa                                  <  1:141425/249‑1 (MQ=255)
                                                                                                                                                                                                                  aCAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGGCACGTAACCAAAACTGGCTCGCTGGCGGGCCCCAAAggg                               >  1:81246/1‑249 (MQ=255)
                                                                                                                                                                                                                         cAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGa                          <  1:411433/249‑1 (MQ=255)
                                                                                                                                                                                                                         cAATATGTAGTCGGAAACCAGCATCGAGCAGAGCTGCCTGATTGCCGAAGAAGAG‑CATCCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTTCAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGTCGATTGTCATGGTTGGACTCGTAACCAAAGATGGCTCGCTGGCGGGTCCTAAAGTGCTg                            <  1:56120/247‑1 (MQ=255)
                                                                                                                                                                                                                         cAATATGGTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTg                            <  1:418521/247‑1 (MQ=255)
                                                                                                                                                                                                                             atGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGCTCCAGGTGATGCATATGGCGGATGTACAGTCCTCGCCTTGCAGCGTGGCGCAGGTGCGTGAAACTGCGGCTTATATGACGCGCTTCGCCCAAACCCGCGGTGTGGCGCTTTTCATTGTTTGGCCCCTAACCCAAGATGGCTCCCTGGCGGGCCCTCAAGTGTTGGTACACtt                    >  2:38037/1‑250 (MQ=255)
                                                                                                                                                                                                                                 tgtCGGAAACCGGCATCGATCAGATCTGCCTCGTTACCGAAGAAGAG‑CAACCGAAGCTGATGGTAAGTGACTCTGTCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAATGTGCTGGAACACTGTAt                 <  2:217997/250‑1 (MQ=255)
                                                                                                                                                                                                                                               aTCGAGCAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGAAGCGCTTCGCCCAAACGCGCGGTGTGGCGGTTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGTCGGGCCCTAAAGGGCTTGAAAACTTTATAGGCTGTTCGGTGCt  >  1:204899/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
CCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAG‑CAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCT  >  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: ‑