Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,574,0510GT42.9% 45.1 / 11.8 35V384V (GTC→GTA) pqqLputative periplasmic M16 family zinc metalloendopeptidase
Reads supporting (aligned to +/- strand):  ref base G (3/17);  new base T (15/0);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 2.93e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGCGC  >  NC_000913/1573904‑1574201
                                                                                                                                                   |                                                                                                                                                      
aCGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGAca                                                                                                                                                     <  2:657171/151‑1 (MQ=255)
    ttCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGAgatca                                                                                                                                                  >  1:198468/1‑150 (MQ=255)
      cGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACatcatc                                                                                                                                                <  1:648125/150‑1 (MQ=255)
                        ttCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGTTCTTTAGCAGaaaa                                                                                                                              >  2:342260/1‑150 (MQ=255)
                         tCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGaaaa                                                                                                                              <  2:246067/149‑1 (MQ=255)
                            ggCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCAt                                                                                                                         <  2:252686/151‑1 (MQ=255)
                            ggCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCAt                                                                                                                         <  1:522045/151‑1 (MQ=255)
                                    gAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATc                                                                                                                  <  2:118627/150‑1 (MQ=255)
                                       cGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGTTCTTCAGCAGAAAAAACATGCTGATCaa                                                                                                                >  2:242820/1‑149 (MQ=255)
                                              ttattTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGGTTTTCAGCAGAAAAAACAAGCTGCTCAAAGGGTGcc                                                                                                        >  2:213194/1‑150 (MQ=255)
                                                 tttAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTAAATCATCGAGTTCTTCAGCAGAAAAAACATGCTGATCAATGGTTGCCaa                                                                                                      >  2:475945/1‑149 (MQ=255)
                                                  ttAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGGGTTCTTTAGCAGAAAAAACATGCTGATCAATTGTTTCCaa                                                                                                      >  2:534908/1‑148 (MQ=255)
                                                    aaTGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTGACATCATCGAGTTTTTCCGCAGAAAAAACATGCTGATCAATGGGTGCCAActc                                                                                                   >  1:506820/1‑149 (MQ=255)
                                                           tGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGAATTTACATCCTCGAGGTTTTTAGCAGAAAAAACATGCTGATCCATGGTTGCCCAATCTGCCCtt                                                                                            >  2:561958/1‑149 (MQ=255)
                                                                   aGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGTTCTTTCGCAGAAAAAACATGCTGATCAATTGTTGCCAACTCTGCCCTTTATGCAtt                                                                                    >  2:527348/1‑149 (MQ=255)
                                                                   aGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCATCGAGTTTTTCAGCAGAAAAACCAAGCTGGTCAAAGGGTGCCAACCCTTCCCTTAATTCAtt                                                                                    >  1:273707/1‑149 (MQ=255)
                                                                       gACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAAACATGCTGATCAATGGGTGCCAACTCTGCCATTAATGCATTCgca                                                                                >  2:624995/1‑149 (MQ=255)
                                                                              cAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTg                                                                        <  1:132554/150‑1 (MQ=255)
                                                                                 cATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCata                                                                    <  2:280763/151‑1 (MQ=255)
                                                                                 cATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCATCGAGGTTTTTATCAGAAAAAACATGCTTGTCAAAGGGTGCCAACTCTGCCCTTAATGCATTTGCAGCATCCTGc                                                                       >  1:73266/1‑148 (MQ=255)
                                                                                    aCGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCAtatta                                                                 <  2:579459/151‑1 (MQ=255)
                                                                                                 tCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCa                                                     <  2:473814/150‑1 (MQ=255)
                                                                                                  cAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCCAGTTCTTTAGCCGAAAAACCATGCTGATCCAAGGGTGCCCACTCTTCCCTTAATGCATTCGCAGCCTCCCTCATATTATCGTCTCGTGCCtt                                                   >  2:577032/1‑151 (MQ=255)
                                                                                                   aGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCAttt                                                  <  2:330226/151‑1 (MQ=255)
                                                                                                        gCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACg                                             <  1:440734/151‑1 (MQ=255)
                                                                                                            aTCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGTTCTTTAGCAGAAAAAACAAGCTGGTCAATGGTTGCCCACCCTTCCATTAATTCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGaa                                          >  2:314177/1‑150 (MQ=255)
                                                                                                                     aTTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGTTCTTTCGCAGAAAAACCAAGCTGATCAATGGTTTCCAACCCTGCCATTAATGCATTTGCAGCATCCTGCATATTATTGTCTCCTGCATTTACACGGgaaaaaa                                     >  2:193748/1‑144 (MQ=255)
                                                                                                                       ttttCAGCCAGGTGAGGCGGGGAGAATTTACATCATCGAGTTCTTTCGCAGAAAAAACATGCTGATCAATGGTTTCCAACTCTGCCATTAATGCATTTGCAGCATCCTGCCTATTATCGTCTCGTGCATTTTCACGGaaaaaaaacg                                  >  1:460663/1‑142 (MQ=255)
                                                                                                                         ttCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGa                             <  1:252091/150‑1 (MQ=255)
                                                                                                                         ttCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGa                             <  1:547478/150‑1 (MQ=255)
                                                                                                                              ccAGGTGAGGCGGGGAGATTTTACATCCTCGGGGTCTTTAGCAGAAAAAACATGCTGATCAATGGGTGCCAACTCTTCCCTTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTTCACGGAAAAAAAACGACTGATaa                          >  1:227576/1‑148 (MQ=255)
                                                                                                                                 ggTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGt                    <  1:413284/151‑1 (MQ=255)
                                                                                                                                    gAGGCGGGGAGATTTTACATCATCGAGGTCTTTAGCAGAAAAACCATGCTGATCAATGGGTGCCCACTCTGCCATTTATGCATTTGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGGAAAACAGCGACTGGTAATCGGGTg                   >  1:216129/1‑149 (MQ=255)
                                                                                                                                                  tGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTgc    <  2:450830/150‑1 (MQ=255)
                                                                                                                                                  tGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTgc    <  1:134479/150‑1 (MQ=255)
                                                                                                                                                   gACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTgcgc  <  2:53382/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
ACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGCGC  >  NC_000913/1573904‑1574201

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