Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,853,602 A→T 5.8% intergenic (+13/‑93) 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,6020AT 5.8% 814.7 / 27.5 250intergenic (+13/‑93)PP_4263/hemNcytochrome biogenesis protein/oxygen‑independent coproporphyrinogen III dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (117/118);  new base T (6/9);  total (123/127)
Fisher's exact test for biased strand distribution p-value = 5.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

Only 100 of 266 total aligned reads displayed.

CAACCGGGCTGGCTGCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTACGATCTCGCCGG  >  NC_002947/4853453‑4853750
                                                                                                                                                     |                                                                                                                                                    
cAACCGGGCTGGCTGCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCaaa                                                                                                                                                      <  1:518514/150‑1 (MQ=255)
 aaCCGGGCTGGCTGCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGGTGCAAag                                                                                                                                                     >  1:2497270/1‑150 (MQ=255)
 aaCCGGGCTGGCTGCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGAGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGAAGGGCCAGTGACGGGGGCTGCAAag                                                                                                                                                     >  2:4869803/1‑150 (MQ=255)
      ggctggctGCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCTGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGcc                                                                                                                                                >  1:3387400/1‑150 (MQ=255)
       gctggctgCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCTTTGCAGccc                                                                                                                                               <  1:5759303/150‑1 (MQ=255)
       gctggctgCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCTTTGCAGccc                                                                                                                                               <  1:2769701/150‑1 (MQ=255)
         tggctgCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGGGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCTTTGCAGccccc                                                                                                                                             <  1:655677/150‑2 (MQ=255)
           gctgCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCTTTGCAGcccccgt                                                                                                                                           <  2:4752726/150‑4 (MQ=255)
           gctgCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAAt                                                                                                                                           <  2:2572902/150‑1 (MQ=255)
                cGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCt                                                                                                                                      >  1:6223829/1‑150 (MQ=255)
                 gAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCtt                                                                                                                                     >  2:3482046/1‑150 (MQ=255)
                 gAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCtt                                                                                                                                     >  1:464132/1‑150 (MQ=255)
                  aaCGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTg                                                                                                                                    >  2:4226728/1‑150 (MQ=255)
                  aaCGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAACCAGCCCCAATGcccgtt                                                                                                                                    >  1:3525091/1‑147 (MQ=255)
                   aCGGATGACCGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGAGGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGa                                                                                                                                   <  1:4127007/150‑1 (MQ=255)
                   aCGGACGCACGCCCTGTTGCGGCGCCGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCGCACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGa                                                                                                                                   <  2:1963665/150‑1 (MQ=255)
                     ggATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATa                                                                                                                                 >  1:1031602/1‑150 (MQ=255)
                      gATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATAc                                                                                                                                <  2:5440882/150‑1 (MQ=255)
                      gATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCATGATAc                                                                                                                                >  2:2904652/1‑150 (MQ=255)
                         gAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACaaa                                                                                                                             <  2:4539012/150‑1 (MQ=255)
                         gAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCATGATACaaa                                                                                                                             >  1:5967531/1‑150 (MQ=255)
                          aaCGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCCATGCCCGTGATACAAAt                                                                                                                            >  1:5262069/1‑150 (MQ=255)
                          aaCGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAAt                                                                                                                            >  1:6245969/1‑150 (MQ=255)
                             gCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCCAAGCAGCCCCAATGCCCATGATACAAAACaa                                                                                                                         >  1:6397036/1‑150 (MQ=255)
                             gCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCaa                                                                                                                         <  2:517038/150‑1 (MQ=255)
                                   ttGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGt                                                                                                                   >  1:3316826/1‑150 (MQ=255)
                                       ggcgACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTcc                                                                                                               <  1:1075770/150‑1 (MQ=255)
                                        gcgACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCGGCAAAGCAGCCCCAAATCCCGTGATACAAAAAAACACAGTTTCCt                                                                                                              >  1:5880380/1‑150 (MQ=255)
                                         cgACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTa                                                                                                             <  2:4988922/150‑1 (MQ=255)
                                           aCGCAGCGTGCGTGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCGGCCGGGCCCACACCAGCACTGGCTGATGGGCCTTTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACa                                                                                                           <  1:6660228/150‑1 (MQ=255)
                                             gCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGa                                                                                                         <  1:4699872/150‑1 (MQ=255)
                                               aGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCATGATACAAAGCCACACAGTTTCCTACAgccc                                                                                                       >  2:2341516/1‑147 (MQ=255)
                                                 cGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGc                                                                                                     <  1:6332919/150‑1 (MQ=255)
                                                  gTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGcc                                                                                                    >  2:1808785/1‑150 (MQ=255)
                                                  gTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGcc                                                                                                    >  2:5271297/1‑150 (MQ=255)
                                                   tGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCa                                                                                                   <  1:6811732/150‑1 (MQ=255)
                                                     cGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATa                                                                                                 <  1:6063231/150‑1 (MQ=255)
                                                          ggCAGGGGGTGTGCTGGTGATGCTATTCGGTATCTGGCCCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTc                                                                                            <  1:2458788/150‑1 (MQ=255)
                                                           gCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCAACAGACGGCCATAGACTcc                                                                                           >  2:5825279/1‑150 (MQ=255)
                                                            cAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGAAGGCCATAGACTCCg                                                                                          >  2:3619164/1‑150 (MQ=255)
                                                              ggCGGTGTGCGGGTGATGCTATTCGGTATCTGGACCCGGCCGGGCCCACCCCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTGTCCTCCAGACGGCCATAGACTCCGGt                                                                                        <  2:605734/150‑1 (MQ=255)
                                                               gCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTg                                                                                       <  1:4226728/150‑1 (MQ=255)
                                                               gCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTg                                                                                       <  2:6838318/150‑1 (MQ=255)
                                                                  gtgtGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCAc                                                                                    <  2:2652023/150‑1 (MQ=255)
                                                                  gtgtGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAAAGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCAc                                                                                    >  2:177930/1‑150 (MQ=255)
                                                                    gtgCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCActg                                                                                  <  2:6446690/150‑1 (MQ=255)
                                                                     tgCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCActgc                                                                                 <  2:3783104/150‑1 (MQ=255)
                                                                        tGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAAAGCCCGTGTTACAAATAAACACACTTTCCGACAGACGGCCCTAGACTCCGGTGCACtgctgc                                                                              >  1:5854108/1‑150 (MQ=255)
                                                                           tGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCGGGAAAGCAGCCCCAAAGCCCGTGTTACAAAACAACAAACTTTACTACAGACGGCCCAATACGCCGGTGCCCAGCTGctct                                                                           >  1:2285321/1‑150 (MQ=255)
                                                                               gCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTcc                                                                       >  2:3370337/1‑150 (MQ=255)
                                                                                 tATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCgg                                                                     <  2:1143646/150‑1 (MQ=255)
                                                                                 tATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCgg                                                                     <  1:5207526/150‑1 (MQ=255)
                                                                                 tATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCgg                                                                     <  2:5083250/150‑1 (MQ=255)
                                                                                 tATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCGTGATACAAAAAAAAACACTTTCCTACAGACGGCCAAAGACACCGGTGCACTGCTGCTCGTTACgg                                                                     >  1:1400833/1‑150 (MQ=255)
                                                                                 tATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAAAGCCCTTGATAAAAATCAACACAGGTTCCTACAGACGGCCAAAGACTCCGGTGCCCTGCTGCTCTTTCCgg                                                                     >  1:358270/1‑150 (MQ=255)
                                                                                 tATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCGGCAAACCACCCCCAATGCCCGTGGTACCAAACAACACAGTTTCCAACCGACGGCCAGAGACTCCGGTGCACTGCGGCTCTTTCCgg                                                                     >  2:321532/1‑150 (MQ=255)
                                                                                  aTTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCggg                                                                    <  2:1031602/150‑1 (MQ=255)
                                                                                  aTTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCggg                                                                    <  2:4906835/150‑1 (MQ=255)
                                                                                      ggTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCGGGAAAGCAGCCCCCAAGCCCATGATAAAAAAAAAAACAGGTTACAACAGACCGCCCAAGGCACCGGGGCCCAGCTGCTCGTTCCgggggcg                                                                >  1:8579/1‑147 (MQ=255)
                                                                                          tCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTAccc                                                            >  2:3399559/1‑150 (MQ=255)
                                                                                           cTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCt                                                           <  2:4036944/150‑1 (MQ=255)
                                                                                            tGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAAAACAGTTTCCTACAGACGGCCAAAGACTCCGGTGCCCTGCTGCTCTTTACGGGGACTAcccgc                                                          >  2:3064423/1‑148 (MQ=255)
                                                                                             ggACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCACAGCAGCCCCAAAGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCa                                                         >  2:5573755/1‑150 (MQ=255)
                                                                                             ggACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAAAGCCCATGATACAAATCAAAACACTTTACTAAAGACGGCCAAAGACACCCGGGCCCAGCGGCTCTTTACGGGGACTACCcaca                                                         >  1:2015984/1‑150 (MQ=255)
                                                                                                 ccTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTACTACACACGGCCATAGACTCCGGTGCACTGCTGCTCTTTACGGGGACAACCCACATGcg                                                     >  2:4330684/1‑149 (MQ=255)
                                                                                                 ccTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATAAACACAGTTTCCTACAGACGGCCATAGACGCCGGTGCACGGCTGCTCTTTCCGGGGACAACCCTCATGCt                                                     >  2:871518/1‑150 (MQ=255)
                                                                                                 ccTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAAACAACACAGTTTCCTACAGACGGCCCTAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACGACCCACAAGca                                                     >  2:6496926/1‑149 (MQ=255)
                                                                                                 ccTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAAAGCCCTTGATACAAATAAACACACTTTCCTACAGACGGCCCTAGACTCCGGTGCACAGCGGCTCTTTCCGGGGCCTACCCTCAAGCt                                                     >  2:3686190/1‑150 (MQ=255)
                                                                                                  cTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTc                                                    <  2:1582694/150‑1 (MQ=255)
                                                                                                  cTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTc                                                    <  2:3836499/150‑1 (MQ=255)
                                                                                                     ccGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCCAAGCCCGTGGTACAAATCAAAACACTTTCCTACAGACGGCCCTAGACTCCGGTGCACTGCTGCTCTTTACGGGGCCTACCCGCAGGCTcgac                                                 >  2:764952/1‑150 (MQ=255)
                                                                                                     ccGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTcgac                                                 <  1:1658462/150‑1 (MQ=255)
                                                                                                     ccGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTcgac                                                 <  2:4702230/150‑1 (MQ=255)
                                                                                                     ccGGGCCCACACCAGCACTGGCTGATGGGCAATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTcgac                                                 <  2:5261312/150‑1 (MQ=255)
                                                                                                       gggCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTcgacga                                               <  2:6090844/150‑1 (MQ=255)
                                                                                                           ccACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTACTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGCCTACCCTCATGCTCGACGACCTc                                           >  2:6367538/1‑150 (MQ=255)
                                                                                                             acacCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAAACAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCACAAGCGCGACGACCTCCg                                         >  2:928114/1‑150 (MQ=255)
                                                                                                             acacCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCGTGATACAAATCAAAACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCACATGCGCGACGAccgctg                                         >  2:1334972/1‑146 (MQ=255)
                                                                                                               cccAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGtt                                       <  2:4634840/149‑1 (MQ=255)
                                                                                                                ccAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTg                                      <  2:4184310/150‑1 (MQ=255)
                                                                                                                 cAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTgg                                     <  2:2522569/150‑1 (MQ=255)
                                                                                                                 cAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAAACAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGAACTCCGGTgg                                     >  1:6776898/1‑150 (MQ=255)
                                                                                                                   gCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGa                                   <  1:4664177/150‑1 (MQ=255)
                                                                                                                       tGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATAcc                               >  2:3916005/1‑150 (MQ=255)
                                                                                                                         gCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGa                             <  1:6146397/150‑1 (MQ=255)
                                                                                                                         gCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGa                             >  2:5619964/1‑150 (MQ=255)
                                                                                                                              tGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGa                        <  1:1849179/150‑1 (MQ=255)
                                                                                                                                    aTTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATccgcc                  <  2:3525091/150‑1 (MQ=255)
                                                                                                                                    aTTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATccgcc                  <  2:1643710/150‑1 (MQ=255)
                                                                                                                                        attggggCTGCTTTGCAGCCCCCATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTa              >  1:3775411/4‑150 (MQ=255)
                                                                                                                                        attggggCTGCTTTGCAGCCCCCATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTa              >  1:3246689/4‑150 (MQ=255)
                                                                                                                                        attggggCTGCTTTGCAGCCCCCATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTa              >  1:4934839/4‑150 (MQ=255)
                                                                                                                                        attggggCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTa              >  1:5503618/4‑150 (MQ=255)
                                                                                                                                          gggggCTGCAAAGCAGCCCCAAGGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTACg            <  1:4842449/150‑1 (MQ=255)
                                                                                                                                           ggggCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTACGa           <  2:2043365/150‑1 (MQ=255)
                                                                                                                                              gCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTACGAtct        <  1:2904652/150‑1 (MQ=255)
                                                                                                                                                tGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTACGATCTCg      >  1:4401945/1‑150 (MQ=255)
                                                                                                                                                tGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTACGATCTCg      >  1:195873/1‑150 (MQ=255)
                                                                                                                                                 gCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTACGATCTCGc     >  1:3380034/1‑150 (MQ=255)
                                                                                                                                                    aaGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTACGATCTCGCCgg  >  1:5439973/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAACCGGGCTGGCTGCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCCACACCAGCACTGGCTGATGGGCCATTGACGGGGGCTGCAAAGCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACCGACCTGATCCGCCGCTACGATCTCGCCGG  >  NC_002947/4853453‑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: