Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091329,4960TA44.7% 50.3 / 19.0 38intergenic (+301/‑155)dapB/carAdihydrodipicolinate reductase/carbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base T (18/3);  new base A (0/17);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 4.02e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAT  >  NC_000913/29346‑29635
                                                                                                                                                      |                                                                                                                                            
tttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAatat                                                                                                                                              >  1:375991/1‑151 (MQ=255)
 ttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAatat                                                                                                                                              <  2:602791/150‑1 (MQ=255)
  tAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAAtata                                                                                                                                             >  2:121768/1‑150 (MQ=255)
   aaTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaa                                                                                                                                           <  1:245173/151‑1 (MQ=255)
      gTAAATTTTGACCATTTGGTCCACTTTTTTCTGCACGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaaaa                                                                                                                                         <  2:129568/150‑1 (MQ=255)
        aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAAt                                                                                                                                       <  2:155964/150‑1 (MQ=255)
                aCCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAAT‑CCCGCCAtt                                                                                                                             >  1:310942/1‑151 (MQ=255)
                  cATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAAT‑CCCGCCATTa                                                                                                                            >  1:65845/1‑150 (MQ=255)
                    tttGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAAT‑CCCGCCATTAAGt                                                                                                                         >  2:102566/1‑151 (MQ=255)
                      tGGTCCATTTTTTTCTGCTTTTTTTTTTTTAATACATTTTTCTTGCTGCGAAAGCGTTTTCAAGAACAGTTTAAATAATTTTTTTGTCGTTTAATCCCTAAAAAACATGCATGAGCCAAAAAATAATAAAAAAAAT‑CCCGCCATTAAGtt                                                                                                                        <  1:614794/150‑1 (MQ=255)
                        gTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAACAT‑CCCGCCATTAAGTTGAc                                                                                                                     >  2:535321/1‑151 (MQ=255)
                             tttttttttGTTTTTTTTTATTTCATGAAATTTTTTTGCTGCGAAAGTTTTTCCAAAAAAAGGTAAAATAATTTTTTTGTCCTTTAATCCATAAAAAAATTGAATGGCCCAAAAAATAATAAAAAAAAT‑CCCGCCATTAAGTTGACttt                                                                                                                  <  1:501146/148‑1 (MQ=255)
                                   tctgttttttttttttTAATGCATTTTTCTTGCTGCGAAAGTTTTTCCAAAAAAAGGTTAGATGATTTTTTTGTCGTTTAATCCCTAAAAAACTTGCTTGACCCAAAAAATAATAAAAAAACT‑CCCGCCATTAAGTTGACTTTTAGCGcc                                                                                                           <  1:322727/142‑1 (MQ=255)
                                   tCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGccc                                                                                                          >  1:522938/1‑151 (MQ=255)
                                        tcttttttATTTCATCCAATCTTTTTGCTGCGAAAGTTTTTTCAAGAACAGTTAAGATGTTTTTTTTGTCGCTTAATGCTTTTAAAACATCCATGACCAAAAAAAAAATAAAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGCCCa                                                                                                         <  2:624585/144‑1 (MQ=255)
                                               tATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGa                                                                                               >  2:381392/1‑150 (MQ=255)
                                                   tCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgcc                                                                                          >  2:453874/1‑151 (MQ=255)
                                                        cAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGt                                                                                     >  1:404677/1‑151 (MQ=255)
                                                           tcttctTGCTGCGCAAGCGTTTCCAAGAACAGTTTAGATGATCTTTTTGTCGTTTAATGCCTGAAAAACATGCATGAGCCAAAAAAAAATAAAAAAAAC‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGccgccg                                                                                      <  2:595773/147‑1 (MQ=255)
                                                                ttGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                >  1:425533/1‑148 (MQ=255)
                                                                              ttttCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTa                                                                 >  1:576585/1‑149 (MQ=255)
                                                                                    aaaaCAGTTTAGATTATCTTTTTGTCGTTTAATGCCTGAAAAACATGCATGACCCAAAAAAAAAAAAAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAg                                                            <  2:469428/146‑1 (MQ=255)
                                                                                    aGAACGGGTTAGTTGTTTTTTTTGTCGCTTAATGCTTGAAAAACATGCATGAGCCAAAAAAAAAAAAAAAAAAC‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGa                                                           <  2:13757/149‑1 (MQ=255)
                                                                                                 tgatTTTTTTGTCTTTAAATGCCTGAAAAACATGAAGGACCCAAAAAAAAATAAAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTa                                             <  2:81612/150‑1 (MQ=255)
                                                                                                 tgatTTTTTTGTCGTTTAATGCTTGAAAAACATGCATGAGCAAAAAAAAAATAAAAAAAAC‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAtt                                               <  2:404714/148‑1 (MQ=255)
                                                                                                       ttttGTCGTTTAATGCCTAAAAAACAGGCATGAGCCAAAAAAAAATAAAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcat                                       <  1:268181/150‑1 (MQ=255)
                                                                                                             cctttatggcctaaaaaaCATGCATGACCAAAAAAATAAAAAAAAAAAC‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtg                                  <  1:282102/135‑1 (MQ=255)
                                                                                                               cTTAATGCCTGTAAAACATGCATGAGCCACAAAATACTATAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGa                               >  2:314099/1‑150 (MQ=255)
                                                                                                                ttAATGCTTGAAAAAAAGGCATGA‑CCAAAAAAAAAAAAAAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCAtt                                   <  2:576622/145‑1 (MQ=255)
                                                                                                                 tAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAtt                            >  2:611602/1‑151 (MQ=255)
                                                                                                                   aTGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTa                           >  2:56490/1‑150 (MQ=255)
                                                                                                                          aaaaaCAGGCATGACCAAAAAAATAAAAAAAAAATC‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                      <  1:169276/147‑1 (MQ=255)
                                                                                                                           aaaaCATGCATGAGCCACAAAATAATATAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaa                   >  2:23458/1‑150 (MQ=255)
                                                                                                                             aaCATGCATGAGCCACAAAATAATATAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                 >  2:622027/1‑150 (MQ=255)
                                                                                                                                    atgACCCAAAAAATAAAAAAAAAAAC‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg           <  1:487391/149‑1 (MQ=255)
                                                                                                                                      gAGCCACAAAATAATATAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagtg       >  1:316815/1‑151 (MQ=255)
                                                                                                                                      gAGCAAAAAAAAAATAAAAAAAAC‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg           <  1:473050/147‑1 (MQ=255)
                                                                                                                                         ccACAAAATAATATAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat    >  2:266573/1‑151 (MQ=255)
                                                                                                                                         ccACAAAATAATATAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat    >  2:424055/1‑151 (MQ=255)
                                                                                                                                           aaaaaaaaATAAAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGaa     <  1:585770/146‑1 (MQ=255)
                                                                                                                                           aaaaaaaaATAAAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat  <  2:65865/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                            
TTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAAT‑CCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAT  >  NC_000913/29346‑29635

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