Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,569,9420G.69.0% 39.0 / 34.0 29coding (548/1536 nt)gadCL‑glutamate:4‑aminobutyrate antiporter
Reads supporting (aligned to +/- strand):  ref base G (4/5);  new base . (11/9);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCATTACCCGGAATGACCATCGCAATAGACAAACCACCAACAGAGCTTAAGCAGATTGCCGCCACCATCAGCAGTAACATAGCCAACGGATAGTCGCGCCCTGGGTTGCTCATTTCATTGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAAAATAAATGCAGGTAACAGGATACCGGCGAAGAAGCCAACTTTAGCAATTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCA‑ATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAAC  >  NC_000913/1569699‑1570185
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
tCATTACCCGGAATGACCATCGCAATAGACAAACCACCAACAGAGCTTAAGCAGATTGCCGCCACCATCAGCAGTAACATAGCCAACGGATAGTCGCGCCCTGGGTTGCTCATTTCATTGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGG‑CACCgg                                                                                                                                                                                                                                                >  1:147531/1‑249 (MQ=255)
     aCCCGGAATGACCATCGCAATAGACAAACCACCAACAGAGCGTAAGCAGATTGCCGCCACCATCAGCAGTAACATAGCCAACGGATAGTCGCGCCCTGGGTTGCTCATTTCATTGACGTTGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCTCCGGAGTGc                                                                                                                                                                                                                                           <  2:487535/250‑1 (MQ=255)
      cccGGAATGACCATCGCAATAGACAAACCACCAACAGAGCTTAAGCAGATTGCCGCCTCCATCAGCAGTAACATAGCCAACGGATAGTCGCGCCCTGGGTTGCTCATTTCATTGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGG‑CACCGGAGTGCAGa                                                                                                                                                                                                                                        <  1:43793/251‑1 (MQ=255)
         ggAATGACCATCGCAATAGACAAACCACCAACAGAGCTTAAGCAGATTGCCGCCACCATCAGCAGTAACATAGCCAACGGATAGTCGCGCCCTGGGTTGCTCATTTCATTGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGG‑CACCGGAGTGCAGATaa                                                                                                                                                                                                                                     <  1:357487/251‑1 (MQ=255)
               aCCAGCGCAATAGACAACCCACCAACAGAGCTTAAGCAGAGGTCCGCCACCATCAGCAGTAACATAGCCAACGGATAGTCGCGCCCTGGGTTGCTCATTTCATTGGCGTGGGTTGCGGATGCTTCTACTCCCATATAACTCAAAATGAACGCAACAACTACTACCAGGGAGCCCATTTTAGAGAAGTCAGTGAAGAAGGTCTTCCAAGCCATTTCGATATCAACGAGGGCACCGGACCGCAGATAAATAGc                                                                                                                                                                                                                                <  2:250713/251‑1 (MQ=255)
                    cGCAATAGACAAACCACCAACAGAGCTTAAGCAGATTGCCGCCACCATCAGCAGTAACATAGCCATCGGATAGTCGCGCCCTGGGTTGCTCATTTCATTGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGG‑CACCGGAGTGCAGATAAATAGCCGCTaa                                                                                                                                                                                                                          <  1:447953/251‑1 (MQ=255)
                    cGCAATAGACAAACCACCAACAGAGCTTAAGCAGATTGCCGCCACCATCAGCAGTAACATAGCCAACGGATAGTCGCGCCCTGGGTTGCTCATTTCATTGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCt                                                                                                                                                                                                                            >  1:242522/1‑250 (MQ=255)
                       aaTAGACAAACCACCAACAGAGCTTAAGCAGATTGCCGCCACCATCAGCAGTAACATAGCCAACGGATAGTCGCGCCCTGGGTTGCTCATTTCATTGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATg                                                                                                                                                                                                                        <  1:195114/251‑1 (MQ=255)
                                accacaAACAGAGCTTAAGCAGATTGCCGCCACCATCAGCAGTAACATAGCCAACGGATAGTCGCGCCCTGGGTTGCTCATTTCATTGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGG‑CACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAcaat                                                                                                                                                                                                              >  1:504102/1‑251 (MQ=255)
                                   accaACAGAGCTTAAGCAGATTGCCGCCACCATCAGCAGTAACATAGCCAACGGATAGTCGCGCCCTGGGTTTCTCATTTCATTGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCTCTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAaaataa                                                                                                                                                                                                            <  1:551601/251‑1 (MQ=255)
                                          gagCTTAAGCAGATTGCCGCCACCATCAGCAGTAACATAGCCAACGGATAGTCGCGCCCTGGGTTGCTCATTTCATTGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGG‑CACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAAAATAAATGCAGGt                                                                                                                                                                                                    <  2:501887/251‑1 (MQ=255)
                                                   cAGATTGCCGCCACCATCAGCAGTAACATAGCCAACGGATAGTCGCGCCCTGGGTTGCTCATTTCATTGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGG‑CACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAAAATAAATGCAGGTAACAGGAt                                                                                                                                                                                            >  1:335360/1‑250 (MQ=255)
                                                            gccACCATCAGCAGTAACCTAGCCAACGGATAGTCGCGCCCTGGGTTGCTCCTTTCATTGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGG‑CCCCGGGGTGCAGATAACTACCCGCTCATGCGATCAAAATAAATGCAGGAAACAGGATACCGGcaa                                                                                                                                                                                    >  2:403361/1‑247 (MQ=255)
                                                                     agcagTAACATAGCCAACGGATCGTCGCGCCCTGCGTTGCTCATTTCATTGACGAGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCCACAAATACTAACAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGG‑CACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAAAATAAATGCAGGTAACAGGATACCGGCGAAGCAGCCAACt                                                                                                                                                                         >  2:285916/1‑251 (MQ=255)
                                                                                         aTAGTCGCGCCCTGGGTTGCTCATTTCATTGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGG‑CACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAAAATAAATGCAGGTAACAGGATACCGGCGAAGAAGCCAACTTTAGCAATTCGCGCCgtgt                                                                                                                                                      >  1:242846/1‑250 (MQ=255)
                                                                                                 gcCCTGGGTTGCTCATTTCATTGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGG‑CACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAAAATAAATGCAGGAAACAGGATACCGGCGAAGAAGCAAACTTTAGCAATTCGCGCCGTGTATTTCGTGc                                                                                                                                             >  2:529989/1‑251 (MQ=255)
                                                                                                        gTTGCTCATTTCATTGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGACCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGG‑CACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAAAATAAATGCAGGTAACAGGATACCGGCGAAGAAGCCAACTTTAGCAATTCGCGCCGTGTATTTCGTGCCACCaaa                                                                                                                                      <  2:335360/251‑1 (MQ=255)
                                                                                                                       gACGGGGGTTGCGGATGCTTCTACGCCCAAATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTGTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGG‑CACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAAAATAAATGCAGGTAACAGGATACCGGCGAAGAAGCCAACTTTAGCAATTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCa                                                                                                                        <  2:989/250‑1 (MQ=255)
                                                                                                                                               gCCCATATAACTCAAAATGAAGGCAACACATACTACCAGGGTGCACACTTTAGAGACGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGG‑CACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAAAATAAATGCAGGTAACAGGATACCGGCGAAGAAGCCAACTTTAGCAATTCGCGCCGTGTATTTCGTGCCACCAAACTGCGATAATGCCAGCGCCCAAAGAATGATGAGtgctgc                                                                                               >  2:139375/1‑251 (MQ=255)
                                                                                                                                                   aaataACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAAAATAAATGCATGTAACAGGATACCGGCGAAGAAGCCAACTTTAGCAATTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCA‑At                                                                                           <  2:538637/249‑1 (MQ=255)
                                                                                                                                                           cAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAAGTCAGGGAAGCAGGTCTTCGAATCCATTTCGATAGCAACGGGG‑CACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAAAATAAATGCAGGTAACAGGATACCGGCGAAGAAGCCACCTTTAGCAATTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCA‑ATAGTTTTGGt                                                                                  >  2:414089/1‑251 (MQ=255)
                                                                                                                                                                        acaaATACTACCAGGGTGCCCACTTTAGAGAAGTCAGGGAAGCAGGTCTTCGAATCCATTTCGATCGCAACGGGG‑CACCGGAGTGCAGATAAATAGCCGCTAATGCGATCCAAATAAATGCAGGTAACAGGATACCGGCGAAGAAGCCAACTTTAGCAATTCGCGCCGTGTATTTCGTGCCACCAACCTGCGTTAATGCCATCGCCCAAAGAATGATTAGTGCTGCACAT‑GCTTTGGTAATGGGATCttct                                                                     >  2:3691/1‑250 (MQ=255)
                                                                                                                                                                           aaTACTACCAGGGTGACCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAAAATAAATGCAGGTAACAGGATACCGGCGAAGAAGCCAACTTTAGCAATTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCGGCA‑ATAGTTTTGGTAATGGGGTCttcatt                                                                   >  1:482013/1‑251 (MQ=255)
                                                                                                                                                                                  ccAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGG‑CACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAAAATAAATGCAGGTAACAGGATACCGGCGAAGAAGCCAACTTTAGCAATTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCA‑ATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTg                                                           <  1:414089/251‑1 (MQ=255)
                                                                                                                                                                                   cAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAAAATAAATGCAGGTAACAGGATACCGGCGAAGAAGCCAACTTTAGCAATTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCA‑ATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTg                                                           >  1:256914/1‑251 (MQ=255)
                                                                                                                                                                                                          cAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGG‑CACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAAAATAAATGCAGGTAACAGGATACCGGCGAAGAAGCCAACTTTAGCAATTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCA‑ATAGTTTGGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGt                                    <  2:266580/250‑1 (MQ=255)
                                                                                                                                                                                                              gaagaaGGTCTTCGAATCCATTTCGATAGCAACGGGG‑CACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAAAATAAATGCAGGTAACAGGATACCGGCGAAGAAGCCAACTTTAGCAATTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCA‑ATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGcccc                               >  2:552232/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aTAGCAACGGGG‑CACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAAAATAAATGCAGGTAACAGGATACCGGCGAAGAAGCCAACTTTAGCAATTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCA‑AGAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATa      <  1:285916/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aaCGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAAAATAAATGCAGGTAACAGGATACCGGCGAAGAAGCCAACTTTAGCAATTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCA‑ATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAAc  >  1:527291/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
TCATTACCCGGAATGACCATCGCAATAGACAAACCACCAACAGAGCTTAAGCAGATTGCCGCCACCATCAGCAGTAACATAGCCAACGGATAGTCGCGCCCTGGGTTGCTCATTTCATTGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAAAATAAATGCAGGTAACAGGATACCGGCGAAGAAGCCAACTTTAGCAATTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCA‑ATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAAC  >  NC_000913/1569699‑1570185

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