Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091329,4860CA48.9% 47.5 / 34.5 45intergenic (+291/‑165)dapB/carAdihydrodipicolinate reductase/carbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base C (17/6);  new base A (0/22);  total (17/28)
Fisher's exact test for biased strand distribution p-value = 1.15e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATT  >  NC_000913/29340‑29636
                                                                                                                                                   |                                                                                                                                                      
tgttgtTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACaaa                                                                                                                                                     <  1:272580/150‑1 (MQ=255)
 gttgttTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACaaaa                                                                                                                                                    <  1:927707/150‑1 (MQ=255)
 gttgttTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACaaaa                                                                                                                                                    <  1:201255/150‑1 (MQ=255)
 gttgttTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACaaaa                                                                                                                                                    <  1:295247/150‑1 (MQ=255)
     ttttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAaataat                                                                                                                                                >  1:172361/1‑150 (MQ=255)
       ttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAatat                                                                                                                                              >  2:489701/1‑150 (MQ=255)
            gTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaaaa                                                                                                                                         <  2:798582/150‑1 (MQ=255)
              aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATc                                                                                                                                      >  1:416374/1‑151 (MQ=255)
                  tttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCg                                                                                                                                   <  2:65160/150‑1 (MQ=255)
                  tttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGc                                                                                                                                  >  2:600533/1‑151 (MQ=255)
                          tttGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAg                                                                                                                           >  1:268928/1‑150 (MQ=255)
                               tCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACt                                                                                                                     >  1:357300/1‑151 (MQ=255)
                                     tttttCTCCTTTTTTTTTTTTCATGCAATTTTTTTGCGGCAAAGGCTTTTCCAAAAAAAGTTTAAATGATTTTTTTTTCATTAAATGCTTGAAAAACATGCATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAg                                                                                                                <  1:697536/150‑1 (MQ=255)
                                           tGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATACCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCat                                                                                                         >  1:31695/1‑151 (MQ=255)
                                                     tATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGa                                                                                                >  1:153646/1‑150 (MQ=255)
                                                       tttCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATg                                                                                             >  2:912259/1‑151 (MQ=255)
                                                             gCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGccgccg                                                                                       >  1:823408/1‑151 (MQ=255)
                                                                  cttcttGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGc                                                                                  >  1:3486/1‑151 (MQ=255)
                                                                   tttttGCTGCGCAAGCTTTTACAAAAACAGGTAAAATTATTTTTTTGTCTTTTAATGCCTAAAAAAAATGCATGACCCAAAAAATAAAAAAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgcc                                                                                        <  1:518551/141‑1 (MQ=255)
                                                                        gccgcgaAAGCTTTTTCCAAAAAGGTTTAGTTAATTTTTTTGCCGTTAAGGCCTTTAAAAACTTGAAGGAGCAAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGaa                                                                             <  1:122465/147‑1 (MQ=255)
                                                                             gcAGGCTTTTCCCAAAACAGGTTAGATGTTTTTTTTGTCGTTTAATGCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCa                                                                                <  2:210793/142‑1 (MQ=255)
                                                                               aaGCTTTTTCAAGAACAGTTTAAATTATTTTTTTGTCGTTTAATGCTTAAAAAACATGCATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcg                                                                        <  1:385747/148‑1 (MQ=255)
                                                                               aaGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg                                                                     >  2:488450/1‑151 (MQ=255)
                                                                                    ttttCCAAAACAGGTTAGATGTTTTTTTTGTCGCTTAATGCCTTAAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGt                                                                   <  2:104531/148‑1 (MQ=255)
                                                                                    tttacaaGAACAGGTAAGATATTTTTTTTGCCGTTTATTCCTTAAAAACCTTGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtc                                                                      <  1:175430/139‑1 (MQ=255)
                                                                                     tttCCAGAATGGGTTAGATGTTTTTTTTGGCGATAAATGCCTGAAAAACATGAATGAGCAAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAg                                                                <  1:742165/150‑1 (MQ=255)
                                                                                           aaaCAGTTTAGATGTTTTTTTTGTCGTTTAATGCTTTAAAAACATGCAGGACCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAttt                                                         <  2:131617/150‑1 (MQ=255)
                                                                                            aaCAGGTTAGATGATTTTTTTGTCGCTTAATGCTTGTAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTg                                                        <  2:416384/151‑1 (MQ=255)
                                                                                            aaCAGGTAAGATAATTTTTTTTTCTCTTAATCCTTATAAAACATGATTGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAttt                                                         <  1:842319/150‑1 (MQ=255)
                                                                                                gttaagttgatTTTTTTGTGGTTTAAGGCGTGAAAAACATGCATGACCAAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCa                                                      <  2:876911/142‑1 (MQ=255)
                                                                                                 gTTAGAAGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCAGTTTGCCAGaaa                                                                            >  2:933968/1‑126 (MQ=255)
                                                                                                  ttAGATGTTTTTTTTGTCGTTAAATGCTTGAAAAACATGCATGACCCAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCa                                                      <  1:409301/147‑1 (MQ=255)
                                                                                                      atgatCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAttt                                               >  2:704531/1‑150 (MQ=255)
                                                                                                             ttttttCGTTAAATCCCTAAAAAACATGAATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcat                                        <  1:600521/150‑1 (MQ=255)
                                                                                                               ttGTCCCTTATTTCTTGTAACATTTCCATGGCCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatca                                      <  2:72758/150‑1 (MQ=255)
                                                                                                                 gTCTCTTATCGCTTGAAAAACAGGCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCAtt                                    <  1:586077/150‑1 (MQ=255)
                                                                                                                   cGCTTAATGCCTGAAAACCATGCATGCGCCAAAAAATAATAAAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgt                                  <  2:263727/150‑1 (MQ=255)
                                                                                                                       ttattccctaaaaaaCATCCAGGAGCAAAAAAAAAAAAAAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGa                                <  1:213872/138‑1 (MQ=255)
                                                                                                                              tGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGc                      >  1:315197/1‑151 (MQ=255)
                                                                                                                                  aaaaCATGCATGACCAAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                       <  2:440701/147‑1 (MQ=255)
                                                                                                                                         gcatgACCAAAAAAATAATAAAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg            <  1:527241/151‑1 (MQ=255)
                                                                                                                                            ttacccAAAAAATAATAAAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAag              <  2:172366/142‑1 (MQ=255)
                                                                                                                                             gacccAAAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagt         <  2:883407/147‑1 (MQ=255)
                                                                                                                                              aGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagtg        >  2:435049/1‑150 (MQ=255)
                                                                                                                                              aGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagtg        >  1:842759/1‑150 (MQ=255)
                                                                                                                                                  acaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat   >  2:185178/1‑151 (MQ=255)
                                                                                                                                                  acaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAtt  >  1:311059/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATT  >  NC_000913/29340‑29636

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