Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913310,9340AC29.3% 80.3 / 10.7 42F135V (TTC→GTC) ecpRDNA‑binding transcriptional dual regulator MatA
Reads supporting (aligned to +/- strand):  ref base A (20/9);  new base C (0/12);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 5.32e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCGTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCC  >  NC_000913/310867‑311001
                                                                   |                                                                   
gCCGTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGaa                                                                     >  1:3247112/1‑68 (MQ=255)
gCCGTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGaa                                                                     >  1:251150/1‑68 (MQ=255)
gCCGTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGaa                                                                     >  1:3232857/1‑68 (MQ=255)
gCCGTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGaa                                                                     >  1:747513/1‑68 (MQ=255)
 ccGTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTNAGGAAc                                                                    <  1:2036981/68‑1 (MQ=255)
 ccGTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAAc                                                                    >  1:395726/1‑68 (MQ=255)
 ccGTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAAc                                                                    >  1:635555/1‑68 (MQ=255)
          ggATAATTTGCATTTCCCGGTCCGTTTTTTTATTTTTCTTAATATCTTTTCTCAGGACCCTGCCATTa                                                           <  1:115134/68‑1 (MQ=255)
          ggATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTa                                                           >  1:2822276/1‑68 (MQ=255)
          ggATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTa                                                           >  1:2921662/1‑68 (MQ=255)
            aTAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAAt                                                         <  1:2968038/68‑1 (MQ=255)
               ttttccatttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCACAGAACCCTGCCATTAATcac                                                      <  1:2967182/67‑1 (MQ=255)
               atttccatttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAGCCTGCCATTAATcac                                                      <  1:213090/68‑1 (MQ=255)
               atttccatttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATcac                                                      <  1:1394867/68‑1 (MQ=255)
               atttccattcccCGGTCCGTGTTTTTATCTTTCTTAATATCTTTTCTGAGGAACCTGCCATTAATcac                                                      <  1:3308702/68‑1 (MQ=255)
               atttccattcccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATcac                                                      <  1:73988/68‑1 (MQ=255)
                  tccatttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATa                                                   >  1:2442634/1‑68 (MQ=255)
                    catttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGa                                                 <  1:1100179/68‑1 (MQ=255)
                    catttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGa                                                 <  1:2628506/68‑1 (MQ=255)
                    catttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGa                                                 <  1:3277982/68‑1 (MQ=255)
                    catttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGa                                                 <  1:2600805/68‑1 (MQ=255)
                    catttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGa                                                 <  1:1303434/68‑1 (MQ=255)
                    catttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCCCAGGAACCTGCCATTAATCACATAGa                                                 <  1:3255848/68‑1 (MQ=255)
                    cattcccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGa                                                 <  1:1131319/68‑1 (MQ=255)
                    cattcccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGa                                                 <  1:2732803/68‑1 (MQ=255)
                     atttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGAc                                                >  1:1568541/1‑68 (MQ=255)
                      tttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGAcc                                               >  1:544945/1‑68 (MQ=255)
                         cccGGTCCGTGGTTTTTTTTTTTTTAATTTTTTTTCCCAGGAACCTGCCATTAATCACATAGACCAGt                                            <  1:414444/68‑1 (MQ=255)
                         cccGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGt                                            <  1:839919/68‑1 (MQ=255)
                         cccGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGt                                            <  1:1540296/68‑1 (MQ=255)
                          ccGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGtt                                           >  1:3240178/1‑68 (MQ=255)
                                    aTTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTCTGTAat                                 >  1:637419/1‑68 (MQ=255)
                                     ttttATCTTTCTTAATAGCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAata                                <  1:1966711/68‑1 (MQ=255)
                                      tttATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAatat                               >  1:2726359/1‑68 (MQ=255)
                                       ttATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATc                              <  1:3451893/68‑1 (MQ=255)
                                         aTCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCAc                            >  1:372760/1‑68 (MQ=255)
                                         aTCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCAc                            >  1:2511698/1‑68 (MQ=255)
                                          tctttcttAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACg                           <  1:562958/68‑1 (MQ=255)
                                          tctttcttAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACg                           >  1:415289/1‑68 (MQ=255)
                                            ttttttAATATCTTTCCCAAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGAc                         <  1:3151944/68‑1 (MQ=255)
                                               tttAATTTCTTTCCCAAGAACCCTGCTATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCa                      <  1:2722629/67‑1 (MQ=255)
                                               cttAATATTTTTTCTCAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCa                      <  1:718122/68‑1 (MQ=255)
                                                 tAATTTCTTTTCTCAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt                    <  1:2795896/68‑1 (MQ=255)
                                                 tAATATTTTTTCTAAGAAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt                    <  1:1621973/68‑1 (MQ=255)
                                                 tAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt                    <  1:2272449/68‑1 (MQ=255)
                                                 tAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt                    <  1:1903088/68‑1 (MQ=255)
                                                 tAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt                    <  1:443749/68‑1 (MQ=255)
                                                 tAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt                    <  1:2195209/68‑1 (MQ=255)
                                                 tAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt                    <  1:2965193/68‑1 (MQ=255)
                                                 tAATATCTTTTCTCAGCAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt                    <  1:321438/68‑1 (MQ=255)
                                                  aaTATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTc                   >  1:1251122/1‑68 (MQ=255)
                                                     tttttttCACAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg                <  1:1833877/65‑1 (MQ=255)
                                                     ttCTTTTGTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg                <  1:247638/67‑1 (MQ=255)
                                                     atCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg                >  1:171583/1‑68 (MQ=255)
                                                     atCTTTTCTCAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg                <  1:1932396/68‑1 (MQ=255)
                                                      ttttttCTCGGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc               <  1:1433800/66‑1 (MQ=255)
                                                      ttttttCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTTAGTCCAGc               <  1:919826/66‑1 (MQ=255)
                                                      ttttttCTCAGGAACCTGCCATTAATCACATATACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc               <  1:2668298/66‑1 (MQ=255)
                                                      ttttttCCCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc               <  1:2067950/66‑1 (MQ=255)
                                                      tCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc               <  1:1564821/68‑1 (MQ=255)
                                                      tCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc               >  1:2805995/1‑68 (MQ=255)
                                                      tCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc               >  1:2642181/1‑68 (MQ=255)
                                                      tCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc               >  1:1787494/1‑68 (MQ=255)
                                                      tCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc               >  1:1357106/1‑68 (MQ=255)
                                                      tCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc               >  1:1215143/1‑68 (MQ=255)
                                                      tCTTTTCTCAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc               <  1:3008478/68‑1 (MQ=255)
                                                      tCTTTTCCCAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc               <  1:1539402/68‑1 (MQ=255)
                                                      tCTTTCCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc               <  1:1979935/68‑1 (MQ=255)
                                                      tCTTTCCTAAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc               <  1:1120002/68‑1 (MQ=255)
                                                        ttttCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGt             <  1:3248939/68‑1 (MQ=255)
                                                        ttttCCCAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGt             <  1:419009/68‑1 (MQ=255)
                                                                ggAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACg     >  1:2632218/1‑68 (MQ=255)
                                                                   aCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGccc  >  1:225860/1‑68 (MQ=255)
                                                                   |                                                                   
GCCGTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCC  >  NC_000913/310867‑311001

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