Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,355,9990GA33.8% 170.6 / 144.9 136W425* (TGG→TGAgltBglutamate synthase subunit GltB
Reads supporting (aligned to +/- strand):  ref base G (40/50);  new base A (22/24);  total (62/74)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 158 total aligned reads displayed.

GCAGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGATG  >  NC_000913/3355952‑3356048
                                               |                                                 
gCAGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtggat                                                 >  1:2680228/1‑50 (MQ=255)
gCAGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtggat                                                 <  1:1968069/50‑1 (MQ=255)
gCAGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtgaat                                                 >  2:2360260/1‑50 (MQ=255)
 cAGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAACGAGtgaatg                                                >  2:4851207/1‑50 (MQ=38)
 cAGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtggatg                                                <  2:4363964/50‑1 (MQ=255)
 cAGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtggatg                                                >  2:2721418/1‑50 (MQ=255)
  aGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtggatg                                                >  1:2912551/1‑49 (MQ=255)
  aGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtgaatgg                                               <  1:1255553/50‑1 (MQ=255)
  aGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtgaatgg                                               <  1:757022/50‑1 (MQ=255)
   gAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtggatg                                                >  1:4004712/1‑48 (MQ=255)
   gAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtggatg                                                >  1:3361319/1‑48 (MQ=255)
    aaaCCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtggatg                                                >  1:3511147/1‑47 (MQ=255)
    aaaCCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGAATGgag                                             <  2:1971678/50‑1 (MQ=255)
      aCCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGAATGgaga                                            >  2:2645847/1‑49 (MQ=255)
       ccGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGaaa                                          <  2:546922/50‑1 (MQ=255)
       ccGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGaaa                                          <  1:2862259/50‑1 (MQ=255)
       ccGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGAATGGAGaa                                           >  2:690102/1‑49 (MQ=255)
       ccGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGAATGGAGaa                                           >  2:3976835/1‑49 (MQ=255)
        cGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGaaaa                                         <  2:3100200/50‑1 (MQ=255)
        cGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGaaa                                          >  2:5048490/1‑49 (MQ=255)
         gATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGaaaaa                                        >  1:3794605/1‑50 (MQ=255)
         gATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGaaaaa                                        >  2:2200456/1‑50 (MQ=255)
         gATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGAATGGAGaaaaa                                        >  1:2778516/1‑50 (MQ=255)
         gATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGAATGGAGaaaaa                                        <  1:5084613/50‑1 (MQ=255)
          aTGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAAc                                       <  1:1559175/50‑1 (MQ=255)
          aTGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAAc                                       <  1:4696304/50‑1 (MQ=255)
          aTGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAAc                                       <  2:2910054/50‑1 (MQ=255)
          aTGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAAc                                       <  2:4457782/50‑1 (MQ=255)
             aCGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGt                                     <  1:921381/49‑1 (MQ=255)
              cGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTcc                                   <  1:650945/50‑1 (MQ=255)
              cGATCTGAAAAACCGCCATGCATATAAAGAGAGAATGTAGAAAAACGTcc                                   <  2:3674926/50‑1 (MQ=255)
               gATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTccg                                  >  2:789761/1‑50 (MQ=255)
                aTCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTccgc                                 >  2:3696564/1‑50 (MQ=255)
                 tCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTccgcc                                >  2:2337978/1‑50 (MQ=255)
                  cTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTccgcc                                >  1:2536629/1‑49 (MQ=255)
                  cTGAAAAGCCGCCATCCATATAAAGAGTGAATGGAGAAAAACGTccgcc                                >  1:3730993/1‑49 (MQ=255)
                   tGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGa                              >  2:4993290/1‑50 (MQ=255)
                   tGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGa                              >  1:3139916/1‑50 (MQ=255)
                   tGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGa                              >  1:5195561/1‑50 (MQ=255)
                    gAAAAGCCGCCATCCATATAAAGAGTGAATGGAGAAAAACGTCCGCCGAc                             >  1:3644231/1‑50 (MQ=255)
                     aaaaGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTccgccg                               <  2:1195790/47‑1 (MQ=255)
                     aaaaGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTccgcc                                <  2:2173833/46‑1 (MQ=255)
                     aaaaGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTccgc                                 <  2:3924431/45‑1 (MQ=255)
                     aaaaGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTccg                                  <  1:4159062/44‑1 (MQ=255)
                     aaaaGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGa                              <  1:3373958/48‑1 (MQ=255)
                     aaaaGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGa                              <  1:1583924/48‑1 (MQ=255)
                     aaaaGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGACt                            <  1:316710/50‑1 (MQ=255)
                     aaaaGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGACt                            <  2:1929681/50‑1 (MQ=255)
                     aaaaGCCGCCATCCATATAAAGAGTGAATGGAGAAAAACGTCCGCCGACt                            <  1:2320082/50‑1 (MQ=255)
                      aaaGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTg                           <  1:1889455/50‑1 (MQ=255)
                      aaaGCCGCCATCCATATAAAGAGTGAATGGAGAAAAACGTCCGCCGACTg                           >  1:393840/1‑50 (MQ=255)
                       aaGCCGCCATCCATATAAAGAGTGAATGGAGAAAAACGTCCGCCGACTgg                          <  1:3201071/50‑1 (MQ=255)
                        aGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTgg                          >  1:5060138/1‑49 (MQ=255)
                        aGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTgg                          >  1:661274/1‑49 (MQ=255)
                        aGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGt                         <  1:1895701/50‑1 (MQ=255)
                         gccgccATCCATATAAAGAGTGAATGGAGAAAAACGTCCGCCGACTGGTa                        <  1:2445497/50‑1 (MQ=255)
                          ccgccATCCATATAAAGAGTGAATGGAGAAAAACGTCCGCCGACTGGTAc                       <  2:3076305/50‑1 (MQ=255)
                           cgccATCCATATAAAGAGTGAATGGAGAAAAACGTCCGCCGACTGGTAcc                      >  1:1069122/1‑50 (MQ=255)
                           cgccATCCATATAAAGAGTGAATGGAGAAAAACGTCCGCCGACTGGTAcc                      >  2:4237562/1‑50 (MQ=255)
                            gccATCCATATAAAGAGTGAATGGAGAAAAACGTCCGCCGACTGGTACCg                     >  2:1425974/1‑50 (MQ=255)
                            gccATCCATATAAAGAGTGAATGGAGAAAAACGTCCGCCGACTGGTACCg                     >  2:3461282/1‑50 (MQ=255)
                             ccatccatATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGt                    >  1:5050838/1‑50 (MQ=255)
                             ccatccatATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGt                    >  1:4281028/1‑50 (MQ=255)
                             ccatccatATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGt                    >  2:2989734/1‑50 (MQ=255)
                              catccatATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGtt                   <  1:1800883/50‑1 (MQ=255)
                              catccatATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGtt                   >  2:717923/1‑50 (MQ=255)
                               atccatATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGttt                  >  2:5171302/1‑50 (MQ=255)
                               atccatATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGttt                  <  1:426343/50‑1 (MQ=255)
                                tccatATAAAGAGTGAATGGAGAAAAACGTCCGCCGACTGGTACCGttt                  >  2:5139081/1‑49 (MQ=255)
                                 ccatATAAAGAGTGAATGGAGAAAAACGTCCGCCGACTGGTACCGTTTg                 >  1:4136067/1‑49 (MQ=255)
                                  catATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTg                 >  2:5198442/1‑48 (MQ=255)
                                  catATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTg                 >  2:3096389/1‑48 (MQ=255)
                                   atatAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTg                 >  2:4497427/1‑47 (MQ=255)
                                   atatAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTg                 >  2:4586511/1‑47 (MQ=255)
                                   atatAAAGAGTGAATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAg              <  1:251957/50‑1 (MQ=255)
                                   atatAAAGAGGGGATGGAGAAAAACGTTCGCCGACTGGTTCCGTTTGAAg              <  2:3543301/50‑1 (MQ=255)
                                     ataAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTg                 >  1:3573645/1‑45 (MQ=255)
                                      taAAGTGTGGCTGGAGCAAAACGTCCGCCGACTAGTACCGTTTGAAGATc           >  2:4020800/1‑50 (MQ=255)
                                      taAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGAt            <  2:1370190/49‑1 (MQ=255)
                                      taAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGAt            <  1:4985756/49‑1 (MQ=255)
                                      taAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATc           <  2:1626354/50‑1 (MQ=255)
                                      taAAGAGTGAATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGAt            <  2:1857331/49‑1 (MQ=255)
                                        aaGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCt          <  1:291480/49‑1 (MQ=255)
                                        aaGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTg         <  1:702219/50‑1 (MQ=255)
                                        aaGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTg         <  2:5344458/50‑1 (MQ=255)
                                        aaGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTg         >  2:4016607/1‑50 (MQ=255)
                                        aaGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTg         <  1:1701149/50‑1 (MQ=255)
                                         agagTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGc        <  2:1189914/50‑1 (MQ=255)
                                         agagTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGc        >  1:521686/1‑50 (MQ=255)
                                          gagTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGcc       >  2:3190460/1‑50 (MQ=255)
                                          gagTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGcc       >  2:2768335/1‑50 (MQ=255)
                                           agTGGATGGAGAAAAACGTCCGCCGTCTGGTACCGTTTGAAGATCTGccc      >  2:2835974/1‑50 (MQ=255)
                                           agTGAATGGAGAAAAACGTCCGCCGACTGGCACCGTTTGAAGATCTGccc      >  2:4721675/1‑50 (MQ=39)
                                            gTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGccc      >  1:2212020/1‑49 (MQ=255)
                                            gTGAATGGAGAAAAACGTCCGTCGACTGGTACCGTTTGAAGATCTGCCCg     <  2:3302219/50‑1 (MQ=39)
                                            gTGAATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGccc      >  1:4338106/1‑49 (MQ=255)
                                            gTGAATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCg     <  1:1098603/50‑1 (MQ=255)
                                              ggatggaGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGAt   >  1:1207053/1‑50 (MQ=255)
                                               gatggaGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGAt   >  2:4348835/1‑49 (MQ=255)
                                               gatggaGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGATg  <  2:1758408/50‑1 (MQ=255)
                                               |                                                 
GCAGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGATG  >  NC_000913/3355952‑3356048

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