Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091312,0860CT47.4% 14.4 / 121.6 76intergenic (‑300/‑77)yaaI/dnaKDUF2541 domain‑containing protein YaaI/chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (19/21);  new base T (18/18);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGAT  >  NC_000913/12039‑12135
                                               |                                                 
ttAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAg                                                 >  1:703753/1‑50 (MQ=255)
ttAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAg                                                 >  1:116869/1‑50 (MQ=255)
 tagtagTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAg                                                 <  1:3781588/49‑1 (MQ=255)
 tagtagTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAg                                                 <  1:1225185/49‑1 (MQ=255)
  agtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGtt                                               >  1:415151/1‑50 (MQ=255)
   gtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGtt                                               >  1:535477/1‑49 (MQ=255)
   gtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGtt                                               <  1:3883150/49‑1 (MQ=255)
   gtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTg                                              <  1:4715461/50‑1 (MQ=255)
   gtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGtt                                               <  1:3033935/49‑1 (MQ=255)
    tagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGa                                             <  1:1855079/50‑1 (MQ=255)
    tagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTg                                              >  1:3177893/1‑49 (MQ=255)
    tagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGa                                             <  1:605730/50‑1 (MQ=255)
     agtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGa                                             >  1:5253962/1‑49 (MQ=255)
      gtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGaa                                            <  1:4480482/49‑1 (MQ=255)
      gtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGaaa                                           <  1:2573579/50‑1 (MQ=255)
       tCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAAc                                          <  1:5479755/50‑1 (MQ=255)
        cAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAAc                                          >  1:2931064/1‑49 (MQ=255)
         aaCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAAcc                                         <  1:3552215/49‑1 (MQ=255)
         aaCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAAc                                          >  1:3601256/1‑48 (MQ=255)
           ccGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGa                                      <  1:4869618/50‑1 (MQ=255)
            cGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGAc                                     <  1:267095/50‑1 (MQ=255)
             gCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACg                                    <  1:4271546/50‑1 (MQ=255)
             gCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACg                                    <  1:1204324/50‑1 (MQ=255)
             gCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACg                                    >  1:3481493/1‑50 (MQ=255)
             gCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACg                                    >  1:4509632/1‑50 (MQ=255)
               agtgagtgAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGtt                                  >  1:5251455/1‑50 (MQ=255)
               agtgagtgAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGtt                                  <  1:1591663/50‑1 (MQ=255)
                gtgagtgaGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGttt                                 <  1:3186715/50‑1 (MQ=255)
                gtgagtgaGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGttt                                 >  1:1190224/1‑50 (MQ=255)
                gtgagtgaGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGttt                                 >  1:3389210/1‑50 (MQ=255)
                 tgagtgagTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTc                                <  1:4141022/50‑1 (MQ=255)
                 tgagtgagTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTc                                <  1:4809074/50‑1 (MQ=255)
                  gagtgagtCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCg                               >  1:1120316/1‑50 (MQ=255)
                  gagtgagtCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCg                               >  1:1221289/1‑50 (MQ=255)
                  gagtgagtCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCg                               >  1:4661603/1‑50 (MQ=255)
                   agtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCg                               <  1:4143191/49‑1 (MQ=255)
                    gtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcc                             >  1:243934/1‑50 (MQ=255)
                     tgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGccc                            >  1:3512050/1‑50 (MQ=255)
                     tgagtCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGccc                            >  1:2316181/1‑50 (MQ=255)
                       agtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcccc                           <  1:1161916/49‑1 (MQ=255)
                        gtCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTa                         >  1:4738007/1‑50 (MQ=255)
                         tCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTAt                        >  1:5540471/1‑50 (MQ=255)
                          cTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTAt                        >  1:4714726/1‑49 (MQ=255)
                            gCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                      >  1:4182047/1‑49 (MQ=255)
                            gCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                      >  1:2173208/1‑49 (MQ=255)
                            gCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTa                      >  1:3362878/1‑49 (MQ=255)
                             cAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                      <  1:1086511/48‑1 (MQ=255)
                             cAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                      <  1:1418268/48‑1 (MQ=255)
                             cAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTAtt                       <  1:622755/47‑1 (MQ=255)
                             cAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTAtt                       <  1:2922035/47‑1 (MQ=255)
                             cAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTAtt                       <  1:610374/47‑1 (MQ=255)
                             cAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTa                      <  1:3604555/48‑1 (MQ=255)
                             cAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTa                      <  1:3955063/48‑1 (MQ=255)
                             cAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACa                    <  1:2924049/50‑1 (MQ=255)
                              aaaaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                      >  1:712523/1‑47 (MQ=255)
                              aaaaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                      >  1:3141619/1‑47 (MQ=255)
                                aaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGAc                 >  1:4287892/1‑50 (MQ=255)
                                aaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGAc                 <  1:618477/50‑1 (MQ=255)
                                aaaaaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGAc                 >  1:3980509/1‑50 (MQ=255)
                                 aaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt                >  1:5145073/1‑50 (MQ=255)
                                 aaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt                >  1:2799422/1‑50 (MQ=255)
                                 aaaaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt                <  1:626242/50‑1 (MQ=255)
                                 aaaaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt                <  1:4368633/50‑1 (MQ=255)
                                 aaaaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt                <  1:3742013/50‑1 (MQ=255)
                                    aTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTca              >  1:1159969/1‑49 (MQ=255)
                                     tGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTcaca            <  1:3066806/50‑1 (MQ=255)
                                      gAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACaa           >  1:5553113/1‑50 (MQ=255)
                                       aaaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAc          >  1:2172832/1‑50 (MQ=255)
                                       aaaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAc          <  1:4867923/50‑1 (MQ=255)
                                         aTTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACca        <  1:453423/50‑1 (MQ=255)
                                          ttGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACca        >  1:4196039/1‑49 (MQ=255)
                                           tGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACcaca      >  1:1563878/1‑50 (MQ=255)
                                            gggTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACcaca      <  1:4804723/49‑1 (MQ=255)
                                             ggCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACatg    <  1:1652028/50‑1 (MQ=255)
                                             ggCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACat     >  1:2462903/1‑49 (MQ=255)
                                              gCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACatga   <  1:180970/50‑1 (MQ=255)
                                               cAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACatgat  >  1:2387615/1‑50 (MQ=255)
                                               cAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACatgat  <  1:3021174/50‑1 (MQ=255)
                                               cAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACatgat  <  1:4386622/50‑1 (MQ=255)
                                               |                                                 
TTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGAT  >  NC_000913/12039‑12135

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