Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009974542,1160AG69.8% 39.7 / 26.0 43intergenic (‑377/+152)argD/RPPX_02300acetylornithine aminotransferase/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/13);  new base G (6/24);  total (6/37)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGGCGATGTTGATGTGTTGGCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTCCCTCCTGTGCCTGGCGCGCGCGACCGCCTTGTTTGTACGACGCGATATCGAGGCATGCGCCAAGGCGTTCGCGCGCAAATGCCACAGGCATAACTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATTTGCGCCCCAACATCAACCTCAAGGCATACACCCCG  >  CP009974/541971‑542252
                                                                                                                                                 |                                                                                                                                        
cGGCGATGTTGATGTGTTGGCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTgcgc                                                                                                                                     <  2:1366399/151‑1 (MQ=17)
 ggCGATGTTGATGTGTTGGCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGcc                                                                                                                                    <  1:2797393/151‑1 (MQ=14)
     gtGTTGATGTGTTGGCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAgg                                                                                                                                <  2:2715209/150‑1 (MQ=14)
     atgtggatgtggTGGCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAgg                                                                                                                                <  2:4529291/151‑1 (MQ=11)
     aTGTTGATGTGTTGGCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAgg                                                                                                                                <  1:4268932/151‑1 (MQ=14)
     aTGTTGATGTGTTGGCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAgg                                                                                                                                <  2:2517038/151‑1 (MQ=14)
      tGTTGATGTGTTGGCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCCGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAGGc                                                                                                                               <  2:397704/151‑1 (MQ=14)
          gATGTGTTGGCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAGGCGTTc                                                                                                                           <  2:4055769/151‑1 (MQ=14)
          gATGTGTTGGCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAGGCGTTc                                                                                                                           <  2:1370783/151‑1 (MQ=14)
          gATGTGTGGGCGGGGAGCTCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAACCCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAGGCGTTc                                                                                                                           <  2:4785766/151‑1 (MQ=11)
           aTGTGTTGGCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAGGCGTTcg                                                                                                                          <  2:2559733/151‑1 (MQ=14)
            tgtgTTGGCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAGGCGTTcgc                                                                                                                         <  1:1727486/151‑1 (MQ=14)
             gtgtTGGCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGcc                                                                                                                                    >  1:1409302/1‑139 (MQ=14)
             gtgtTGGCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGcc                                                                                                                                    <  2:1409302/139‑1 (MQ=14)
             gtgtTGGCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAGGCGTTcgcg                                                                                                                        <  2:806478/151‑1 (MQ=11)
              tgtTGGCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAGGCGTTcgcgc                                                                                                                       <  2:4721496/151‑1 (MQ=11)
              tgtTGGCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAGGCGTTcgcgc                                                                                                                       <  1:4468498/151‑1 (MQ=11)
              tgtTGGCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTGCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAGGCGTTcgcgc                                                                                                                       <  1:4473157/151‑1 (MQ=11)
               gtTGGCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAGGCGTTcgcgcg                                                                                                                      <  1:434010/151‑1 (MQ=11)
                   gCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAGGCGTTCGCGCGCaaa                                                                                                                  >  1:2224529/1‑151 (MQ=11)
                   gCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGCGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGTGATATCGAGGTGTGCGCCAAGGCGTTCGCGCGCaaa                                                                                                                  <  2:4253501/151‑1 (MQ=11)
                    cGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAGGCGTTCGCGCGCAAAt                                                                                                                 <  1:5026609/151‑1 (MQ=11)
                    cGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAGGCGTTCGCGCGCAAAt                                                                                                                 <  2:463876/151‑1 (MQ=11)
                    cGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAGGCGTTCGCGCGCAAAt                                                                                                                 <  1:4556285/151‑1 (MQ=11)
                    cGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAGGCGTTCGCGCGCAAAt                                                                                                                 <  1:2045426/151‑1 (MQ=11)
                     ggggAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGCTTGTTTCGGGCCAGTCCCTCCTGCGCCCGGCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAGGCGTTCGCGCGCAAATg                                                                                                                <  1:939785/151‑1 (MQ=11)
                          gaccGAGCGCCGCCCGCGCGGCGCTGCGCGGGTCAACCCGCCCCTCGGTTTGTTTCGGGCCAGTCCCTCCTGCACCAGTCGCGCGCGACCGCCTTGTTTGTACGCCGCGATATCGAGACATGCGCCAAGGCGGTCGCGCGCAAATCCcaca                                                                                                           <  2:2112582/148‑1 (MQ=1)
                             cgggcgcCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTCCATCCTGTGCCTGGCGCGCGCGACCGCCTTGTTTGTACGACGCGATATCGAGGCATGCACCAAGGCGTTCGCGCGCAAATCCCACAGGc                                                                                                        <  1:1069479/148‑1 (MQ=0)
                                gccccgcccgcGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTCCATCCTGTGCCTGGCGCGCGCGACCGCCTTGTTTGTACGACGCGATATCGAGGCATGCACCAAGGCGTTCGCGCGCAAATCCCACAGGCATa                                                                                                     <  2:1826651/148‑1 (MQ=0)
                                        cgagcggcgCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTCCCTCCTGCACCAGTCGCGCGCGACCGCCTTGTTTGTACGACGCGATATCGAGACATGCGCCAAGGCGGTCGCGCGCAAATCCCACAGGAATAATTGGccc                                                                                             <  2:1608207/148‑1 (MQ=1)
                                        cgagcggcgCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTCCCTCCTGCACCAGTCGCGCGCGACCGCCTTGTTTGTACGACGCGATATCGAGACATGCGCCAAGGCGGTCGCGCGCAAATCCCACAGGAATAATTGGccc                                                                                             <  1:4900378/148‑1 (MQ=1)
                                           gccgcgcTTGGCGGGTAGTCGCGCGCCTTGGTTTGTTCTGGGCCAGTCCCTCCGGCACCAGTCGCGCGGGACCGCCTTGTTTGTACGACGCGATATCGAGACATGCGCCAAGGCGGTCGCGCGCAAATCCCACAGGAATAATTGGCCCGaa                                                                                          <  2:1051559/148‑1 (MQ=0)
                                           gccgcgcTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTCCCTCCTGCACCAGTCGCGCGCGACCGCCTTGTTTGTACGACGCGATATCGAGACATGCGCCAAGGCGGTCGCGCGCAAATCCCACAGGAATAATTGGCCCGaa                                                                                          <  1:1121807/148‑1 (MQ=1)
                                                             acgcgatCCTACGTTTGTTTCGGGCCAGCGCCTCCTGTGCCTGGCGCGCGCGACCGCCTTGTTGGTACGACGCGATATCGAGACAGGCACCAAGGCGTTCGCGCGCGAATCCCACAGGCATAATTGGCCCGAAACAAACGTAGGAGTGGGt                                                                        <  2:4482635/145‑1 (MQ=1)
                                                                 gccccTACGTTTGTTGCGGGCCGGGCCCTCCTGCACCAGTCGCGCGGGACCGCCTGGTTTGTACGACGCAGTAGCGAGACATGCGCCAAGGCGGGCGCGCGCAAATCCCACAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTAc                                                                    <  2:3361807/148‑1 (MQ=1)
                                                                         gtgtgtttGGGGCCAGTCCATCCTGGGCCTGCCGCGCGCGACCGCCTTGTTTGTACGACGCGAGATCGAGGCATGCACCAAGGCGTTGGCGCGCAAACCCCACAGGCATAATTGGCCTGAAACAAACGTAGGAGCGGGTTTACCCGCGAAg                                                            <  2:4937566/148‑1 (MQ=1)
                                                                           ttttttCGGGCCAGTCCCTCCTGTGCCTGGCGCTCGCGACCGCCTTGTTGGTACGACGCGATATCGAGACAGGCACCAAGGCGTTCGCGCGCGAATCCCACAGGCATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAgcg                                                          <  1:3095526/148‑1 (MQ=1)
                                                                           ttggtacGGTCCAGTCCCCCCTGTGCCTGGCGCGCGCGCCCGCCTTGGTGGTACAACGCGAAAGCGAGACAGGCACCAGGGCGTTCGCGCGCGAATCCCACAGGCATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAgcg                                                          <  2:1843621/145‑1 (MQ=1)
                                                                             gtgtCGGGCCAGTCCCTCCTGTGCCTGGCGCGCGCGACCGCCTTGTTGGTACGACGCGATATCGAGACAGGCACCAAGGCGTTCGCGCGCGAATCCCACAGGCATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGcgcc                                                        <  1:2331502/148‑1 (MQ=1)
                                                                                                tgcgcccggCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAGGCGTTCGCGCGCAAATGCCACAGGCATAACTGGCCCGAAACAAACGTAGGAGCAGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtt                                     >  2:3452595/8‑151 (MQ=1)
                                                                                                tgcgcccggCGCGCGCGACCGCCTTGTTTGTTCGACGGGATATCGAGGTGTGCGCCAAGGCGTTCGCGCGCAAATGCCACAGGCATAACTGGCCCGAAACAAACGTAGGAGCAGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtt                                     >  2:3634354/8‑151 (MQ=1)
                                                                                                                             gttcgacgGGATATCGAGGTGTGCGCCAAGGCGTTCGCGCGCAAATGCCACAGGCATAACTGGCCCGAAACAAACGTAGGAGCAGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATTTGCGCCCCAACATCAACCTCAAGGCATac        >  2:2809562/4‑151 (MQ=1)
                                                                                                                                   cgggATATCGAGGTGTGCGCCAAGGCGTTCGCGCGCAAATGCCACAGGCATAACTGGCCCGAAACAAACGTAGGAGCAGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATTTGCGCCCCAACATCAACCTCAAGGCATACACCCCg  >  1:4989991/4‑151 (MQ=0)
                                                                                                                                                 |                                                                                                                                        
CGGCGATGTTGATGTGTTGGCGGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTCCCTCCTGTGCCTGGCGCGCGCGACCGCCTTGTTTGTACGACGCGATATCGAGGCATGCGCCAAGGCGTTCGCGCGCAAATGCCACAGGCATAACTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATTTGCGCCCCAACATCAACCTCAAGGCATACACCCCG  >  CP009974/541971‑542252

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑