Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,875,0501.C35.0% 135.3 / 143.4 80coding (527/822 nt)nimRDNA‑binding transcriptional repressor NimR
Reads supporting (aligned to +/- strand):  ref base . (30/22);  new base C (14/14);  total (44/36)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 106 total aligned reads displayed.

GGCAAGATTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCACTCGACAGGCCCTTTCGCCATCATCTCCACCATCGTGCGG  >  NC_000913/1875001‑1875097
                                                  |                                               
ggCAAGATTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATg                                                  >  1:5139671/1‑50 (MQ=255)
cgCAAGATTGCGTTCACTCATCGCATAAAAGCCAGCCCATTGCCCCAATg                                                  >  2:4730854/2‑50 (MQ=255)
 gCAAGATTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGc                                                >  1:5364855/1‑50 (MQ=255)
 gCAAGATTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGc                                                >  2:294482/1‑50 (MQ=255)
 gCAAGATTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGc                                                >  1:1801578/1‑50 (MQ=255)
  cAAGATTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGcc                                               >  2:4119776/1‑50 (MQ=255)
  cAAGATTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGcc                                               <  1:623730/50‑1 (MQ=255)
  cAAGATTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGcc                                               <  2:3952755/50‑1 (MQ=255)
  cAAGATTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGcc                                               >  2:5049654/1‑50 (MQ=255)
   aaGATTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGccc                                              >  1:48866/1‑50 (MQ=255)
   aaGATTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGccc                                              >  2:2612781/1‑50 (MQ=255)
   aaGATTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGccc                                              <  1:4924069/50‑1 (MQ=255)
   aaGATTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGccc                                              >  1:5080011/1‑50 (MQ=255)
   aaGATTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGccc                                              >  2:4665159/1‑50 (MQ=255)
   aaGATTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGccc                                              >  1:1824668/1‑50 (MQ=255)
    aGATTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGcccc                                             >  1:2364275/1‑50 (MQ=255)
    aGATTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGcccc                                             >  1:2951734/1‑50 (MQ=255)
    aGATTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGcccc                                             >  1:3119939/1‑50 (MQ=255)
    aGATTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGcccc                                             >  2:211414/1‑50 (MQ=255)
    aGATTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGcccc                                             <  1:3893380/50‑1 (MQ=255)
      aTTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGcccccc                                           >  2:4523976/1‑50 (MQ=255)
      aTTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGccccc                                            >  2:247486/1‑49 (MQ=255)
       ttGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGccccccc                                          <  1:3687366/50‑2 (MQ=255)
       ttGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGcccccc                                           >  1:1747393/1‑49 (MQ=255)
          cGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCACTc                                       <  1:1153718/50‑1 (MQ=255)
          cGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCACTc                                       >  2:1355766/1‑50 (MQ=255)
          cGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCACTc                                       >  1:3140635/1‑50 (MQ=255)
            ttCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCACTCGa                                     >  1:4600041/1‑50 (MQ=255)
             tCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCACTCGa                                     >  1:3001768/1‑49 (MQ=255)
              cACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCACTCGa                                     >  2:4522780/1‑48 (MQ=255)
              cACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCACGCGAc                                    <  2:4133579/49‑1 (MQ=255)
               aCTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCCACTCGa                                     >  1:3821083/1‑48 (MQ=38)
               aCTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCACTCGACa                                   <  1:1384167/49‑1 (MQ=255)
               aCTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCACTCGACAg                                  <  2:612923/50‑1 (MQ=255)
                cTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCACTCGACAgg                                 >  1:1192222/1‑50 (MQ=255)
                 tCATCGCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCACTCGACAGGc                                >  1:4215849/1‑50 (MQ=255)
                   aTCGCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCCACTCGACAGGcc                               >  2:3159100/1‑50 (MQ=39)
                   aTCGCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCACTCGACAGGccc                              >  2:2580864/1‑50 (MQ=255)
                   aTCGCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCACTCGACAGGcc                               <  1:657611/49‑1 (MQ=255)
                    tCGCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCCACTCGACAGGccc                              <  2:1434993/50‑1 (MQ=39)
                    tCGCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCACTCGACAGGCCCt                             >  2:2719555/1‑50 (MQ=255)
                      gCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCACTCGACAGGCCCttt                           <  2:5318586/50‑1 (MQ=255)
                      gCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCACTCGACAGGCCCtt                            <  2:4426979/49‑1 (MQ=255)
                      gCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCACTCGACAGGCCCtt                            <  1:3047341/49‑1 (MQ=255)
                       cAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCCACTCGACAGGCCCttt                           >  2:680822/1‑50 (MQ=39)
                       cAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCACTCGACAGGCCCTTTc                          <  2:3440007/50‑1 (MQ=255)
                        aaaaaaGCCAGCCCATTGCCCCAATGCCCCCCCACTCGACAGGCCCTTTc                          <  1:1839304/50‑1 (MQ=39)
                        aaaaaaGCCAGCCCATTGCCCCAATGCCCCCCACTCGACAGGCCCTTTCg                         >  1:3498369/1‑50 (MQ=255)
                        aaaaaaGCCAGCCCATTGCCCCAATGCCCCCCACTCGACAGGCCCTTTCg                         <  1:3242084/50‑1 (MQ=255)
                         aaaaaGCCAGCCCATTGCCCCAATGCCCCCCACTCGACAGGCCCTTTCGc                        <  1:1798300/50‑1 (MQ=255)
                          aaaaGCCAGCCCATTGCCCCAATGCCCCCCCACTCGACATGCCCTTACGc                        >  2:3729280/1‑50 (MQ=37)
                          aaaaGCCAGCCCATTGCCCCAATGCCCCCCCACTCGACAGGCCCTTTCGc                        >  2:544327/1‑50 (MQ=39)
                          aaaaGCCAGCCCATTGCCCCAATGCCCCCCCACTCGACAGGCCCTTTCGc                        >  1:3858866/1‑50 (MQ=39)
                          aaaaGCCAGCCCATTGCCCCAATGCCCCCCACTCGACAGGCCCTTTCGc                        >  2:3883683/1‑49 (MQ=255)
                           aaaGCCAGCCCATTGCCCCAATGCCCCCCCACTCGACAGGCCCTTTCGc                        >  2:610400/1‑49 (MQ=38)
                            cagccagccCATTGCCCCAATGCCCCCCCACTCGACAGGCCCTTTCGCca                      >  2:3166215/2‑50 (MQ=38)
                            aaGCCAGCCCATTGCCCCAATGCCCCCCCACTCGACAGGCCCTTTCGCca                      >  2:501443/1‑50 (MQ=39)
                            aaGCCAGCCCATTGCCCCAATGCCCCCCACTCGACAGGCCCTTTCGCcat                     <  1:3591690/50‑1 (MQ=255)
                             agccagccCATTGCCCCAATGCCCCCCCACTCGACAGGCCCTTTCGCcat                     <  2:1376074/50‑1 (MQ=39)
                             agccagccCATTGCCCCAATGCCCCCCCACTCGACAGGCCCTTTCGCcat                     <  1:1677225/50‑1 (MQ=39)
                             agccagccCATTGCCCCAATGCCCCCCACTCGACAGGCCCTTTCGCcat                     >  1:3916820/1‑49 (MQ=255)
                              gccagccCATTGCCCCAATGCCCCCCCACTCGACAGGCCCTTTCGccct                     >  2:3910152/1‑47 (MQ=38)
                              gccagccCATTGCCCCAATGCCCCCCACTCGACAGGCCCTTTCGCcatca                   <  2:4896085/50‑1 (MQ=255)
                                cagccCATTGCCTCAATGCCCCCCACTCGACAGGCCCTTTCGCCatcatc                 >  2:2765465/1‑50 (MQ=255)
                                cagccCATTGCCCCAATGCCCCCCACTCGACAGGCCCTTTCGCCatcatc                 >  1:5223824/1‑50 (MQ=255)
                                 agccCATTGCCCCAATGCCCCCCACTCGACAGGCCCTTTCGCCATCAtct                <  2:1756189/50‑1 (MQ=255)
                                 agccCATTGCCCCAATGCCCCCCACTCGACAGGCCCTTTCGCCATCAtct                >  2:3889516/1‑50 (MQ=255)
                                  gccCATTGCCCCAATGCCCCCCCACTCGACAGGCCCTTTCGCCATCAtct                <  2:4751386/50‑1 (MQ=39)
                                  gccCATTGCCCCAATGCCCCCCCACTCGACAGGCCCTTTCGCCATCAtct                <  2:2901984/50‑1 (MQ=39)
                                   cccATTGCCCCAATGCCCCCCCACTCGACAGGCCCTTTCGCCATCAtctc               <  1:5280387/50‑1 (MQ=39)
                                   cccATTGCCCCAATGCCCCCCACTCGACAGGCCCTTTCGCCATCAtctc               >  1:1042631/1‑49 (MQ=255)
                                    ccATTGCCCCAATGCCCCCCACTCGACAGGCCCTTTCGCCATCATCTcca             >  1:3365702/1‑50 (MQ=255)
                                     cATTGCCCCAATGCCCCCCCACTCGACAGGCCCTTTCGCCATCATCTcca             >  2:4318290/1‑50 (MQ=255)
                                     cATTGCCCCAATGCCCCCCCACCCGACAGGCCCTTTCGCCATCAGCTcca             <  2:2185124/50‑1 (MQ=255)
                                     cATTGCCCCAATGCCCCCCACTCGACAGGCCCTTTCGCCATCATCTccac            >  1:429169/1‑50 (MQ=255)
                                     cATTGCCCCAATGCCCCCCACTCGACAGGCCCTTTCGCCATCATCTccac            <  1:4675070/50‑1 (MQ=255)
                                     cATTGCCCCAATGCCCCCCACTCGACAGGCCCTTTCGCCATCATCTccac            >  2:3975366/1‑50 (MQ=255)
                                     cATTGCCCCAATGCCCCCCACTCGAAAGCCCCTTTCGCCATCATCTccac            >  1:3382948/1‑50 (MQ=39)
                                      aTTGCCCCAATGCCCCCCACTCGACAGGCCCTTTCGCCATCATCTccac            >  1:1926186/1‑49 (MQ=255)
                                       ttGCCCCAATGCCCCCCACTCGACAGGCCCTTTCGCCATCATCTccacca          <  1:1837168/50‑1 (MQ=255)
                                         gCCCCAATGCCCCCCACTCGACAGGCCCTTTCGCCATCATCTCCACCAt         <  1:99302/49‑1 (MQ=255)
                                          ccccAATGCCCCCCCACTCGACAGGCCCTTTCGCCATCATCTCCACCAt         <  1:2638020/49‑1 (MQ=39)
                                          ccccAATGCCCCCCCACTCGACAGGCCCTTTCGCCATCATCTCCACCAt         <  1:3947320/49‑1 (MQ=39)
                                          ccccAATGCCCCCCCACTCGACAGGCCCTTTCGCCATCATCTCCACCATc        >  2:2472499/1‑50 (MQ=255)
                                          ccccAATGCCCCCCCACTCGACAGGCCCTTTCGCCATCATCTCCACCATc        <  1:3685583/50‑1 (MQ=255)
                                          ccccAATGCCCCCCCACTCGACAGGCCCTTTCGCCATCATCTCCACCATc        <  2:2877408/50‑1 (MQ=255)
                                          ccccAATGCCCCCCCACTCGACAGGCCCTTTCGCCATCATCTCCACCATc        >  2:5273884/1‑50 (MQ=255)
                                          ccccAATGCCCCCCCACTCGACAGGCCCTTGCGCCATCATCTCCACCATc        <  1:2608966/50‑1 (MQ=255)
                                          ccccAATGCCCCCCACTCGACAGGCCCTTTCGCCATCATCTCCACCATCg       >  2:2581132/1‑50 (MQ=255)
                                          ccccAATGCCCCCCACTCGACAGGCCCTTTCGCCATCATCTCCACCATCg       <  1:3353648/50‑1 (MQ=255)
                                          ccccAATGCCCCCCACTCGACAGGCCCTTTCGCCATCATCTCCACCATCg       >  1:573005/1‑50 (MQ=255)
                                          ccccAATGCCCCCCACTCGACAGGCCCTTTCGCCATCATCTCCACCATCg       >  2:2834350/1‑50 (MQ=255)
                                           cccAATGCCCCCCACTCGACAGGCCCTTTCGCCATCATCTCCACCATCGt      <  2:3374483/50‑1 (MQ=255)
                                            ccAATGCCCCCCACTCGACAGGCCCTTTCGCCATCATCTCCACCATCGTg     >  2:2804037/1‑50 (MQ=255)
                                             cAATGCCCCCCCACTCGACAGGCCCTTTCGCCATCATCTCCACCATCGTg     <  2:1944259/50‑1 (MQ=255)
                                               aTGCCCCCCACTCGACAGGCCCTTTCGCCATCATCTCCACCATCGTGCgg  <  1:4095300/50‑1 (MQ=255)
                                               aTGCCCCCCACTCGACAGGCCCTTTCGCCATCATCTCCACCATCGTGCgg  >  2:4684115/1‑50 (MQ=255)
                                               aTGCCCCCCACTCGACAGGCCCTTTCGCCATCATCTCCACCATCGTGCgg  <  1:2615031/50‑1 (MQ=255)
                                               aTGCCCCCCACTCGACAGGCCCTTTCGCCATCATCTCCACCATCGTGCg   >  1:4069282/1‑49 (MQ=255)
                                               aTGCCCCCCACTCGACAGGCCCTTTCGCCATCATCTCCACCATCGTGCg   >  2:4046680/1‑49 (MQ=255)
                                                  |                                               
GGCAAGATTGCGTTCACTCATCGCAAAAAAGCCAGCCCATTGCCCCAATGCCCCCCACTCGACAGGCCCTTTCGCCATCATCTCCACCATCGTGCGG  >  NC_000913/1875001‑1875097

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