Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,853,601 A→T 5.4% intergenic (+12/‑94) PP_4263 → / → hemN cytochrome biogenesis protein/oxygen‑independent coproporphyrinogen III dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,853,6010AT 5.4% 873.4 / 26.5 259intergenic (+12/‑94)PP_4263/hemNcytochrome biogenesis protein/oxygen‑independent coproporphyrinogen III dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (121/123);  new base T (6/8);  total (128/131)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 270 total aligned reads displayed.

GCAACCGGGCTGGCTGCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTACGATCTCGCCGG  >  NC_002947/4853452‑4853750
                                                                                                                                                     |                                                                                                                                                     
gCAACCGGGCTGGCTGCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCaa                                                                                                                                                       >  1:5061611/1‑150 (MQ=255)
 cAACCGGGCTGGCTGCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCaaa                                                                                                                                                      <  2:6009640/150‑1 (MQ=255)
  aaCCGGGCTGGCTGCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGGTGCAAag                                                                                                                                                     >  1:2497270/1‑150 (MQ=255)
  aaCCGGGCTGGCTGCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAag                                                                                                                                                     >  2:4664177/1‑150 (MQ=255)
  aaCCGGGCTGGCTGCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAag                                                                                                                                                     <  1:3512252/150‑1 (MQ=255)
  aaCCGGGCTGGCTGCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGAGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGAAGGGCCAGTGACGGGGGCTGCAAag                                                                                                                                                     >  2:4869803/1‑150 (MQ=255)
    ccGGGCTGGCTGCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAagca                                                                                                                                                   >  1:4184310/1‑150 (MQ=255)
    ccGGGCTGGCTGCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAagca                                                                                                                                                   <  1:5229497/150‑1 (MQ=255)
    ccGGGCTGGCTGCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAagca                                                                                                                                                   <  2:1669327/150‑1 (MQ=255)
       ggctggctGCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCTGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGcc                                                                                                                                                >  1:3387400/1‑150 (MQ=255)
        gctggctgCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCTTTGCAGccc                                                                                                                                               <  1:2769701/150‑1 (MQ=255)
         ctggctgCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCTTTGCAGcccc                                                                                                                                              <  1:6226478/150‑1 (MQ=255)
            gctgCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCTTTGCAGcccccgt                                                                                                                                           <  2:4738426/150‑4 (MQ=255)
            gctgCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAAt                                                                                                                                           <  2:2573433/150‑1 (MQ=255)
              tgCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGc                                                                                                                                         >  1:5353931/1‑150 (MQ=255)
                 cGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAAAGCCCt                                                                                                                                      >  1:4653847/1‑150 (MQ=255)
                  gAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCtt                                                                                                                                     <  2:1432242/150‑1 (MQ=255)
                    aCGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGTTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCTTTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGa                                                                                                                                   <  1:2020328/150‑1 (MQ=255)
                    aCGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGa                                                                                                                                   <  2:1994623/150‑1 (MQ=255)
                      ggATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATa                                                                                                                                 >  1:1031602/1‑150 (MQ=255)
                      ggATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATa                                                                                                                                 <  2:4221040/150‑1 (MQ=255)
                          taaCGCCCTGTTGCGGCGACGCAGGGTGGGTGTGGCTGGGGGTGTGCTGTTGATGTTATTGGTTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTTATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACaaa                                                                                                                             <  1:6618855/149‑1 (MQ=255)
                            aCGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATc                                                                                                                           <  1:3865293/150‑1 (MQ=255)
                             cGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAAACa                                                                                                                          >  2:4904662/1‑150 (MQ=255)
                              gCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCCAAGCAGCCCCAATGCCCATGATACAAAACaa                                                                                                                         >  1:6397036/1‑150 (MQ=255)
                               cccTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCTTTGCAGcccc                                                                                                                                              <  2:5244358/128‑1 (MQ=255)
                               cccTGTGGGGGCGACGCTGCGGGCGAGTGGCTGGGGGTGTGGTGGTGTTGCTTTTGGGTATCTGGACCCTGCTGGGCCCACACCAGCAGGGGCTGATGGGCCATTGAGGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAac                                                                                                                        <  1:2443215/150‑1 (MQ=255)
                                ccTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAACCGCCCCAAAGCCCAAGATACaaaaaaaaa                                                                                                                       >  1:289920/1‑144 (MQ=255)
                                  tGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAAcaca                                                                                                                     <  2:3812344/150‑1 (MQ=255)
                                  tGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAAcaca                                                                                                                     <  1:1666590/150‑1 (MQ=255)
                                  tGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAAcaca                                                                                                                     <  2:747761/150‑1 (MQ=255)
                                    ttGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGt                                                                                                                   >  1:3316826/1‑150 (MQ=255)
                                       cggcgACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTc                                                                                                                <  2:4308019/150‑1 (MQ=255)
                                       cggcgACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTc                                                                                                                <  1:6015772/150‑1 (MQ=255)
                                        ggcgACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTcc                                                                                                               <  1:1075770/150‑1 (MQ=255)
                                            aCGCAGCGTGCGTGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCGGCCGGGCCCACACCAGCACTGGCTGATGGGCCTTTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACa                                                                                                           <  1:6660228/150‑1 (MQ=255)
                                              gCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGa                                                                                                         <  2:4972557/150‑1 (MQ=255)
                                                aGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAATCAGCCCCAAAGCCCTTGATACAAATAAACACAGTTTCCTACAGaag                                                                                                       >  1:4875762/1‑148 (MQ=255)
                                                aGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCGTGATACAAATCAACACACTTTCCTACAGACg                                                                                                       >  2:2341922/1‑150 (MQ=255)
                                                  cGTGCGAGGGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGTTGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGc                                                                                                     <  1:6332673/150‑1 (MQ=255)
                                                   gTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGcc                                                                                                    >  2:5271297/1‑150 (MQ=255)
                                                      cGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATa                                                                                                 <  1:6063295/150‑1 (MQ=255)
                                                      cGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATa                                                                                                 <  2:3423867/150‑1 (MQ=255)
                                                       gAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCATGATACAAATCAACACAGTTTCCTACAGACGGCCCTAg                                                                                                >  2:3506916/1‑150 (MQ=255)
                                                        aGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGa                                                                                               <  1:5862438/150‑1 (MQ=255)
                                                         gTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGAc                                                                                              <  1:5476244/150‑1 (MQ=255)
                                                            gCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCAACAGACGGCCATAGACTcc                                                                                           >  2:5825279/1‑150 (MQ=255)
                                                             cAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCg                                                                                          >  2:2257602/1‑150 (MQ=255)
                                                                gCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTg                                                                                       <  2:6838318/150‑1 (MQ=255)
                                                                         tGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCATGATAAAAAACAACACAGTTTACTAAAGAAGGCCCTAGACTCCGGTGCACtgctgc                                                                              >  1:5853875/1‑150 (MQ=255)
                                                                          ggTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTgctgct                                                                             <  2:3239767/150‑1 (MQ=255)
                                                                            tGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCGGGAAAGCAGCCCCAAAGCCCGTGTTACAAAACAACAAACTTTACTACAGACGGCCCAATACGCCGGTGCCCAGCTGctct                                                                           >  1:2285321/1‑150 (MQ=255)
                                                                              aTGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAAACAACACAGTTTCCTACAGACGGCCCTAGACTCCGGTGGACTGCTGCTCttt                                                                         >  2:1157914/1‑150 (MQ=255)
                                                                               tGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTc                                                                        <  1:468181/150‑1 (MQ=255)
                                                                                gCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCGGGAAAGCAGCCCCAAAGCCCTTGATACAAAAAAAAACAGATTACAAAAGAACGCCATAGAAAACGGTGCACAGCAGCTCgttac                                                                       >  1:355135/1‑148 (MQ=255)
                                                                                  tATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCgg                                                                     <  2:5083250/150‑1 (MQ=255)
                                                                                  tATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCGTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCgg                                                                     >  2:4182908/1‑150 (MQ=255)
                                                                                  tATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCGTGATACAAAAAAAAACACTTTCCTACAGACGGCCAAAGACACCGGTGCACTGCTGCTCGTTACgg                                                                     >  1:1400833/1‑150 (MQ=255)
                                                                                  tATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCGGCAAAGCAGCCCCAAAGCCCGTGATACAAAACAACACACTTTACTAAAGAAGGCCCTATACAACCGTGCACAGCTGCTCATTACgg                                                                     >  1:5194636/1‑150 (MQ=255)
                                                                                   aTTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCggg                                                                    <  2:1031602/150‑1 (MQ=255)
                                                                                   aTTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCGTTCCggg                                                                    >  2:3536906/1‑150 (MQ=255)
                                                                                    ttCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATAAACACAGTTTCCTAAAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCgggg                                                                   >  2:1807408/1‑150 (MQ=255)
                                                                                    ttCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCGGCAAAGCAGCCCCCAAGCCCGTGATAAAAAACAAAACACTTTGCTAAAGACGGCCCAAGGCACCCGTGCACAGCTGCTCGTTACgggg                                                                   >  1:6345115/1‑150 (MQ=255)
                                                                                     tCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGa                                                                  >  2:4432134/1‑150 (MQ=255)
                                                                                       ggTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCGGGAAAGCAGCCCCCAAGCCCATGATAAAAAAAAAAACAGGTTACAACAGACCGCCCAAGGCACCGGGGCCCAGCTGCTCGTTCCgggggcg                                                                >  1:8579/1‑147 (MQ=255)
                                                                                        tggtgtGGGCCCGGCAGGGTCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGGAAAGCAGCCCCAAAGCCCTTGATAAAAAAAAACACCGTTTCCTACAGACGGCCATAGACACCCGTGCACAGCTGCTCTTTCCGGGGAcga                                                               >  1:3655648/6‑148 (MQ=255)
                                                                                              ggACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATTCCCTTGATAAAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCa                                                         >  2:2952083/1‑150 (MQ=255)
                                                                                                  ccTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCTAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCt                                                     >  2:2781443/1‑150 (MQ=255)
                                                                                                  ccTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATAAACACAGTTTCCTACAGACGGCCATAGACGCCGGTGCACGGCTGCTCTTTCCGGGGACAACCCTCATGCt                                                     >  2:871518/1‑150 (MQ=255)
                                                                                                  ccTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAAACAACACAGTTTCCTACAGACGGCCCTAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACGACCCACAAGca                                                     >  2:6496926/1‑149 (MQ=255)
                                                                                                   cTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTc                                                    <  2:1582694/150‑1 (MQ=255)
                                                                                                   cTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTc                                                    <  2:4605560/150‑1 (MQ=255)
                                                                                                      ccGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTcgac                                                 <  2:5401285/150‑1 (MQ=255)
                                                                                                        gggCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTcgacga                                               <  2:6090844/150‑1 (MQ=255)
                                                                                                        gggCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCATGATACAAAAAAACACAGTTTACTACAGACGGCCAAAGACACCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCAAGCTcgacga                                               >  1:6156724/1‑150 (MQ=255)
                                                                                                            ccACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTACTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGCCTACCCTCATGCTCGACGACCTc                                           >  2:6367538/1‑150 (MQ=255)
                                                                                                              acacCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAAACAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCACAAGCGCGACGACCTCCg                                         >  2:928114/1‑150 (MQ=255)
                                                                                                                acCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGtt                                       <  2:3872901/150‑1 (MQ=255)
                                                                                                                 ccAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTg                                      <  2:3654260/150‑1 (MQ=255)
                                                                                                                  gggCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCCTAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTAACCTCAAGCTCGACGACCTCCGTTgg                                     >  2:5472180/3‑150 (MQ=255)
                                                                                                                  cAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTgg                                     <  2:3129743/150‑1 (MQ=255)
                                                                                                                  cAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTgg                                     <  2:3936773/150‑1 (MQ=255)
                                                                                                                  cAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAAACAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGAACTCCGGTgg                                     >  1:6776898/1‑150 (MQ=255)
                                                                                                                          gCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGa                             >  2:5619964/1‑150 (MQ=255)
                                                                                                                           cTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGAc                            <  2:2530557/150‑1 (MQ=255)
                                                                                                                            tGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCATGATACAAAAAAACACAGTTTCCTACAGACGGCCAAAGACTACCGTGCACTGCTGCTCTTTACGGGGACTAACCACATGCTCGACCACCACCCTTGGGATACCCAcc                           >  1:430727/1‑150 (MQ=255)
                                                                                                                               tGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGa                        <  1:1849179/150‑1 (MQ=255)
                                                                                                                                         attggggCTGCTTTGCAGCCCCCATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTa              >  1:3246689/4‑150 (MQ=255)
                                                                                                                                         aCGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTa              <  2:5165885/150‑1 (MQ=255)
                                                                                                                                          cGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTAc             >  2:182045/1‑150 (MQ=255)
                                                                                                                                           gggggCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTACg            >  2:3373081/1‑150 (MQ=255)
                                                                                                                                           gggggCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTACg            >  1:4790280/1‑150 (MQ=255)
                                                                                                                                           gggggCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACCGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTACg            >  1:5893149/1‑150 (MQ=255)
                                                                                                                                            ggggCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTACGa           >  2:3765707/1‑150 (MQ=255)
                                                                                                                                             gggCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTACGAt          <  1:4574639/150‑1 (MQ=255)
                                                                                                                                                 tGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTACGATCTCg      >  1:195873/1‑150 (MQ=255)
                                                                                                                                                 tGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTACGATCTCg      >  1:6872366/1‑150 (MQ=255)
                                                                                                                                                 tGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTACGATCTCg      >  2:6382144/1‑150 (MQ=255)
                                                                                                                                                  gCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTACGATCTCGc     >  1:3379869/1‑150 (MQ=255)
                                                                                                                                                    aaaGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTACGATCTCGCCg   <  1:5537065/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAACCGGGCTGGCTGCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTACGATCTCGCCGG  >  NC_002947/4853452‑4853750

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