Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi639,7240TG23.1% 66.4 / 26.9 39intergenic (+216/‑29)ahpC/ahpFalkyl hydroperoxide reductase, C22 subunit/alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base T (15/15);  new base G (5/4);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.88e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGAAATCAAGGAACTGTTG  >  NC_000913_3_bme_pgi/639581‑639872
                                                                                                                                               |                                                                                                                                                    
cccGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCatgatgat                                                                                                                                                <  2:392210/150‑1 (MQ=255)
   ggggATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGtt                                                                                                                                             >  4:125713/1‑150 (MQ=255)
     ggATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTa                                                                                                                                           <  3:152811/150‑1 (MQ=255)
           acaGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGccc                                                                                                                                     <  4:254323/150‑1 (MQ=255)
            caGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCa                                                                                                                                    <  1:111877/150‑1 (MQ=255)
                    ccgccTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATaa                                                                                                                            <  1:318434/150‑1 (MQ=255)
                        cTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACAt                                                                                                                        >  1:341435/1‑150 (MQ=255)
                         ttGCTCTGCCGCCAAATACTTCGGAAATTCACCTAATTCTTCGTGTGCTGCGGCACCCGCTTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATg                                                                                                                       <  3:162728/150‑1 (MQ=255)
                            ctctGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTc                                                                                                                    <  2:341435/150‑1 (MQ=255)
                                    cgAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGga                                                                                                                                 >  4:260084/1‑129 (MQ=255)
                                    cgAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGga                                                                                                                                 <  3:260084/129‑1 (MQ=255)
                                                   aTTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACCAATATGAAaactcaact                                                                                             <  2:99226/150‑1 (MQ=255)
                                                                 tcGGGTGCCGCAGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTg                                                                               >  1:462843/1‑150 (MQ=255)
                                                                   gggTGCCGCAGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTgag                                                                             >  1:473798/1‑150 (MQ=255)
                                                                    ggTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTgaga                                                                            >  3:26542/1‑150 (MQ=255)
                                                                        ctgcGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAAtt                                                                        <  3:125713/150‑1 (MQ=255)
                                                                                 ccGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAgg                                                                                        >  3:170194/1‑125 (MQ=255)
                                                                                 ccGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAgg                                                                                        <  4:170194/125‑1 (MQ=255)
                                                                                 ccGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGcc                                                               >  3:130091/1‑150 (MQ=255)
                                                                                 ccGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGcc                                                               >  2:155148/1‑150 (MQ=255)
                                                                                  cGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCt                                                              <  2:473798/150‑1 (MQ=255)
                                                                                    aTTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCt                                                                                                                     <  1:98871/93‑1 (MQ=255)
                                                                                    aTTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCt                                                                                                                     >  2:98871/1‑93 (MQ=255)
                                                                                                     atgatgCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGc                                           >  1:279946/1‑150 (MQ=255)
                                                                                                        atgCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTgg                                        >  2:465238/1‑150 (MQ=255)
                                                                                                            aaGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGa                                    <  2:462843/150‑1 (MQ=255)
                                                                                                              gCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACa                                  <  2:279946/150‑1 (MQ=255)
                                                                                                              gCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACa                                  >  3:267957/1‑150 (MQ=255)
                                                                                                                  cATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAgcgc                              >  1:177722/1‑150 (MQ=255)
                                                                                                                   aTCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAaactcaa                                                                                               >  2:299341/1‑84 (MQ=255)
                                                                                                                   aTCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAaactcaa                                                                                               <  1:299341/84‑1 (MQ=255)
                                                                                                                             ccGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAg                   >  2:269290/1‑150 (MQ=255)
                                                                                                                             ccGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAg                   >  2:269124/1‑150 (MQ=255)
                                                                                                                               acAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCt                                                                                 <  3:204557/85‑1 (MQ=255)
                                                                                                                               acAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCt                                                                                 >  4:204557/2‑86 (MQ=255)
                                                                                                                                      ccGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGAAATCAAgg          >  4:257743/1‑150 (MQ=255)
                                                                                                                                        gCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTg                                                                       >  1:101478/1‑87 (MQ=255)
                                                                                                                                        gCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTg                                                                       <  2:101478/87‑1 (MQ=255)
                                                                                                                                            gCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCCACACAAATATGAAAACTCAACTCAAGGCTTACCTTGCGAAATTGACCAAGCCTGTTGAGTTAATTGCCCCGCTGGATGACAGCGCTAAATCGGCAGAAATCAAGGAACTGt    <  1:4877/150‑1 (MQ=255)
                                                                                                                                             cATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGAAATCAAGGAACTGtt   <  4:130091/150‑1 (MQ=255)
                                                                                                                                              atgatgatgTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGAAATCAAGGAACTGTTg  <  4:267957/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGAAATCAAGGAACTGTTG  >  NC_000913_3_bme_pgi/639581‑639872

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