Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,908,5260AG27.4% 127.0 / 16.1 62intergenic (‑160/+23)pstB/pstAphosphate ABC transporter ATP binding subunit/phosphate ABC transporter membrane subunit PstA
Reads supporting (aligned to +/- strand):  ref base A (32/13);  new base G (0/17);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.62e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAA  >  NC_000913/3908477‑3908574
                                                 |                                                 
aaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGa                                                   <  1:5304857/50‑1 (MQ=255)
 aaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                  <  1:772486/50‑1 (MQ=255)
  aGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAgg                                                 >  1:1174155/1‑50 (MQ=255)
  aGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                 <  1:2422530/50‑1 (MQ=255)
  aGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                 <  1:2136341/50‑1 (MQ=255)
   ggCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg                                                <  1:5322885/50‑1 (MQ=255)
   ggCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg                                                <  1:1995152/50‑1 (MQ=255)
   ggCGCTCGGTCGGTCACCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg                                                >  1:2367687/1‑50 (MQ=255)
    gcgcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg                                                >  1:2895262/1‑49 (MQ=255)
     cgcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgggg                                               >  1:4411181/1‑49 (MQ=255)
       cTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGCGaa                                            >  1:4432698/1‑49 (MQ=255)
          ggtcggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAcc                                        >  1:4269662/1‑50 (MQ=40)
          ggtcggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAcc                                        >  1:3950943/1‑50 (MQ=40)
           gtcggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCg                                       >  1:1840680/1‑50 (MQ=40)
             cggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCgtg                                     >  1:6023329/1‑50 (MQ=255)
              ggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCgtgt                                    >  1:5966340/1‑50 (MQ=255)
               gtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCgtgt                                    >  1:2946545/1‑49 (MQ=255)
                     tCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCaa                              >  1:1824762/1‑49 (MQ=255)
                     ccGCCCCTTGGGGGGGGGGGTTGGGGGGGGGGGAAAACCGTGTCAGCAat                             <  1:6756278/49‑1 (MQ=12)
                      cGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAata                            >  1:6691655/1‑50 (MQ=255)
                       gCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAatat                           >  1:5564766/1‑50 (MQ=255)
                       gCCCCTCTGGGGGGGGGTTGGGGGGGGGGGAAAAACCGTGTCAGCAatat                           <  1:1818865/50‑1 (MQ=12)
                       gCCCCTCGGGGGAGGGGGTTGGGGGGGGGGGAAAACCGTGTCAGCAatat                           <  1:1512445/50‑1 (MQ=21)
                        ccccTGGGGGGAGAGGGTTGGGGTGGGGGGAAAACCGTGTCAGCAATATc                          <  1:6867385/50‑1 (MQ=21)
                         cccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATc                          >  1:700110/1‑49 (MQ=255)
                         cccTCGGGGGGGAGGTTTAGGGTGGGGGGAAAACCGTGTCAGCAATATCa                         <  1:5305188/50‑1 (MQ=37)
                         cccTCGGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCa                         <  1:4838164/50‑1 (MQ=255)
                         cccTCGGGGAAGAGGGTGGGGGAGGGGGAAAAACCGTGTCAGCAATATCa                         <  1:2120252/50‑1 (MQ=255)
                          cttctGGGGAGGGGTTTAGGGGGAGGGGAAAACCGTGTCAGCAATATCaa                        <  1:5741696/48‑1 (MQ=25)
                           ctctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                       >  1:1243778/1‑50 (MQ=255)
                           ctcgGGGGGGGGGGTAGGGGTGGGGGGAAAACCGTGTCAGCAATATCAAc                       <  1:4927742/50‑1 (MQ=255)
                            tctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAcc                      >  1:921158/1‑50 (MQ=255)
                              tGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgt                    >  1:2362623/1‑50 (MQ=255)
                              gggggAGGGGGTGGGGGGGGGGGGAAAACCGTGTCAGCAATATCAACCgt                    <  1:4303276/49‑1 (MQ=255)
                              gggggAGAGGGTTGGGGTGGGGGGAAAACCGTGTCAGCAATATCAACCg                     <  1:5415176/48‑1 (MQ=37)
                                gggAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgtgt                  >  1:2373689/1‑50 (MQ=255)
                                gggAGAGGGTAGGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgtg                   <  1:2466003/49‑1 (MQ=37)
                                 ggAGGGGTTTGGGGTGGGGGGAAAACCGTGTCAGCAATATCAACCGTGtt                 <  1:4888258/50‑1 (MQ=255)
                                 ggAGAGGGTTGGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGtt                 <  1:593756/50‑1 (MQ=255)
                                 ggAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGtt                 >  1:5783674/1‑50 (MQ=255)
                                 ggAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGtt                 >  1:3744899/1‑50 (MQ=255)
                                  gagaGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGtt                 >  1:4474224/1‑49 (MQ=255)
                                  gagaGGGTAAGGGGGGGGGGAAAACCGTGTCAGCAATATCAACCGTGttt                <  1:5834866/50‑1 (MQ=38)
                                   agagGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTa               >  1:7281831/1‑50 (MQ=255)
                                   agagGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTa               >  1:7450684/1‑50 (MQ=255)
                                   agagGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTa               >  1:1292278/1‑50 (MQ=255)
                                   agagGGTGGGGGGGGGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTa               <  1:5832734/50‑1 (MQ=255)
                                    gagGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAt              >  1:5138177/1‑50 (MQ=255)
                                      gggTTAGGGGGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttc            <  1:1448772/50‑1 (MQ=255)
                                       ggTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttct           >  1:6738654/1‑50 (MQ=255)
                                       ggTGAGGGGGAGGGAAAAACCGTGTCAGCAATATCAACCGTGTTTAttct           <  1:5415071/50‑1 (MQ=255)
                                        gTTGGGGTGGGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttctt          <  1:3480655/50‑1 (MQ=255)
                                        gTTGGGGTGGGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttctt          <  1:1347526/50‑1 (MQ=255)
                                          tAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg        >  1:3223357/1‑50 (MQ=255)
                                          tAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg        >  1:2974169/1‑50 (MQ=255)
                                           ggggTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg        <  1:561749/48‑1 (MQ=255)
                                           aGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg        >  1:2369542/1‑49 (MQ=255)
                                            gggTGGGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCa      <  1:4711295/50‑1 (MQ=255)
                                             ggTGGGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaa     <  1:418770/50‑1 (MQ=255)
                                             ggTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCa      >  1:2494984/1‑49 (MQ=255)
                                              gTGGGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaa    <  1:642722/50‑1 (MQ=255)
                                              gTGAGGGAAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCgaa    <  1:200673/50‑1 (MQ=255)
                                               tGAGGGGAAAACCGTGTCAGCAATATCAACTGTGTTTATTCTTCGCaaa    >  1:6467316/1‑49 (MQ=255)
                                                 aGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaaaa  <  1:2089604/49‑1 (MQ=255)
                                                 |                                                 
AAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAA  >  NC_000913/3908477‑3908574

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