Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,944,1970TC30.8% 250.6 / 146.4 156*262Q (TAA→CAA) znuBZn(2(+)) ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (54/54);  new base C (31/17);  total (85/71)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 160 total aligned reads displayed.

TGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCC  >  NC_000913/1944149‑1944245
                                                |                                                
tgtgCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCt                                                  <  2:5321559/49‑1 (MQ=255)
tgtgCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCt                                                  >  2:2255777/1‑49 (MQ=255)
tgtgCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCt                                                  >  1:4415827/1‑49 (MQ=255)
  tgCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCCaa                                                <  1:2945643/49‑1 (MQ=255)
    cGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATc                                              <  1:3277854/49‑1 (MQ=255)
    cGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATc                                              >  1:2256240/1‑49 (MQ=255)
    cGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATc                                              <  2:2509291/49‑1 (MQ=255)
    cGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCt                                             <  2:1829374/50‑1 (MQ=255)
     ggCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCCAATCTa                                            <  2:2715823/50‑1 (MQ=255)
      gCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTAt                                           >  1:1369659/1‑50 (MQ=255)
       cACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATc                                          <  1:2370535/50‑1 (MQ=255)
       cACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATc                                          >  1:2606156/1‑50 (MQ=255)
       cACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATc                                          <  1:4631028/50‑1 (MQ=255)
       cACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATc                                          <  1:1101383/50‑1 (MQ=255)
       cACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCCAATCTATc                                          >  1:3138520/1‑50 (MQ=255)
        aCTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCg                                         >  1:5346376/1‑50 (MQ=255)
        aCTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCg                                         <  1:3734548/50‑1 (MQ=255)
        aCTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCg                                         <  2:4248530/50‑1 (MQ=255)
        aCTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCg                                         >  2:4387701/1‑50 (MQ=255)
        aCTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCg                                         >  2:3241339/1‑50 (MQ=255)
        aCTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCg                                         <  1:330765/50‑1 (MQ=255)
        aCTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCCAATCTATCg                                         >  2:3964792/1‑50 (MQ=255)
        aCTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCCAATCTATCg                                         >  2:1029482/1‑50 (MQ=255)
         cTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGc                                        >  1:4923530/1‑50 (MQ=255)
         cTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGc                                        >  2:284068/1‑50 (MQ=255)
         cTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGc                                        <  2:2540420/50‑1 (MQ=255)
         cTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCCAATCTATCGc                                        >  1:4350700/1‑50 (MQ=255)
         cTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCCAATCTATCGc                                        >  2:2137192/1‑50 (MQ=255)
         cTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCCAATCTATCGc                                        >  2:3463706/1‑50 (MQ=255)
          tGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCt                                       >  2:1785588/1‑50 (MQ=255)
          tGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCCAATCTATCGCt                                       >  1:204074/1‑50 (MQ=255)
          tGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCCAATCTATCGCt                                       >  1:4463417/1‑50 (MQ=255)
            ttatttatTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTg                                      >  1:3438647/1‑49 (MQ=255)
            ttatttatTCTCAGTATGATGAAAAAGCAGGCCAGCCAATCTATCGCTGa                                     <  1:3951713/50‑1 (MQ=255)
            ttatttatTCTCAGAATGATGAAAAAGCAGGCCAGCCAATCTATCGCTg                                      >  1:63114/1‑49 (MQ=38)
               tttattCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAaca                                  >  2:5023504/1‑50 (MQ=255)
               tttattCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAaca                                  <  2:3372818/50‑1 (MQ=255)
               tttattCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAac                                   <  1:2697942/49‑1 (MQ=255)
               tttattCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAac                                   <  1:4223564/49‑1 (MQ=255)
               tttattCTCAGTATGATGAAAAAGCAGGCCAGCCAATCTATCGCTGAaca                                  <  2:3890895/50‑1 (MQ=255)
               tttattCTCAGTATGATGAAAAAGCAGGCCAGCCAATCTATCGCTGAaca                                  >  1:78983/1‑50 (MQ=255)
                 tattCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAacac                                 <  1:241598/49‑1 (MQ=255)
                 tattCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAAcaca                                >  2:5346507/1‑50 (MQ=255)
                  attCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACAt                               >  1:2029186/1‑50 (MQ=255)
                  attCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACAt                               >  2:4709983/1‑50 (MQ=255)
                  attCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACAt                               >  2:3260312/1‑50 (MQ=255)
                  attCTCAGTATGAGGAAAAAGCAGGCCAGCTAATCTATCGCTGAAcaca                                <  1:3578385/49‑1 (MQ=255)
                   ttCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACGCAtt                              <  2:3915244/50‑1 (MQ=255)
                   ttCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACAtt                              <  1:694756/50‑1 (MQ=255)
                   ttCTCAGTATGATGAAAAAGCAGGCCAGCCAATCTATCGCTGAACACAtt                              >  2:703954/1‑50 (MQ=255)
                   ttCTCAGTATGATGAAAAAGCAGGCCAGCCAATCTATCGCTGAACACAtt                              >  2:3970670/1‑50 (MQ=255)
                    tctcAGTATGATGAAAAAGCAGGCCAGCCAATCTATCGCTGAACACAttt                             >  2:4897565/1‑50 (MQ=255)
                     ctcAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACAttt                             <  2:739771/49‑1 (MQ=255)
                     ctcAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTg                            >  2:2849576/1‑50 (MQ=255)
                     ctcAGTATGATGAAAAAGCAGGCCAGCCAATCTATCGCTGAACACATTTg                            >  1:3481318/1‑50 (MQ=255)
                        aGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTg                            <  2:4594308/47‑1 (MQ=255)
                        aGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCg                         >  1:4151067/1‑50 (MQ=255)
                        aGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCg                         >  2:3049638/1‑50 (MQ=255)
                        aGTATGATGAAAAAGCAGGCCAGCCAATCTATCGCTGAACACATTTg                            <  1:1266641/47‑1 (MQ=255)
                        aGTATGATGAAAAAGCAGGCCAGCCAATCTATCGCTGAACACATTTGTc                          <  2:360762/49‑1 (MQ=255)
                        aGTATGATGAAAAAGCAGGCCAGCCAATCTATCGCTGAACACATTTGTCg                         <  1:969990/50‑1 (MQ=255)
                        aGTATGATGAAAAAGCAGGCCAGCCAATCTATCGCTGAACACATTTGTCg                         <  1:942598/50‑1 (MQ=255)
                        aGTATGATGAAAAAGCAGGCCAGCCAATCTATCGCTGAACACATTTGTCg                         >  2:2810784/1‑50 (MQ=255)
                        aGTATGATGAAAAAGCAGGCCAGCCAATCTATCGCTGAACACATTTGTCg                         <  2:2966027/50‑1 (MQ=255)
                          tATGATGAAAAAGCAGGCCAGCCAATCTATCGCTGAACACATTTGTCGGa                       >  1:2191609/1‑50 (MQ=255)
                            tgatgaAAAAGCAGGCCAGCCAATCTATCGCTGAACACATTTGTCGGATg                     <  2:594255/50‑1 (MQ=255)
                              atgaAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATgcg                   >  1:4446220/1‑50 (MQ=255)
                              atgaAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATgcg                   >  1:3743575/1‑50 (MQ=255)
                               tgaAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATgcgg                  <  2:528195/50‑1 (MQ=255)
                               tgaAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATgcgg                  <  1:3191100/50‑1 (MQ=255)
                               tgaAAAAGCAGGCCAGCCAATCTATCGCTGAACACATTTGTCGGATgcgg                  >  1:526453/1‑50 (MQ=255)
                                 aaaaaGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATgcggcg                >  2:1360974/1‑50 (MQ=255)
                                 aaaaaGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATgcggcg                <  2:370358/50‑1 (MQ=255)
                                 aaaaaGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATgcggc                 <  2:1668727/49‑1 (MQ=255)
                                 aaaaaGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATgcggc                 <  1:1256016/49‑1 (MQ=255)
                                 aaaaaGCAGGCCAGCCAATCTATCGCTGAACACATTTGTCGGATgcggcg                >  1:4972937/1‑50 (MQ=255)
                                  aaaaGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGgcgc               <  1:5175791/50‑1 (MQ=255)
                                  aaaaGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGgcgc               >  2:3597760/1‑50 (MQ=255)
                                   aaaGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGgcgc               >  2:1274559/1‑49 (MQ=255)
                                    aaGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGgcgc               >  2:3107238/1‑48 (MQ=255)
                                    aaGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGgcgc               >  2:2662681/1‑48 (MQ=255)
                                    aaGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGa             <  1:2334446/50‑1 (MQ=255)
                                    aaGCAGGCCAGCTAATCTATCGCTGAACACATTTGCCGGATGCGgcgc               >  1:5121513/1‑48 (MQ=255)
                                    aaGCAGGCCAGCCAATCTATCGCTGAACACATTTGTCGGATGCGgcgc               >  2:1614476/1‑48 (MQ=255)
                                    aaGCAGGCCAGCCAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGa             <  2:2227960/50‑1 (MQ=255)
                                      gCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAgc           <  2:4531931/50‑1 (MQ=255)
                                      gCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAgc           <  1:337003/50‑1 (MQ=255)
                                        aGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAgcgc         >  2:3433170/1‑50 (MQ=255)
                                        aGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAgcgc         <  2:4951015/50‑1 (MQ=255)
                                        aGGCCAGCCAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAgcgc         >  1:1965250/1‑50 (MQ=255)
                                         ggCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGcc        <  1:404053/50‑1 (MQ=255)
                                         ggCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGcc        >  1:2717192/1‑50 (MQ=255)
                                          gccagccaATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCt       >  1:3894811/1‑50 (MQ=255)
                                            cAGCCAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTa     >  1:4365270/1‑50 (MQ=255)
                                             aGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTAt    >  2:4431838/1‑50 (MQ=255)
                                             aGCCAATCTATCGCTGAACACATTTGCCGGATGCGGCGCGAGCGCCTTAt    >  2:1749795/1‑50 (MQ=25)
                                              gCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATc   <  2:201440/50‑1 (MQ=255)
                                               cTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATcc  <  2:4329088/50‑1 (MQ=255)
                                               cTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATcc  <  2:2082909/50‑1 (MQ=255)
                                               cTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATcc  <  1:1111163/50‑1 (MQ=255)
                                                |                                                
TGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCC  >  NC_000913/1944149‑1944245

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