Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,337,7230AC52.5% 46.9 / 65.3 40G566G (GGT→GGG) gyrADNA gyrase (type II topoisomerase), subunit A
Reads supporting (aligned to +/- strand):  ref base A (18/1);  new base C (0/21);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 1.68e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.67e-08

GAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGAC  >  NC_000913/2337587‑2337871
                                                                                                                                        |                                                                                                                                                    
gAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCacg                                                                                                                                        >  1:476538/1‑151 (MQ=255)
        tCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTgcg                                                                                                                                >  1:78042/1‑151 (MQ=255)
              aaCTGATAACCTTTAATCGAATAGACGCGCCCACGGCTGGGGAAGCACAGAATAGGGTCGGGAGTGTTCGCCCCCAGCATTGGGTCGATAAAGTCTTTTTCTTTAATACGTGCGGAAGTTTTCCCTTGCCCGCCACGACGCTGCGCTTCGt                                                                                                                          <  1:713921/151‑1 (MQ=255)
                      aaCTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCCCCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGaa                                                                                                                  <  2:464196/151‑1 (MQ=255)
                                gAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGAt                                                                                                        >  2:179306/1‑151 (MQ=255)
                                gAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGAt                                                                                                        >  1:218402/1‑151 (MQ=255)
                                  aTAGACGCGACCACGGTGGGGGAAGCAGAGAATATGGTCGTGAGTGTTCGCCCCCAGCAGTCGGTCGATAAATTCTTCTTTTTTAATACGTGGGGAAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATAc                                                                                                      <  1:59818/151‑1 (MQ=255)
                                         cgACCACGCCTGGAGAAGCACAAAATAGGGTCGTGATTGTCCGCCACAAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATGCGTGCGGCAGATTTCCTTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGt                                                                                               <  1:452740/151‑1 (MQ=255)
                                             cACGGCTGGAGAAGCACAGAATATGGTCGTGAGTTTTCCCCCCCAGCTGTCGTTCGATAAACTCTTCTTCTTTAATACGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGct                                                                                           <  1:104859/151‑2 (MQ=255)
                                             cACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGcc                                                                                           <  2:144241/151‑1 (MQ=255)
                                              ccGGCTGGAAAAGCACAGAATATGGTCTTGAGTGTCCCCCACCAGCATTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGCCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGccc                                                                                          <  2:111490/150‑1 (MQ=255)
                                              aCGGCGGGAGAAGCACAGAATAGGGTCGTGAGTGTTCGCCCCCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGccc                                                                                          <  2:510961/151‑1 (MQ=255)
                                                ggCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTTTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGcc                                                                                           >  2:253324/1‑148 (MQ=255)
                                                  cTGGAGAAGCACAGAATTGGGTCGTGAGTGTCCGCCCCCAGCAGTCGGTCGATAAAGTCTTCTTTTTTAATCCGTGCGGCAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTgg                                                                                       <  1:566722/150‑1 (MQ=255)
                                                   tGGAGAGGCAAAGATTTGGTTTGTGAGTGTTCCCCCCCAGGAGTGGCTCGAAAAAGTCTTCTGTTTTAATGCGTGGGGAAGATTGCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTg                                                                                     <  1:532466/151‑1 (MQ=255)
                                                   tGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGCCGGTCGATAAAGTTTTCTTTTTTAATCCGGGCGGCAGATTTCCCTTTCCCGCCCCGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTg                                                                                     <  2:571721/151‑1 (MQ=255)
                                                     gagaAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgag                                                                                   >  2:407061/1‑151 (MQ=255)
                                                        aaGCACAGAATATGGTCGTGATTGTCCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATCCGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgagaga                                                                                <  1:158689/151‑1 (MQ=255)
                                                          gCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGagagag                                                                               >  1:274886/1‑150 (MQ=255)
                                                          gCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGagagag                                                                               >  2:474134/1‑150 (MQ=255)
                                                               gAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCa                                                                          >  1:736119/1‑150 (MQ=255)
                                                               gAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCa                                                                          >  1:665713/1‑150 (MQ=255)
                                                               gAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCa                                                                          >  1:575757/1‑150 (MQ=255)
                                                               gAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCa                                                                          >  1:374242/1‑150 (MQ=255)
                                                                 atatGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGa                                                                       >  1:283691/1‑151 (MQ=255)
                                                                     ggTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACca                                                                    >  1:451046/1‑150 (MQ=255)
                                                                          ggAGTGTCCCCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCt                                                               <  1:277215/149‑1 (MQ=255)
                                                                                 tCGCCACCAGAAGTCGGTCGATAAATTCTTCTTCTTAAATACGTGGGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGt                                                       <  2:424009/151‑1 (MQ=255)
                                                                                     ccccaGCAGTCGGTCGATAAAGTCTTCTTTTTTAATCCGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATc                                                   <  2:538655/149‑1 (MQ=255)
                                                                                     cacaaGCAGTGGGTCGATAAATCCTTCTTCTTTAATACGGGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATc                                                   <  1:211081/151‑1 (MQ=255)
                                                                                       ccaGCAGTCGGTCGATAAATTCTTCTTCTTTAATCCGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAg                                                 <  2:82133/151‑1 (MQ=255)
                                                                                           cattcggtcgAAAAAGTCTTCTTTTTTAATACGTGCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATc                                              <  1:77956/147‑1 (MQ=255)
                                                                                               cggtcgATAAAGCTTTCTTCTTAAATGCGTGCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCa                                         <  1:302170/151‑1 (MQ=255)
                                                                                                  tcgAAAAAGTCTTTTTCTTTAATACGTGCGGAAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGt                                      <  2:314312/151‑1 (MQ=255)
                                                                                                   cgATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGtt                                     >  2:364521/1‑151 (MQ=255)
                                                                                                         aGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTTAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGt                                >  2:244989/1‑150 (MQ=255)
                                                                                                              ttttcttAAATCCGGGCGGCATTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTgcgc                          <  2:546295/149‑1 (MQ=255)
                                                                                                                 tcttTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTgcg                           >  2:585209/1‑147 (MQ=255)
                                                                                                                                     tttACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCCGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTacga   >  1:354800/1‑151 (MQ=255)
                                                                                                                                      ttACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTacgac  <  2:562632/151‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                    
GAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGAC  >  NC_000913/2337587‑2337871

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: ‑