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

TTTTTGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAT  >  NC_000913/29336‑29635
                                                                                                                                                         |                                                                                                                                                     
tttttGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCcac                                                                                                                                                       >  2:124369/1‑151 (MQ=255)
         ttttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAaataat                                                                                                                                               >  2:661046/1‑150 (MQ=255)
          tttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAataata                                                                                                                                              <  1:1093/150‑1 (MQ=255)
              aTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaa                                                                                                                                          >  1:577989/1‑150 (MQ=255)
                gTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaaaa                                                                                                                                        >  2:18715/1‑150 (MQ=255)
                    aTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATccc                                                                                                                                   >  2:433988/1‑151 (MQ=255)
                    aTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATccc                                                                                                                                   >  1:613654/1‑151 (MQ=255)
                      tttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCg                                                                                                                                  <  1:124332/150‑1 (MQ=255)
                                 ggTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTg                                                                                                                       >  2:793037/1‑150 (MQ=255)
                                    cctcTTTTTTCTCCTTTTTTTTATTTCATCCAATCTTCTTGCTGCGCAAGCTTTTCCAAAAACGGTTTAAATTATTTTTTTGTCCCTAAATCCCTTAAAAACATGCATGACCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGAc                                                                                                                     <  2:139244/146‑1 (MQ=255)
                                    ccATTTTTTTCTTTTTTTTTTTATTTCATGCAATCTTCTTGCTGCGAAAGCTTTTCCAAAAACGGTTAAAATAATTTTTTTGTCTTTTAATCCCTTTAAAACATGCATGAGCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACt                                                                                                                    <  1:282949/150‑1 (MQ=255)
                                        ttttttttGCTCTTTTTTTTTTAATGCATTTTTTTTGCTCCAAATGCTTTTTCAAAAAAAGGTAAGATAATTTTTTTGTCGTTTATTCCCGAAAAAAGTTGCTTGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACtttt                                                                                                                 <  1:169831/149‑1 (MQ=255)
                                          ttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgcg                                                                                                             >  2:52126/1‑151 (MQ=255)
                                                    gTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                   >  1:279304/1‑151 (MQ=255)
                                                     tttttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                   >  1:754187/1‑150 (MQ=255)
                                                      ttttATTTCATGCATTTTTTTTGCTGCGAAAGCTTTTCCAAAAACAGGTTAGATATTTTTTTTGCCGTTTATTCCCTAAAAAACATGAATGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                  <  1:784151/150‑1 (MQ=255)
                                                             tCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgcc                                                                                          >  2:780930/1‑151 (MQ=255)
                                                               aTGCATTCTTCTTGCGGCGAAGGCTTTTACAAAAAAGTGTTAAAAAATTTTTTTGACGTTTAGGCCCTAAAAAACGTGTTGGAGCCAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgcc                                                                                          <  1:660987/149‑1 (MQ=255)
                                                                tGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgcc                                                                                       >  2:231453/1‑151 (MQ=255)
                                                                       ttcttGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                >  2:132716/1‑151 (MQ=255)
                                                                            gctgcGCAAGCGTTTTCCAGTACAGGTTAGATGATCTTTTTGTGGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGa                                                                             >  1:96022/1‑149 (MQ=255)
                                                                              tgcGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAtt                                                                         >  1:359277/1‑151 (MQ=255)
                                                                                   aagcttttcaaaaaaCAGTTTAAATATTTTTTTTGTCGTTAATTCCCTGAAAAACATTCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtc                                                                     <  1:158349/137‑1 (MQ=255)
                                                                                   aaGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg                                                                    >  1:750671/1‑151 (MQ=255)
                                                                                   aaGCGTTTCCAAAAACGGGTTAAATGATCTTTTTGTCGTTTAATGCCTGAAAAACATGCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgt                                                                      <  2:565673/149‑1 (MQ=255)
                                                                                         tttCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                              >  1:94899/1‑151 (MQ=255)
                                                                                          ttCCAGAACAGGTTAGATGATCTTTTTGGCGCTTAAAGCCTGTAACACCTGCATAGGCCACCAAATACTCTAAAAAATCCCGCCATTAAGTTAACTGTTAGCGCCCATATCTCCCGAATGCCGCCGTTTGCAACAAATTCGTCGGTAAGc                                                              >  1:529128/1‑150 (MQ=255)
                                                                                            ccAAAACAGGTAAGTTGATTTTTTGTTCTTTTATGGCCTTAAAAAATTGCAGGAACAAAAAAATAAAAAAAAACACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCa                                                             <  1:354375/149‑1 (MQ=255)
                                                                                                          aaattatttttttGTCGCTTAATGCTTAAAAAACAGGAATGAGCCAAAAAATAAAAAAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGa                                                 <  1:106500/141‑1 (MQ=255)
                                                                                                            atgtttttttttGCGCTAAATGCCTTAAAAACAGGTAGGAGCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGa                                                 <  2:485951/141‑1 (MQ=255)
                                                                                                             taatttttttGGCGAATATTGCCTTAAAAACAAGCAGGACCAAAAAAAAAAAATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTa                                             <  1:708301/145‑1 (MQ=255)
                                                                                                                tttttttGTCATTTATTCCTTGAAAAACATGAATGAGCAAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGt                                          <  2:790300/148‑1 (MQ=255)
                                                                                                                tCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGTAGATTTGCATTGATTTACGtc                                         >  2:423739/1‑151 (MQ=255)
                                                                                                                 ttttttGTCGTTTAATCCTTGAAAAACTTGCATGACCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                         <  2:614755/149‑1 (MQ=255)
                                                                                                                  tttttGTCGTTTAATGCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtca                                        <  2:776181/150‑1 (MQ=255)
                                                                                                                    tttGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatc                                      >  2:455373/1‑150 (MQ=255)
                                                                                                                      tGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCAtt                                   >  1:106089/1‑151 (MQ=255)
                                                                                                                                     ttAAAACTTTCATTGACCCAAAAAATAATATAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                      <  1:45468/149‑1 (MQ=255)
                                                                                                                                       aaacaatgcttgAGCAAAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaaa                  <  1:11599/146‑1 (MQ=255)
                                                                                                                                          aCAGGCAGGACCAACAAAATAATAAAAAAAAACCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATa                <  1:540004/150‑1 (MQ=255)
                                                                                                                                           catgcatgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaa               >  2:435888/1‑150 (MQ=255)
                                                                                                                                           catgcatgAACCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaa               <  2:105881/150‑1 (MQ=255)
                                                                                                                                               catgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg           >  1:680790/1‑150 (MQ=255)
                                                                                                                                                atgACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtga          <  2:297622/150‑1 (MQ=255)
                                                                                                                                                      ccACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGTTTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAgtgagtga      >  1:511518/1‑149 (MQ=255)
                                                                                                                                                      ccACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat    >  2:301258/1‑151 (MQ=255)
                                                                                                                                                        acaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat  >  2:681467/1‑151 (MQ=255)
                                                                                                                                                         |                                                                                                                                                     
TTTTTGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAT  >  NC_000913/29336‑29635

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