Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,079,1740TA38.2% 33.1 / 39.8 34intergenic (‑371/+26)mgtE/PP_t56magnesium transporter/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base T (10/11);  new base A (6/7);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCTATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGC  >  NC_002947/5079027‑5079322
                                                                                                                                                   |                                                                                                                                                    
gCTATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAAt                                                                                                                                                    <  2:237953/150‑1 (MQ=255)
          cgacgaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCt                                                                                                                                          >  2:1084347/1‑150 (MQ=255)
             cgaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCGCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCtttt                                                                                                                                       <  1:1084347/150‑1 (MQ=255)
                  gttgtGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATg                                                                                                                                  >  2:1063099/1‑150 (MQ=255)
                             gAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCgg                                                                                                                       >  1:559444/1‑150 (MQ=255)
                                        tGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGa                                                                                                            >  2:504685/1‑150 (MQ=255)
                                               tGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAAc                                                                                                     <  1:1063099/150‑1 (MQ=255)
                                                 aGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCttttt                                                                                                                                      >  1:1225516/1‑115 (MQ=255)
                                                 aGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCttttt                                                                                                                                      <  2:1225516/115‑1 (MQ=255)
                                                     cGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGa                                                                                               <  1:504685/150‑1 (MQ=255)
                                                      gCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATggg                                                                                                                                >  1:533272/1‑116 (MQ=255)
                                                      gCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATggg                                                                                                                                <  2:533272/116‑1 (MQ=255)
                                                          aGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCt                                                                                          >  1:1115028/1‑150 (MQ=255)
                                                            cGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCg                                                                                        <  2:1171441/150‑1 (MQ=255)
                                                                                            gagCGTAGATGTACATGCTCAAGAATTCTTTGGCGTAGAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGGAGGATCACAACCTGCGACCGCTCGGTGCGGAGCCGAGTACTCTATCCAGCTGAGCt                                                        <  1:167973/150‑1 (MQ=25)
                                                                                                gTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGAGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGa                                                    <  1:554902/150‑1 (MQ=255)
                                                                                                       gACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTg                                             >  1:989944/1‑150 (MQ=255)
                                                                                                          aGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCtttt                                                                                                                                       <  1:556114/57‑1 (MQ=255)
                                                                                                          aGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCtttt                                                                                                                                       >  2:556114/1‑57 (MQ=255)
                                                                                                             cTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATgggg                                                                                                                               >  2:389682/1‑62 (MQ=255)
                                                                                                             cTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATgggg                                                                                                                               <  1:389682/62‑1 (MQ=255)
                                                                                                                    aTTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAg                                <  2:559444/150‑1 (MQ=255)
                                                                                                                          tGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAgg                          <  2:93535/150‑1 (MQ=255)
                                                                                                                              gcACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTAc                      >  1:887316/1‑150 (MQ=255)
                                                                                                                                cccatcaaaAAGCCCGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCa                                                                <  2:635107/101‑1 (MQ=21)
                                                                                                                                cccatcaaaAAGCCCGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCa                                                                >  1:635107/6‑106 (MQ=21)
                                                                                                                                  catcaaaAAGCCCGCTCAATAGAGCGGGCTTTTTGTTGGGGTATGGCGGACTCAGCAGGATTCGAACCTGc                                                                                                 <  1:831701/68‑1 (MQ=18)
                                                                                                                                  catcaaaAAGCCCGCTCAATAGAGCGGGCTTTTTGTTGGGGTATGGCGGACTCAGCAGGATTCGAACCTGc                                                                                                 >  2:831701/4‑71 (MQ=18)
                                                                                                                                   atcaaaAAGCCCGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGa                                            >  2:1019289/3‑123 (MQ=255)
                                                                                                                                   atcaaaAAGCCCGCTCAATAGAGCGGGCTCATTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGa                                            <  1:1019289/121‑1 (MQ=18)
                                                                                                                                        aaaGCCCGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCg                                                                                        >  1:689865/1‑74 (MQ=25)
                                                                                                                                        aaaGCCCGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCg                                                                                        <  2:689865/74‑1 (MQ=25)
                                                                                                                                            cccGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCt                                                        <  1:1128698/102‑1 (MQ=255)
                                                                                                                                            cccGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCt                                                        >  2:1128698/1‑102 (MQ=255)
                                                                                                                                                  ctATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGc  <  2:989944/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCTATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGC  >  NC_002947/5079027‑5079322

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: ‑