Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091312,0860CT48.4% 23.7 / 197.2 126intergenic (‑300/‑77)yaaI/dnaKDUF2541 domain‑containing protein YaaI/chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (27/37);  new base T (24/36);  total (52/74)
Fisher's exact test for biased strand distribution p-value = 8.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 127 total aligned reads displayed.

TTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGAT  >  NC_000913/12038‑12135
                                                |                                                 
tttAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCa                                                  <  1:4426296/50‑1 (MQ=255)
  tagtagTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAg                                                 <  2:3229599/49‑1 (MQ=255)
  tagtagTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGt                                                >  1:2811039/1‑50 (MQ=255)
  tagtagTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGt                                                >  2:484938/1‑50 (MQ=255)
   agtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGt                                                <  1:1845230/49‑1 (MQ=255)
    gtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTg                                              <  2:2499012/50‑1 (MQ=255)
    gtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTg                                              <  2:1779156/50‑1 (MQ=255)
    gtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGtt                                               <  2:3218423/49‑1 (MQ=255)
    gtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGtt                                               <  1:4190701/49‑1 (MQ=255)
     tagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGa                                             <  1:1567888/50‑1 (MQ=255)
     tagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGa                                             <  1:1170536/50‑1 (MQ=255)
       gtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGGAGTTGaaa                                           <  2:1302270/50‑1 (MQ=255)
       gtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGaaa                                           <  2:2046884/50‑1 (MQ=255)
       gtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGaaa                                           <  1:4266677/50‑1 (MQ=255)
       gtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGaaa                                           <  1:1415478/50‑1 (MQ=255)
       gtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGaa                                            <  2:4192190/49‑1 (MQ=255)
          aaCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAAcc                                         <  1:2155161/49‑1 (MQ=255)
            ccGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGa                                      >  1:437830/1‑50 (MQ=255)
            ccGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGa                                      <  1:1800396/50‑1 (MQ=255)
             cGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGAc                                     <  2:448503/50‑1 (MQ=255)
             cGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGAc                                     <  2:805675/50‑1 (MQ=255)
              gCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACg                                    <  2:3678457/50‑1 (MQ=255)
              gCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACg                                    >  1:3736577/1‑50 (MQ=255)
              gCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACg                                    >  2:2177096/1‑50 (MQ=255)
              gCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACg                                    <  1:3751833/50‑1 (MQ=255)
              gCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACg                                    >  2:1332593/1‑50 (MQ=255)
               cAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGt                                   >  1:715046/1‑50 (MQ=255)
               cAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGt                                   >  2:1121072/1‑50 (MQ=255)
               cAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGt                                   >  2:2407294/1‑50 (MQ=255)
                agtgagtgAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGtt                                  >  2:3784469/1‑50 (MQ=255)
                agtgagtgAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGtt                                  >  1:2731161/1‑50 (MQ=255)
                agtgagtgAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGtt                                  >  1:3935891/1‑50 (MQ=255)
                agtgagtgAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGtt                                  <  1:1185283/50‑1 (MQ=255)
                 gtgagtgaGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGttt                                 >  1:4156880/1‑50 (MQ=255)
                 gtgagtgaGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGttt                                 >  2:745930/1‑50 (MQ=255)
                  tgagtgagTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTc                                <  2:4061101/50‑1 (MQ=255)
                  tgagtgagTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTc                                >  2:1486083/1‑50 (MQ=255)
                  tgagtgagTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTc                                <  1:1742388/50‑1 (MQ=255)
                  tgagtgagTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTc                                <  1:1561505/50‑1 (MQ=255)
                  tgagtgagTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTc                                <  1:1468574/50‑1 (MQ=255)
                   gagtgagtCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCg                               >  2:1333941/1‑50 (MQ=255)
                   gagtgagtCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCg                               >  2:1239662/1‑50 (MQ=255)
                    agtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGc                              <  2:4147542/50‑1 (MQ=255)
                    agtgagtCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGc                              >  1:2554867/1‑50 (MQ=255)
                      tgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGccc                            >  2:307047/1‑50 (MQ=255)
                      tgagtCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGccc                            <  1:1559907/50‑1 (MQ=255)
                        agtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcccc                           >  1:809196/1‑49 (MQ=255)
                        agtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCt                          <  1:305652/50‑1 (MQ=255)
                        agtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCt                          <  2:3771619/50‑1 (MQ=255)
                        agtCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGcccc                           >  1:925100/1‑49 (MQ=255)
                         gtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTa                         <  1:1512320/50‑1 (MQ=255)
                          tCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTAt                        >  2:2968106/1‑50 (MQ=255)
                            tGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTa                      >  1:2455074/1‑50 (MQ=255)
                              cAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTAtt                       <  2:44209/47‑1 (MQ=255)
                              cAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTAtt                       <  2:4116612/47‑1 (MQ=255)
                              cAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                      >  2:3230971/1‑48 (MQ=255)
                              cAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                      >  2:4514068/1‑48 (MQ=255)
                              cAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTAc                     <  2:2563034/49‑1 (MQ=255)
                              cAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTAtt                       <  1:1646907/47‑1 (MQ=255)
                              cAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTAtt                       <  1:1334591/47‑1 (MQ=255)
                              cAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTa                      >  2:791110/1‑48 (MQ=255)
                              cAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTAc                     <  1:2031247/49‑1 (MQ=255)
                               aaaaaaaTGAAATTGGGGAGTTGAAACCAGACGTTTCGCCCCTATTa                      >  2:707818/1‑47 (MQ=255)
                               aaaaaaaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTa                      >  2:14028/1‑47 (MQ=255)
                                 aaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGAc                 <  1:3315750/50‑1 (MQ=255)
                                  aaaaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt                <  2:1967049/50‑1 (MQ=255)
                                  aaaaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt                <  1:3787326/50‑1 (MQ=255)
                                  aaaaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt                <  1:548930/50‑1 (MQ=255)
                                  aaaaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt                <  2:1162943/50‑1 (MQ=255)
                                   aaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACta               <  2:1691090/50‑2 (MQ=255)
                                   aaaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt                >  2:3465275/1‑49 (MQ=255)
                                    aaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTca              >  1:668507/1‑50 (MQ=255)
                                     aggAAATTGGGTAGTTGAAACCAGCCGTTTCCCCCCTATTACAGACTcac             <  2:3450481/48‑1 (MQ=255)
                                     aTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTcac             <  1:3549470/50‑1 (MQ=255)
                                     aTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTca              >  2:1990250/1‑49 (MQ=255)
                                      tGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTcaca            <  1:1154740/50‑1 (MQ=255)
                                      tGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTcaca            <  1:112347/50‑1 (MQ=255)
                                      tGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTcaca            <  2:3081016/50‑1 (MQ=255)
                                      tGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTcaca            <  1:2330697/50‑1 (MQ=255)
                                      tGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTcaca            <  1:4139432/50‑1 (MQ=255)
                                      tGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTcaca            <  2:1242810/50‑1 (MQ=255)
                                        aaaTTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACaa           <  2:2972562/49‑1 (MQ=255)
                                        aaaTTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACaa           <  2:2528640/49‑1 (MQ=255)
                                        aaaTTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAc          <  1:1614740/50‑1 (MQ=255)
                                        aaaTTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAc          <  2:1644203/50‑1 (MQ=255)
                                         aaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAc          >  2:3951725/1‑49 (MQ=255)
                                          aTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACca        >  1:3738251/1‑50 (MQ=255)
                                          aTTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACca        <  2:1418133/50‑1 (MQ=255)
                                          aTTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACca        >  2:1021147/1‑50 (MQ=255)
                                           ttGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACcac       <  2:17902/50‑1 (MQ=255)
                                           ttGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACcac       <  2:2309242/50‑1 (MQ=255)
                                           ttGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACca        >  2:3698737/1‑49 (MQ=255)
                                           ttGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCGCAACcac       <  2:1348607/50‑1 (MQ=255)
                                             gggTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACat     >  1:992902/1‑50 (MQ=255)
                                              ggCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACat     <  2:2743985/49‑1 (MQ=255)
                                              ggCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACat     >  2:3424700/1‑49 (MQ=255)
                                              ggCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACat     <  2:3356276/49‑1 (MQ=255)
                                               gCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACatga   <  2:2629896/50‑1 (MQ=255)
                                               gCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACatga   <  2:813480/50‑1 (MQ=255)
                                                cAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACatgat  <  1:4012359/50‑1 (MQ=255)
                                                |                                                 
TTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGAT  >  NC_000913/12038‑12135

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