Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,684,7490AC29.7% 159.6 / 14.0 74D164A (GAC→GCC) tusDNA replication terminus site‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (36/16);  new base C (0/22);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 7.14e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAT  >  NC_000913/1684701‑1684798
                                                |                                                 
cGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGAc                                                  >  1:7591182/1‑50 (MQ=255)
cGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGAc                                                  <  1:6983387/50‑1 (MQ=255)
 ggggCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGAcc                                                 >  1:4548359/1‑50 (MQ=255)
  gggCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGAccc                                                <  1:6261889/50‑1 (MQ=255)
   ggCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGAcccc                                               <  1:7135345/50‑1 (MQ=255)
   ggCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGAcccc                                               >  1:6573951/1‑50 (MQ=255)
    gCTGATCACCCTTAATGCTTACCGCGCGCTCACCGTTCTGCACGACCCCg                                              <  1:5078397/50‑1 (MQ=255)
    gCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCg                                              <  1:1128982/50‑1 (MQ=255)
     cTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGc                                             >  1:713649/1‑50 (MQ=255)
     cTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGc                                             <  1:1363730/50‑1 (MQ=255)
     cTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGc                                             <  1:6840504/50‑1 (MQ=255)
      tGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGc                                             >  1:4249593/1‑49 (MQ=255)
       gATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCa                                           >  1:2597906/1‑50 (MQ=255)
        aTCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCa                                           <  1:2378654/49‑1 (MQ=255)
         tCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACt                                         >  1:2639041/1‑50 (MQ=255)
            cccTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTa                                      >  1:3054415/1‑50 (MQ=255)
             ccTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTAc                                     >  1:7745096/1‑50 (MQ=255)
             ccTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTAc                                     >  1:3074629/1‑50 (MQ=255)
              cTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACg                                    >  1:2825220/1‑50 (MQ=255)
               ttAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGc                                   >  1:4348606/1‑50 (MQ=255)
                tAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCt                                  >  1:4932684/1‑50 (MQ=255)
                  aTGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCttt                                >  1:7740832/1‑50 (MQ=255)
                   tGCTTCCCGCACGCTCCCCTTTCTGCACGCCCCCGCCACTTTACGCTttg                               <  1:2695974/50‑1 (MQ=255)
                    tcTTCCCGCCCCCTCCCCGTTCTGCACGACCCCGCCACTTTACGCTttg                               <  1:7653789/48‑1 (MQ=37)
                      ttACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTttggtt                            >  1:6391845/1‑50 (MQ=255)
                        cccGCACGCTCACCTTTCTCCACGACCCCGCCACTTTACGCTttggtt                            <  1:2802975/47‑1 (MQ=255)
                        aCCGCGCGCTCCCGGTCCTCCACGCCCCCGCCACTTTACGCTttggttgg                          <  1:1085422/50‑1 (MQ=255)
                        aCCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTttggttgg                          >  1:2783611/1‑50 (MQ=255)
                         ccgcccgcTCCCTTTTCGGCACGACCCCGCCACTTTACGCTTTGGTTggg                         <  1:3397074/50‑1 (MQ=255)
                         ccgcccgcTCCCCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTggg                         <  1:3309326/50‑1 (MQ=255)
                         ccGCACGCTCCCCGTTCGGCACGCCCCCGCCACTTTACGCTttggttgg                          <  1:6733991/49‑1 (MQ=255)
                         ccGCACGCTCACCGTTCTGCACGCCCCCGCCACTTTACGCTttggttgg                          <  1:3292602/49‑1 (MQ=255)
                         ccGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTggg                         >  1:4432348/1‑50 (MQ=255)
                          cGCCCGCTCCCCGTTCTGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGc                        <  1:2724197/50‑1 (MQ=255)
                          cGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGc                        >  1:7690690/1‑50 (MQ=255)
                          cGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGc                        >  1:5457674/1‑50 (MQ=255)
                           gCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGc                        >  1:4934384/1‑49 (MQ=255)
                            cccGCTAACCTTTCGGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCta                      <  1:4428160/48‑1 (MQ=255)
                             aCGCTCCCCTTTCTGCCCGACCCCGCCACTTTACGCTTTGGTTGGGCtaa                     <  1:8432665/50‑1 (MQ=255)
                             aCGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCtaa                     >  1:6308283/1‑50 (MQ=255)
                             aCGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCtaa                     >  1:3344008/1‑50 (MQ=255)
                              cccTCCCCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCtaat                    <  1:1887537/48‑1 (MQ=255)
                               gCTCCCCTTTCCGCACGACCCCGCCACTTTACGCTTTGGTTGGGCtaata                   <  1:6472752/50‑1 (MQ=255)
                               gCTCACCTTTCTCCACGACCCCGCCACTTTACGCTTTGGTTGGGCtaata                   <  1:7116921/50‑1 (MQ=255)
                               gCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCtaat                    >  1:2836787/1‑49 (MQ=255)
                                cTCCCCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCtaataa                  <  1:1441216/50‑1 (MQ=255)
                                cTCACCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCtaataa                  <  1:2017840/50‑1 (MQ=255)
                                  cACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATaaa                 >  1:8080189/1‑49 (MQ=255)
                                  cACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATaaa                 >  1:2199989/1‑49 (MQ=255)
                                  cACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATaaa                 >  1:855131/1‑49 (MQ=255)
                                   aCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATaaa                 <  1:4889373/48‑1 (MQ=255)
                                   aCCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACa               >  1:902030/1‑50 (MQ=255)
                                   aCCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACa               >  1:2880350/1‑50 (MQ=255)
                                   aCCGTTCGACACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACa               <  1:3350725/50‑1 (MQ=255)
                                     cGTTCTGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACata             <  1:1241842/50‑1 (MQ=255)
                                     cGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACata             >  1:4452679/1‑50 (MQ=255)
                                      gTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACatat            <  1:388185/50‑1 (MQ=255)
                                      gTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACatat            >  1:5180264/1‑50 (MQ=255)
                                        tctcCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAGCATATc           <  1:452092/49‑1 (MQ=255)
                                        tCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACatat            <  1:7031521/48‑1 (MQ=255)
                                         cggCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCa          <  1:3217743/47‑1 (MQ=255)
                                            cACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCAt         <  1:5039994/47‑1 (MQ=255)
                                            cACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCAt         <  1:623038/47‑1 (MQ=255)
                                            cACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCAt         <  1:3376314/47‑1 (MQ=255)
                                            cACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTa       <  1:209670/49‑1 (MQ=255)
                                            cACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTa       <  1:4343589/49‑1 (MQ=255)
                                            cACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTa       <  1:855041/49‑1 (MQ=255)
                                            cACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCAtt        >  1:4998358/1‑48 (MQ=255)
                                            cACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCAtt        >  1:2478315/1‑48 (MQ=255)
                                            cACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCAtt        >  1:651106/1‑48 (MQ=255)
                                            cACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCAtt        >  1:7723820/1‑48 (MQ=255)
                                             aCGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAg     <  1:5411438/50‑1 (MQ=255)
                                               gACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGaa   >  1:2553116/1‑50 (MQ=255)
                                               gACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGaa   >  1:6732412/1‑50 (MQ=255)
                                                aCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAt  >  1:2693483/1‑50 (MQ=255)
                                                |                                                 
CGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAT  >  NC_000913/1684701‑1684798

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