Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,842,9410C.51.8% 87.0 / 129.5 56intergenic (+40/‑257)ccoP‑II/PP_4259cbb3‑type cytochrome c oxidase subunit/iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (14/13);  new base . (13/16);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 7.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCA  >  NC_002947/4842796‑4843089
                                                                                                                                                 |                                                                                                                                                     
gCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGAt                                                                                                                                                  <  1:336645/150‑1 (MQ=255)
  aGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGa                                                                                                                                                   <  2:1013801/148‑1 (MQ=255)
  aGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGa                                                                                                                                                   >  1:1013801/1‑148 (MQ=255)
  aGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATcc                                                                                                                                                <  1:174026/150‑1 (MQ=255)
  aGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATcc                                                                                                                                                >  2:467791/1‑150 (MQ=255)
    aTGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCt                                                                                                                                               >  2:293105/1‑150 (MQ=255)
        cGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCTTACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCtgtgta                                                                                                                                          <  1:414994/150‑6 (MQ=38)
         tgcacAGGAACAACTGCAGGGCAACGCCAAGTTGCACTTTCTGGCGTCTTAAGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAg                                                                                                                                         <  2:744645/146‑1 (MQ=39)
         ggCCCAGGAACAACTGCAGGTCAACGACAAGGTGCACGTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGCCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAg                                                                                                                                         <  1:274264/150‑1 (MQ=25)
           cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGa                                                                                                                                       <  1:293105/150‑1 (MQ=255)
           cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGa                                                                                                                                       <  1:1876235/150‑1 (MQ=255)
            ccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAg                                                                                                                                      <  1:844486/150‑1 (MQ=255)
            ccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAg                                                                                                                                      >  2:844486/1‑150 (MQ=255)
            ccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAg                                                                                                                                      <  1:467791/150‑1 (MQ=38)
               ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGGAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                                   <  2:1477493/150‑1 (MQ=255)
               ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                                   <  2:1711175/150‑1 (MQ=255)
               ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                                   <  1:284969/150‑1 (MQ=255)
               ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                                                                                                   <  2:674853/150‑1 (MQ=255)
               ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                                                                                                   <  2:1060039/150‑1 (MQ=255)
               ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGAAACCTCTGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                                                                                                   <  1:1573569/150‑1 (MQ=37)
                 aacaacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGc                                                                                                                                     >  1:407039/1‑146 (MQ=255)
                 aacaacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGc                                                                                                                                     <  2:407039/146‑1 (MQ=255)
                               aaCGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCggg                                                                                                                                  <  2:553372/135‑1 (MQ=255)
                               aaCGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCggg                                                                                                                                  >  1:553372/1‑135 (MQ=255)
                                                      ggcTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAg                                                                                                                                         >  1:205436/1‑105 (MQ=37)
                                                      ggcTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAg                                                                                                                                         <  2:205436/105‑1 (MQ=37)
                                                            cGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                                   >  2:622177/1‑105 (MQ=255)
                                                            cGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                                   <  1:622177/105‑1 (MQ=255)
                                                             gTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                                                                                                   <  1:1649953/104‑1 (MQ=25)
                                                             gTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                                                                                                   >  2:1649953/1‑104 (MQ=25)
                                                              tCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCtgtgt                                                                                                                                           >  1:1494993/1‑91 (MQ=255)
                                                              tCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCtgtgt                                                                                                                                           <  2:1494993/95‑5 (MQ=255)
                                                                         tCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATcc                                                                                                                                                >  1:664503/1‑80 (MQ=255)
                                                                         tCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATcc                                                                                                                                                <  2:664503/80‑1 (MQ=255)
                                                                                  gCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                                                                                                   >  2:1196766/1‑83 (MQ=18)
                                                                                  gCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTAGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                                                                                                   <  1:1196766/83‑1 (MQ=18)
                                                                                           gCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                                                                                                   <  1:1070347/74‑1 (MQ=17)
                                                                                           gCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                                                                                                   >  2:1070347/1‑74 (MQ=17)
                                                                                             cAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAg                                                                                                                                      >  1:395215/1‑69 (MQ=255)
                                                                                             cAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAg                                                                                                                                      <  2:395215/69‑1 (MQ=255)
                                                                                                                       tctTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGt                                                                                 <  1:166178/96‑1 (MQ=255)
                                                                                                                       tctTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGt                                                                                 >  2:166178/1‑96 (MQ=255)
                                                                                                                       tctTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGCTCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAAcgc                           >  1:903867/1‑150 (MQ=255)
                                                                                                                            gcgGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGCaga                                                                                                                    <  1:369965/56‑1 (MQ=255)
                                                                                                                            gcgGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAaga                                                                                                                    >  2:369965/1‑56 (MQ=255)
                                                                                                                                    aCCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGc              >  2:1055905/1‑150 (MQ=255)
                                                                                                                                     cccGCTCCTACACAGGATCCTTGTAGGAACGGCTCTAGCCGCGAAGAGGCCAGTACTGGCACAACATCTCTCGCAGGCAAGTGACCTGGATCAATTGACACCAGCCATAACACGTTTCATACCACGAACCCAACGCGACGAAAGGTGGCa             >  1:733882/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg            >  2:1580884/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg            >  1:1245659/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCTACAAGGGTCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg            >  1:1901359/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACACGACTAAAGGTCGCAg            >  2:378389/1‑150 (MQ=255)
                                                                                                                                        gCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTaaa                    <  2:1290905/140‑1 (MQ=255)
                                                                                                                                        gCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTaaa                    >  1:1290905/1‑140 (MQ=255)
                                                                                                                                        gCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTg          >  2:467052/1‑150 (MQ=255)
                                                                                                                                        gCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTg          >  1:1664959/1‑150 (MQ=255)
                                                                                                                                          tCCTACAAGGATCCTGTGTAGGAGCGGGATTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCg        >  1:1306172/1‑150 (MQ=255)
                                                                                                                                               cacaGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCa                                                                                    >  2:340618/1‑69 (MQ=255)
                                                                                                                                               cacaGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCa                                                                                    <  1:340618/69‑1 (MQ=255)
                                                                                                                                                acaGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCa                                                                                                 <  2:1143550/55‑1 (MQ=255)
                                                                                                                                                acaGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCa                                                                                                 >  1:1143550/1‑55 (MQ=255)
                                                                                                                                                acaGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCa  >  2:638308/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCA  >  NC_002947/4842796‑4843089

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