Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,047,9800CA26.8% 111.8 / 11.8 56A37S (GCC→TCC) PP_5298glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base C (20/21);  new base A (0/15);  total (20/36)
Fisher's exact test for biased strand distribution p-value = 4.01e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTGGGGGCGCCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCAGGTACC  >  NC_002947/6047847‑6048112
                                                                                                                                     |                                                                                                                                    
cTGGGGGCGCCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGTGCCGGCAAGATCAGTGGCAAGCCGGCaa                                                                                                                                    >  2:808534/1‑136 (MQ=255)
      gcgcCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAgccgg                                                                                                                              >  1:240106/1‑136 (MQ=255)
       cgcCATTGTAATGATGCGGCTCGATATTTGACGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAgccggc                                                                                                                             <  1:762705/136‑1 (MQ=255)
                       cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCATCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTg                                                                                                             >  1:781580/1‑136 (MQ=255)
                       cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTg                                                                                                             >  2:131854/1‑136 (MQ=255)
                       cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTg                                                                                                             >  1:565149/1‑136 (MQ=255)
                       cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTg                                                                                                             >  1:421645/1‑136 (MQ=255)
                       cGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAACCGGCCGCCTCGACATACTTg                                                                                                             >  1:461231/1‑136 (MQ=255)
                            cGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGcc                                                                                                        >  1:12687/1‑136 (MQ=255)
                            cGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGcc                                                                                                        >  1:749111/1‑136 (MQ=255)
                            cGATATTTGAAG‑‑GAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCgc                                                                                                      >  1:288646/1‑135 (MQ=255)
                                 tttGACGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCg                                                                                                   <  2:417070/136‑1 (MQ=255)
                                    taaGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCgtgt                                                                                                <  1:336065/135‑1 (MQ=255)
                                      cggcttgCCGTTAAAACAAATGCGGTGCTGGCGAGCCAGCAGCGTTTGGGTGTCGCGGCCGCCATCCGGGGCCGGCAACATCTGTGGCACGCGGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGAACGTGTGc                                                                                              <  1:760306/130‑1 (MQ=255)
                                         cGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAGGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCgag                                                                                           <  2:785403/136‑1 (MQ=255)
                                          gacccGTTAAAAACAATGCCGTGCAGGCGAGCCAGCTGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCTGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCgagg                                                                                          <  2:512835/133‑1 (MQ=255)
                                                             cGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCTATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACctgctg                                                                       >  1:764705/1‑136 (MQ=255)
                                                             cGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACctgctg                                                                       >  1:495398/1‑136 (MQ=255)
                                                                gCAGGCGAGCCCGCAGCGCTTGCGTGTCGCTGCCGCAATCACGGGCCGGAAAGATCAGTGCCACGCCGGAAAAGCGGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACg                                                                    <  2:495398/136‑1 (MQ=255)
                                                                gCAGGCGAGCCAGCTGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCATGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACg                                                                    <  2:84255/136‑1 (MQ=255)
                                                                  aGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGc                                                                  >  1:196548/1‑136 (MQ=255)
                                                                  aGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGc                                                                  >  1:433159/1‑136 (MQ=255)
                                                                          cagatgcgcTGGGGGGTCGCGGCCCCCATCCGGGGCCGGAATCATCTGTGGCACGCCGCAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTg                                                          <  1:684866/131‑1 (MQ=255)
                                                                           cgcagcGCTTGCGTGTCGCGGCCGCCATCACGGGCCGGCAATATCAGTGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGa                                                         <  2:651944/135‑1 (MQ=255)
                                                                           agcagcGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGa                                                         <  1:450852/136‑1 (MQ=255)
                                                                                   ttGCGTGTCGCTGCCGCCATCTCGGGCCGGCAAGATCAGTGGCATGCCGGCAAAGCCGCCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATc                                                 <  2:287389/136‑1 (MQ=255)
                                                                                     gtgtgtCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCaa                                               <  2:219110/134‑1 (MQ=255)
                                                                                     gCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGCCCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGTCAGGCGCTGCCACCGATCaa                                               >  2:744134/1‑136 (MQ=255)
                                                                                     gCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGGCGGCAGGCGCTGACACCGATCaa                                               >  2:567422/1‑136 (MQ=255)
                                                                                     gCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCaa                                               >  2:678888/1‑136 (MQ=255)
                                                                                     gCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCaa                                               <  2:73804/136‑1 (MQ=255)
                                                                                     gCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCaa                                               <  2:214019/136‑1 (MQ=255)
                                                                                     gCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCaa                                               >  2:320515/1‑136 (MQ=255)
                                                                                                      tccggggCCGGCATGATCTGTGGCACGCGGTAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                              <  1:729671/132‑1 (MQ=255)
                                                                                                      tCAGGGGCCGGCATGATCAGTGGAAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                              <  2:511823/136‑1 (MQ=255)
                                                                                                      tCAGGGGCCGGCAAGATCTGTGGCAAGCCGGAAAAGCCGACCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                              <  1:830976/136‑1 (MQ=255)
                                                                                                      tCACGGGCCGGCAATATCAGTGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                              <  2:538223/136‑1 (MQ=255)
                                                                                                      tCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                              <  2:421645/136‑1 (MQ=255)
                                                                                                      tCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                              <  2:433159/136‑1 (MQ=255)
                                                                                                      tCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                              <  2:792943/136‑1 (MQ=255)
                                                                                                      tCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAACCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                              <  2:461231/136‑1 (MQ=255)
                                                                                                      tCACGGGCCGGCAAGATCAGTGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                              <  1:785681/136‑1 (MQ=255)
                                                                                                      tCACGGGCCGGAATGATCTGTGGCAAGCCGCAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                              <  2:693796/136‑1 (MQ=255)
                                                                                                      tCACGGGCCGGAAATATCAGTGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                              <  2:749111/136‑1 (MQ=255)
                                                                                                      ccccGGGCCGGCAAGACCTGTGGCACGCCGGAAAACCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                              <  1:654844/133‑1 (MQ=255)
                                                                                                      ccaggggCCGCCAATACCTGTGCCACCCCGCAAACGCGGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGTTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                              <  2:802562/132‑1 (MQ=255)
                                                                                                      acACGGGCCGGAAAGATCTGTGCCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                              <  2:565149/135‑1 (MQ=255)
                                                                                                      acACGGCCGGGAAATATCAGTGGAAAGCCGCAATAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTgcgc                              <  2:288646/135‑1 (MQ=255)
                                                                                                        ccGGGCCGGCAATATCAGTGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTc                            <  1:538950/135‑1 (MQ=255)
                                                                                                             ccGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATggg                       <  2:87774/136‑1 (MQ=255)
                                                                                                                gCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGGCGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCtt                    <  1:567422/136‑1 (MQ=255)
                                                                                                                actgtaTCCGTGCCACGCCGCAAACGCGGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCtt                    <  2:777298/131‑1 (MQ=255)
                                                                                                                         gTGGCAAGCCGGAAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCgtgt           <  1:20930/136‑1 (MQ=255)
                                                                                                                            gCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCa        <  2:175686/136‑1 (MQ=255)
                                                                                                                                  cGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCAggctgt  >  1:777298/1‑132 (MQ=255)
                                                                                                                                  cGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCAGGTAcc  >  1:370189/1‑136 (MQ=255)
                                                                                                                                     |                                                                                                                                    
CTGGGGGCGCCATTGTAATGATGCGGCTCGATATTTGAAGGCGAGCCGGTAAAAACAATGCCGTGCAGGCGAGCCAGCAGCGCTTGCGTGTCGCTGCCGCCATCACGGGCCGGCAAGATCAGTGGCAAGCCGGCAAAGCCGGCCGCCTCGACATACTTGTCGCCGACCGTGTGCGAGGAGTTCTTCCCTACCTGCTGACGGCAGGCGCTGACACCGATCAAAGGGACCGCATTTGCGCTCATGGGCTTACACCGTGTGCAGGTACC  >  NC_002947/6047847‑6048112

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