Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1432,066,4450CG40.5% 52.4 / 10.5 42R7P (CGC→CCC) fatA_1Ferric‑anguibactin receptor FatA
Reads supporting (aligned to +/- strand):  ref base C (14/11);  new base G (0/17);  total (14/28)
Fisher's exact test for biased strand distribution p-value = 1.21e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTGCCTGGTAGCCATCATTGGCGTAGGTGTCGGAGACCTGCAATGGGGCCAGCATCACGCTGTCTTCGGCAAGCACAGAGGGGCACTGCAGGGCCAGGCCCAGGCTGAGCAGGCTCAGGGGCAGGCGGGGGTGATGCGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑TAAACCCACTCCTACACCGGGTGGTGATACACCAGATGGTGATCCCCCTGTA  >  ACB143/2066320‑2066564
                                                                                                                             |                                                                                                                        
gTGCGTGGTAGCCATCATTGGCGTAGGGGTCGTAGACCTGCATGGGGGCCAGCATCACGCGGTCTTCGGCAAGCACAGAGGGGCACTGCAGGGCCAGGCCCAGGCTGAGCAGGCGCAGGGGCAGGCGGGGGTGAt                                                                                                                 <  8:137960/135‑1 (MQ=255)
      ggTAGCCATCATTGGCGTAGGTGTCGGAGACCTGCAATGGGGCCAGCATCACGCTGTCTTCGGCAAGCACAGAGGGGCACTGCAGGGCCAGGCCCAGGCTGAGCAGGCTCAGGGGCAGGCGGGGGTGATGCggggg                                                                                                          >  3:42808/1‑136 (MQ=255)
                            gtCGGAGACCTGCAATGGGGCCAGCATCACGCTGTCTTCGGCAAGCACAGAGGGGCACTGCAGGGCCAGGCCCAGGCTGAGCAGGCTCAGGGGCAGGCGGGGGTGATGCGGGGGCATGCGGTGGGGCTCCTGGAt                                                                                     >  2:80125/1‑135 (MQ=255)
                            ggcGGACACCTGCAAGGGGGCCAGCACCCCGCTGTCTGCGGCACCCACCCTGGGCCACTGCCGGCCCCGGCCCTGGCTGAGCCGGCTCAGGGGCCGGCGGGGGTGAGGCGGGGGCATGCGGTGGGGCTCCTGGATg                                                                                    <  4:123061/134‑1 (MQ=255)
                             tCGGCGACCTGCAATGGGGCCAGCACCACGCTGTCTTCGGCACGCACAGAGGGCCACTGCAGGGCCAGGCCCAGGCTGAGCAGGCTCAGGGGCAGGCGGGGGGGATGCGGGGGCATGCGGGGGGGCTCCTGGATg                                                                                    <  7:54806/135‑1 (MQ=255)
                             gcgccgaCCTGCATGGGGGCCAGCAGCCCGCTGTCTGCGGCACCCACAGTGGGCCCCTCCTGGGCCAGGCCCTGGCTGAGCTGGGGCTGGGGCAGGGGGGGGTTTGGCGGGGGCATGCGGTGGGGCTCCTGGATg                                                                                    <  8:42623/130‑1 (MQ=255)
                                  gaCCTGCAATGGGGCCAGCATCACGCTGTCTTCGGCAAGCACAGAGGGGCACTGCAGGGCCAGGCCCAGGCTGAGCAGGCTCAGGGGCAGGCGGGGGTGATGCGGGGGCATGCGGTGGGGCTCCTGGATGcgcg                                                                                >  5:73855/1‑134 (MQ=255)
                                       gCAATGGGGCCAGCATCACGCTGTCTTCGGCAAGCACAGAGGGGCACTGCAGGGCCAGGCCCAGGCTGAGCAGGCTCAGGGGCAGGCGGGGGTGATGCGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCt                                                                          >  5:152525/1‑135 (MQ=255)
                                        cAATGGGGCCAGCATCACGCTGTCTTCGGCAAGCACAGAGGGGCACTGCAGGGCCAGGCCCAGGCTGAGCAGGCTCAGGGGCAGGCGGGGGTGATGCGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCt                                                                          >  5:100517/1‑134 (MQ=255)
                                                ccAGCATCCCGCTGTCTTCGGCAAGCACAGAGGGGCACTGCAGGGCCAGGCCCAGGCTGAGCAGGCTCAGGGGCAGGCGGGGGTGATGCGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGcc                                                                 >  1:328353/1‑135 (MQ=255)
                                                  aGCATCACGCTGTCTTCGGCAAGCACAGAGGGGCACTGCAGGGCCAGGCCCAGGCTGAGCAGGCTCAGGGGCAGGCGGGGGTGATGCGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCt                                                               >  1:203676/1‑135 (MQ=255)
                                                  aGCATCACGCTGTCTTCGGCAAGCACAGAGGGGCACTGCAGGGCCAGGCCCAGGCTGAGCAGGCTCAGGGGCAGGCGGGGGTGATGCGGGGGCATGAGGTGGGCCTCCTGGATGCGCGATACCTGTACTGGCCCt                                                               >  8:369423/1‑135 (MQ=255)
                                                   gCATCACGCTGTCTTCGGCACCCCCAGAGGGGCACTGCAGGGCCAGGCCCAGGCTGCCCCGGCTCTGGGGCCGGGGGGGGGGAGGCGGGGGCATGCGGGGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCtt                                                              <  8:213419/135‑1 (MQ=255)
                                                     aTCACGCTGTCTTCGGCAAGCACAGAGGGGCACTGCAGGGCCAGGCCCAGGCTGAGCAGGCTCAGGGGCAGGCGGGGGTGATGCGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTc                                                            >  2:142782/1‑135 (MQ=255)
                                                     aTCACGCTGTCTTCGGCAAGCACAGAGGGGCACTGCAGGGCCAGGCCCAGGCTGAGCAGGCTCAGGGGCAGGCGGGGGTGATGCGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTc                                                            >  6:234405/1‑135 (MQ=255)
                                                        aCGCTGTCTTCGGCAAGCACAGAGGGGCACTGCCGGGCCCGGCCCAGGCTGAGCCGGCTCTGGGGCAGGCGGGGGTGAGGCGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTcgcg                                                         <  7:54365/135‑1 (MQ=255)
                                                            tGTCTTCGGCAAGCACAGAGGGGCACTGCAGGGCCAGGCACAGGCTGAGCAGGCTCAGGGGCAGGCGGGGGTGATGCGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCgggg                                                      >  6:5747/1‑133 (MQ=255)
                                                            tGTCTTCGGCAACCACCGCGGGCCACTGCTGGGCCAGGCCCAGGCTGCCCCGGCTCCGGGGCAGGGGGGGGTGAGGCGGGGGCATGCGGGGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑Taa                                                   <  3:80196/136‑1 (MQ=255)
                                                            tGTCTGCGGCAACAAATGCGGGGCACGGCAGGGCACGGCCAAGGGTGAGCTGGCTCTGGGTCTGGGGGGGGTGATGGGGGGGCTTGGGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑Taa                                                   <  3:125701/136‑1 (MQ=255)
                                                            tGTCTGCGCCAAGCACTGCGGGGCACTGCAGGGCCAGGCCCAGGCTGAGCAGGCTCAGGGGCAGGGGGGGGTGATGCGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑Taa                                                   <  3:364313/136‑1 (MQ=255)
                                                             tgcTTCGGCACCCACCTCGGGCCACTGCTGGGCATGCCCCAGGCTGTGCAGGCACAGGGGCCGGGGGGGGTGATGGGGGGGCTTGCGGGGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑Taa                                                   <  8:279705/133‑1 (MQ=255)
                                                             ggctgcGGCACGCACGGGGGGGCACTGCTGGGCCTGGCCCTGGCTGCGCTGGGTCTGGGGCGGGGGGGGGTGTGGCGGGGGCTTGGGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑Taa                                                   <  5:31074/130‑1 (MQ=255)
                                                              tCTTCGGCACGCACTGCGGGGCACTGCAGGCCCCGGCCCAGGCTGACAAGGCGCAGGGCCTGGCGGGGGTCATGCGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑Taaa                                                  <  7:197427/135‑1 (MQ=255)
                                                               gtgcGGCAAGCACCGCGGGCCCCGGCTGGGACCGGCCCCGGCTGCCCTGGCTCGGGGGCTGGCGGGGGGGAGGCGGGGGCAGGGGGGGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑TAAAc                                                 <  2:285364/132‑1 (MQ=255)
                                                               attCGGTAAGCACTGCGGGCCACTGCGGGGCCAGGCCCCGGCTGAGCAGGCTACGGGGCAGGGGGGGGGGAGGCGGGGGCTGGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑Taa                                                   <  8:218939/132‑1 (MQ=255)
                                                                 tCGGCAAGCACAGAGGGGCACTGCAGGGCCAGGCCCAGGCTGAGCAGGCTCAGGGGCAGGCGGGGGTGATGCGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑TAAAccc                                               >  1:58135/1‑135 (MQ=255)
                                                                   ggCAAGCACAGAGGGGCACTGCAGGGCCAGGCCCAGGCTGAGCAGGCTCAGGGGCAGGCGGGGGTGATGCGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGGTTAAAcccc                                              >  1:62411/1‑134 (MQ=255)
                                                                            cgcgggcCACTCCAGGCCCCGGCCCTGGCTCCCCAGGCTCAGGGCCAGGGGGGGGTTATGCGGGGCCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑TAAACCCACTCCTACAcc                                    <  2:58135/132‑1 (MQ=255)
                                                                            atcggggCACTGCTGGGCCAGGCCCTGGCTGAGCAGGCTCTGGGGCAGGGGGGGGGGATGCGGGGGCTTGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑TAAACCCACTCCTACAcc                                    <  8:46068/132‑1 (MQ=255)
                                                                                   cACGGCCGGGCCTGGCCCAGGCGGAGCAGGCGCAGGGGCAGGCGGGGGTGATGCGGGGGCATGCGGGGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑TAAACCCACTCCTACACCGggtggt                             <  7:57145/135‑1 (MQ=255)
                                                                                   cACGGCCGGGCCCTGCCCAGGCTGAGCAGGCGCGGGGGCAGGCGGGGGGGATGCGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑TAAACCCACTCCTACACCGggtggt                             <  8:291704/135‑1 (MQ=255)
                                                                                    acggctgggCCTGGCCCTGGGTGAGCTGGCTCTGGGGCTGGGGGGGGTGAGGCGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑TAAACCCACTCCTACACCGGgtggtg                            <  6:152525/129‑1 (MQ=255)
                                                                                                 gccctggcTTAGCAGGCTCAGGGGCAGGGGGGGGTGATGCGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑TAAACCCACTCCTACACCGGGTGGTGATACACCAGATGGt              <  3:111980/136‑1 (MQ=255)
                                                                                                 gccctggcTGAGCAGGCTCCGGGGCCGGGGGGGGGGAGGGGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑TAAACCCACTCCTACACCGGGTGGTGATACACCAGATGGt              <  4:354764/136‑1 (MQ=255)
                                                                                                  ccctggcGGAGCTGGCTCCGGGGCTGGGGGGGGGGAGGGGGGGGCAGGCGGGGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑TAAACCCACTCCTACACCGGGTGGTGATACACCAGATGGt              <  5:314981/135‑1 (MQ=255)
                                                                                                  ccccggcTGAGCAGGCTCTGGGGCAGGGGGGGGTTATGCGGGGGCATGCGGGGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑TAAACCCACTCCTACACCGGGTGGTGATACACCAGATGGt              <  1:105690/135‑1 (MQ=255)
                                                                                                  ccccggcTGAGCAGGCGCTGGGGCAGGCGGGGGTGTTGGGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑TAAACCCACTCCTACACCGGGTGGTGATACACCAGATGGt              <  5:84387/135‑1 (MQ=255)
                                                                                                  cccaggcTGAGCAGGCTCTGGGGCTGGCGGGGGTGTTGCGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑TAAACCCACTCCTACACCGGGTGGTGATACACCAGATGGt              <  2:246426/135‑1 (MQ=255)
                                                                                                  cccaggcTGAGCACGCTCCGGGGCAGGGGGGGGTGATGCGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑TAAACCCACTCCTACACCGGGTGGTGATACACCAGATGGt              <  8:49524/135‑1 (MQ=255)
                                                                                                      ggcttcgCAGGCGCTGGGGCTGGGGGGGGGGAGGCGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑TAAACCCACTCCTACACCGGGTGGTGATACACCAGATGGTGATc          <  2:219216/129‑1 (MQ=255)
                                                                                                         tgcgccggCTCCGGGGCCGGCGGGGGTGAGGCGGGGGCTGGCGGGGGGGCGCCTGGATGCGCGATACCTGTACTGGCCCCTTCGCGGG‑TAAACCCACTCCTACACCGGGTGGTGATACACCAGATGGTGATcccc       <  2:281219/129‑1 (MQ=37)
                                                                                                              aGGCTCAGGGGCAGGCGGGGGTGATGCGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑TTAACCCACACATAACGCGGGTTGTGATTAAACACATGGTGGTgccccagga  >  2:301597/1‑131 (MQ=255)
                                                                                                                             |                                                                                                                        
GTGCCTGGTAGCCATCATTGGCGTAGGTGTCGGAGACCTGCAATGGGGCCAGCATCACGCTGTCTTCGGCAAGCACAGAGGGGCACTGCAGGGCCAGGCCCAGGCTGAGCAGGCTCAGGGGCAGGCGGGGGTGATGCGGGGGCATGCGGTGGGGCTCCTGGATGCGCGATACCTGTACTGGCCCTTTCGCGGG‑TAAACCCACTCCTACACCGGGTGGTGATACACCAGATGGTGATCCCCCTGTA  >  ACB143/2066320‑2066564

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