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

TGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATT  >  NC_000913/29340‑29636
                                                                                                                                                   |                                                                                                                                                      
tgttgtTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCATGAGCCACaaa                                                                                                                                                     <  1:272580/150‑1 (MQ=255)
 gttgttTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCATGAGCCACaaaa                                                                                                                                                    <  1:927707/150‑1 (MQ=255)
 gttgttTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCATGAGCCACaaaa                                                                                                                                                    <  1:201255/150‑1 (MQ=255)
 gttgttTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCATGAGCCACaaaa                                                                                                                                                    <  1:295247/150‑1 (MQ=255)
     ttttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCATGAGCCACAAaataat                                                                                                                                                >  1:172361/1‑150 (MQ=255)
       ttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCATGAGCCACAAAATAatat                                                                                                                                              >  2:489701/1‑150 (MQ=255)
            gTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCATGAGCCACAAAATAATATaaaaa                                                                                                                                         <  2:798582/150‑1 (MQ=255)
              aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATc                                                                                                                                      >  1:416374/1‑151 (MQ=255)
                  tttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCg                                                                                                                                   <  2:65160/150‑1 (MQ=255)
                  tttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGc                                                                                                                                  >  2:600533/1‑151 (MQ=255)
                          tttGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAg                                                                                                                           >  1:268928/1‑150 (MQ=255)
                               tCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACt                                                                                                                     >  1:357300/1‑151 (MQ=255)
                                     tttttCTCCTTTTTTTTTTTTCATGCAATTTTTTTGCGGCAAAGGCTTTTCCAAAAAAAGTTTAAATGATTTTTTTTTCATTAAATGCTTGA‑AAAACATGCATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAg                                                                                                                <  1:697536/150‑1 (MQ=255)
                                           tGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATACCTGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCat                                                                                                         >  1:31695/1‑151 (MQ=255)
                                                     tATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGa                                                                                                >  1:153646/1‑150 (MQ=255)
                                                       tttCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATg                                                                                             >  2:912259/1‑151 (MQ=255)
                                                             gCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGccgccg                                                                                       >  1:823408/1‑151 (MQ=255)
                                                                  cttcttGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGc                                                                                  >  1:3486/1‑151 (MQ=255)
                                                                   tttttGCTGCGCAAGCTTTTACAAAAACAGGTAAAATTATTTTTTTGTCTTTTAATGCCTAA‑AAAAAATGCATGACCCAAAAAATAAAAAAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgcc                                                                                        <  1:518551/141‑1 (MQ=255)
                                                                        gccgcgaAAGCTTTTTCCAAAAAGGTTTAGTTAATTTTTTTGCCG‑TTAAGGCCTTTAAAAACTTGAAGGAGCAAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGaa                                                                             <  1:122465/147‑1 (MQ=255)
                                                                             gcAGGCTTTTCCCAAAACAGGTTAGATGTTTTTTTTGTCGTTTAATGCCTGA‑AAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCa                                                                                <  2:210793/142‑1 (MQ=255)
                                                                               aaGCTTTTTCAAGAACAGTTTAAATTATTTTTTTGTCGTTTAATGCTTAA‑AAAACATGCATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcg                                                                        <  1:385747/148‑1 (MQ=255)
                                                                               aaGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg                                                                     >  2:488450/1‑151 (MQ=255)
                                                                                    ttttCCAAAACAGGTTAGATGTTTTTTTTGTCGCTTAATGCCTTA‑AAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGt                                                                   <  2:104531/148‑1 (MQ=255)
                                                                                    tttacaaGAACAGGTAAGATATTTTTTTTGCCGTTTATTCCTTAA‑AAACCTTGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtc                                                                      <  1:175430/139‑1 (MQ=255)
                                                                                     tttCCAGAATGGGTTAGATGTTTTTTTTGGCGATAAATGCCTGA‑AAAACATGAATGAGCAAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAg                                                                <  1:742165/150‑1 (MQ=255)
                                                                                           aaaCAGTTTAGATGTTTTTTTTGTCGTTTAATGCTTTA‑AAAACATGCAGGACCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAttt                                                         <  2:131617/150‑1 (MQ=255)
                                                                                            aaCAGGTTAGATGATTTTTTTGTCGCTTAATGCTTGT‑AAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTg                                                        <  2:416384/151‑1 (MQ=255)
                                                                                            aaCAGGTAAGATAATTTTTTTTTCTCTTAATCCTTAT‑AAAACATGATTGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAttt                                                         <  1:842319/150‑1 (MQ=255)
                                                                                                gttaagttgatTTTTTTGTGGTTTAAGGCGTGA‑AAAACATGCATGACCAAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCa                                                      <  2:876911/142‑1 (MQ=255)
                                                                                                 gTTAGAAGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCAGTTTGCCAGaaa                                                                            >  2:933968/1‑126 (MQ=255)
                                                                                                  ttAGATGTTTTTTTTGTCGTTAAATGCTTGA‑AAAACATGCATGACCCAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCa                                                      <  1:409301/147‑1 (MQ=255)
                                                                                                      atgatCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAttt                                               >  2:704531/1‑150 (MQ=255)
                                                                                                             ttttttCGTTAAATCCCTAA‑AAAACATGAATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcat                                        <  1:600521/150‑1 (MQ=255)
                                                                                                               ttGTCCCTTATTTCTTGT‑AACATTTCCATGGCCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatca                                      <  2:72758/150‑1 (MQ=255)
                                                                                                                 gTCTCTTATCGCTTGA‑AAAACAGGCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCAtt                                    <  1:586077/150‑1 (MQ=255)
                                                                                                                   cGCTTAATGCCTGA‑AAACCATGCATGCGCCAAAAAATAATAAAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgt                                  <  2:263727/150‑1 (MQ=255)
                                                                                                                       ttattccctaaaaaaCATCCAGGAGCAAAAAAAAAAAAAAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGa                                <  1:213872/138‑1 (MQ=255)
                                                                                                                              tGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGc                      >  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)
                                                                                                                                                  acaAAATAATATAAAAAATCCCGCCATTAAGTTGAC‑TTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAtt  >  1:311059/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATT  >  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: ‑