Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273822,9550TG46.9% 75.7 / 14.2 65T495P (ACC→CCC) ybhFATP‑binding cassette domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (3/31);  new base G (30/0);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 3.39e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAA  >  NZ_CP009273/822871‑823044
                                                                                     |                                                                                          
cTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAggg                                                                                       >  1:2370953/1‑90 (MQ=255)
  tCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGgggggggggg                                                                                     >  1:3181916/1‑90 (MQ=255)
  tCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAggggg                                                                                     >  2:685784/1‑90 (MQ=255)
    aTCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGgggggggggggg                                                                                   >  2:2518789/1‑88 (MQ=255)
         tAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGc                                                                              >  1:1720436/1‑90 (MQ=255)
              tgggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGaa                                                                         >  2:767216/1‑90 (MQ=255)
               gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAg                                                                        >  1:2454490/1‑90 (MQ=255)
               gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAg                                                                        >  2:2517594/1‑90 (MQ=255)
               gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAg                                                                        >  2:3344535/1‑90 (MQ=255)
                  tggtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCg                                                                    <  1:537973/90‑1 (MQ=255)
                      gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGGCGGGtc                                                                >  2:626200/1‑90 (MQ=255)
                      gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGGAGGCgggtt                                                                >  1:2518231/1‑87 (MQ=255)
                       aCCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCCACGGCCGAAGGCGGGtcg                                                               >  1:361552/1‑90 (MQ=255)
                          aTCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtcgtcg                                                            <  1:2483848/90‑1 (MQ=255)
                                gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGAAGGCGGTCCGTCGAGAAAc                                                      >  2:1922295/1‑90 (MQ=255)
                                gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAGGCCCGAAGTCGGTTCGTCGAGaaaa                                                      >  1:1629216/1‑89 (MQ=255)
                                gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTTGAGaaaa                                                      >  2:736836/1‑89 (MQ=255)
                                 tGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACa                                                     <  1:1920557/90‑1 (MQ=255)
                                   aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGa                                                   <  2:2370953/90‑1 (MQ=255)
                                    cGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCTGGAAGTCGGTTCGTCGAGAAACAGaa                                                  <  2:2940759/90‑1 (MQ=255)
                                       cTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGt                                               <  2:3311080/90‑1 (MQ=255)
                                        tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGCAGGCGGGTCGTCGAGGAAAAAAAAGTc                                              >  2:452169/1‑90 (MQ=255)
                                         ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                             <  1:632251/90‑1 (MQ=255)
                                         ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                             <  1:93471/90‑1 (MQ=255)
                                         ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                             <  2:1846310/90‑1 (MQ=255)
                                         ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                             <  2:1827063/90‑1 (MQ=255)
                                             ctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGtt                                         <  2:2554457/90‑1 (MQ=255)
                                                ccATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCAt                                      <  2:1602050/90‑1 (MQ=255)
                                                ccATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACCCCGGAAGGCGGGTCGGGGGGAAAAAAAAAGTCCGGGTCAt                                      >  1:3029294/1‑90 (MQ=255)
                                                 cATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCGACCCCCGAAGGCGGTTCGTCGGGGAAAAAAAAGTCCGGGTTATg                                     >  2:2545143/1‑90 (MQ=255)
                                                  aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc                                    <  2:2612405/90‑1 (MQ=255)
                                                  aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc                                    <  1:557979/90‑1 (MQ=255)
                                                  aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc                                    <  2:3217263/90‑1 (MQ=255)
                                                  aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc                                    <  2:1629216/90‑1 (MQ=255)
                                                  aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGGACGCCGGGAGGCGGGTCGTCGGGAAAAAAAAAGTCCGGGTCATGc                                    >  2:1627951/1‑90 (MQ=255)
                                                  aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGAAGGCGGGTCGTCGGGAAAAAGAAAGTGCGGGTCGTGc                                    >  2:3000799/1‑90 (MQ=255)
                                                  aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGTTCGTCGAGAACCAGAATGTCCGGGTGATGc                                    >  2:2921415/1‑90 (MQ=255)
                                                  aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGAAGGCGGGTCGTCGGGAAAAAAAAAGTCCCGGTCATgg                                    >  1:436955/1‑89 (MQ=255)
                                                  aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCCGGAGTCGGTTTGTCGGGAAAAAAAAAGTCCGGGTTATGc                                    >  1:655271/1‑90 (MQ=255)
                                                    gCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGAAGGCGGGTCGGCGGGGAAAAAAAAGGCCGGGTCGTGCAt                                  >  2:2480/1‑90 (MQ=255)
                                                     cTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGTCGGGGCGGCGAGAAAAAGAAAGTGCGGGTCCTGCATc                                 >  2:2965580/1‑90 (MQ=255)
                                                     cTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGGCGGGTCGTCGGGAAAAAAAAAGTCCGGGTCATGCATc                                 >  1:3296725/1‑90 (MQ=255)
                                                     cTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAAGCGGGTCGTCGGGGAAAAAAAAGGCCGGGTCATGCaac                                 >  2:1681616/1‑88 (MQ=255)
                                                       gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGGAACCCCCGGAGTCGGGTTGTCGAGGAAAAAAAAGTCCGGGTTCTGCATCa                                >  1:2302954/1‑90 (MQ=255)
                                                       gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGGACGCCGGGAGGGGGGGGGTCGAGAAAAAAAAAGGCCGGGTTCTGCATCAg                               >  2:3193129/1‑90 (MQ=255)
                                                            tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                          <  1:489251/90‑1 (MQ=255)
                                                            tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                          <  1:3223458/90‑1 (MQ=255)
                                                            tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                          <  2:3013868/90‑1 (MQ=255)
                                                            tgtgCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGGCGGGTCGTCGAGAAAAAGAAAGTCCGGGTCATGCATCAGCGaaa                          >  2:2288928/1‑89 (MQ=255)
                                                             gtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACa                         <  2:1832840/90‑1 (MQ=255)
                                                             gtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACa                         <  2:552428/90‑1 (MQ=255)
                                                             gtgCAGCCAAAATTCACGGCGGGGGAGGGGGTCAACGCCGGGAGCCGGTGCGGCGAGAAACAGAACGACCAGGTCATGCCTCAGCGaaaa                         >  2:3049508/1‑88 (MQ=255)
                                                               gCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAACAGAAAGTCCGGGTCATGCATCAACGGAAAgg                       >  1:2583650/1‑90 (MQ=255)
                                                               gCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGTCGGTTCGTCGGGAAAAAGAAAGTCCGGTTCATGCACCAGCGGACAgg                       >  1:1489464/1‑90 (MQ=255)
                                                                 aGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgcc                     <  1:718039/90‑1 (MQ=255)
                                                                 aGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGGCGGGTGGTCCAGAAACAGAAAGTCCGGGTCGTGCAGCAGCGAGCAGgcc                     >  2:3015878/1‑90 (MQ=255)
                                                                 aGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGAGAAAAAGAAAGTCCGGGTTATTCATCAGCGAACAAgcc                     >  1:2512919/1‑90 (MQ=255)
                                                                 aGCCAAAATTCACGGCGGGCGAAGGGGAAAACGCCGCAAGTACGGTCCTTCGAAAATCAAAAACTCGGTTCTTGCTTCAGGGAACaaggc                     >  1:1094308/1‑86 (MQ=255)
                                                                  gCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgcca                    <  2:3181916/90‑1 (MQ=255)
                                                                        aTTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagcgcca              <  1:626200/90‑1 (MQ=255)
                                                                             cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAACCAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc         <  1:2921415/90‑1 (MQ=255)
                                                                             cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc         <  1:1032951/90‑1 (MQ=255)
                                                                             cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc         <  2:1489464/90‑1 (MQ=255)
                                                                             cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc         <  1:3061754/90‑1 (MQ=255)
                                                                             cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc         <  1:1627951/90‑1 (MQ=255)
                                                                             cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc         <  1:1576031/90‑1 (MQ=255)
                                                                              ggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCt        <  2:3265153/90‑1 (MQ=255)
                                                                                     tGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaa  <  2:1199960/90‑1 (MQ=255)
                                                                                     |                                                                                          
CTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAA  >  NZ_CP009273/822871‑823044

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: