Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,595,3160CA37.9% 83.8 / 15.9 58F129L (TTC→TTA) PP_3171hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (13/23);  new base A (22/0);  total (35/23)
Fisher's exact test for biased strand distribution p-value = 2.73e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCTGTTCGGTGCTGACCTGGCGGCGA  >  NC_002947/3595181‑3595449
                                                                                                                                        |                                                                                                                                     
gTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTc                                                                                                                                       <  1:196028/136‑1 (MQ=255)
    gtggtggCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGcc                                                                                                                                   >  2:829670/1‑136 (MQ=255)
     gggtggCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCa                                                                                                                                  <  2:352000/135‑1 (MQ=255)
     gcggggCCGAACAGGGTGGCACCCGCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCa                                                                                                                                  <  2:994118/132‑1 (MQ=255)
                 tttttGGCACCCTTAATTCTTTGCACCCTGGTCGTATCGCCTGGGGCGTGGGCGGTGCGCCTGGCACTCCGGGGCCGACGGTTGGGGCTTTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGa                                                                                                                      <  2:994277/134‑1 (MQ=255)
                                aTACATTGCACCCAGGTCGTATCGCCCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCCCCGCGCCAGA‑GCGGGATTTAGGCCGGGACCGCGCCGGGCGGCGTGATTAccg                                                                                                       >  1:239831/1‑135 (MQ=255)
                                aTACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCCCGCCACC‑CCGGGCTTTAGGCCAGGACAACGCCGCGCTGCGTGATTaaca                                                                                                       >  1:473263/1‑133 (MQ=255)
                                         aCCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGGGCGGCGCCGCGCCCCA‑GCGGGCTTTAGGGCCGGCCCTCCCCGCGCGGCGGGATTTCCAGGcccgccc                                                                                              >  1:965966/1‑131 (MQ=255)
                                         aCCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCCGGACATCGCCGAGCTGCGTGATTACCTGGCCGCCaa                                                                                              >  2:504736/1‑136 (MQ=255)
                                         aCCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATAGCCGCGCTGCGGGATTACCAGGGCGACaa                                                                                              >  1:238814/1‑136 (MQ=255)
                                         aCCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTAGGCCCCGACATCGCCGCGCTGCGGGATTTCCAGGCCGCCaa                                                                                              >  2:783202/1‑136 (MQ=255)
                                         aCCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCCCCCCAGA‑GCGGGCTTTAGGCCAGGACAAGGCCGCGCTGGGGGATTAACAGGCCGCCaa                                                                                              >  1:511561/1‑136 (MQ=255)
                                         aCCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCCGCCCCCCCCCACACCCGGCATTACGG‑CAGGGCAAAGCCGGGCTGGGTGATTACCAGGcccccca                                                                                              >  1:506618/1‑131 (MQ=255)
                                          cccAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGCCCGTGCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGTTTTCGGGCCGGACATCGCCGCGCTGCGTGATTTCCTGGCCgcccac                                                                                             >  2:245509/1‑131 (MQ=255)
                                           cccggTCGTAGCGACCGGGGCCTGGGCCGTGCGCCAGGCACGGCCGGGCCGACCGGTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACg                                                                                            <  2:750617/133‑1 (MQ=255)
                                           ccAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGGCCCCGA‑GCGGGGTTTAAGCCAGGCCATAGCCGAGCTGCGTGGTTTCCTGGCCGACAAcc                                                                                            >  1:559831/1‑135 (MQ=255)
                                           ccAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACg                                                                                            <  2:504914/136‑1 (MQ=255)
                                                tCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCACA‑GCGGGATTTAGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCaaa                                                                                       >  2:579616/1‑136 (MQ=255)
                                                tCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCCCGCCAGA‑GCGGGATTTAGGCCAGGACATAGCCGAGCTGCGTGCTTACCTGGCCGACAACGGCaaa                                                                                       >  2:492564/1‑136 (MQ=255)
                                                     tCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGCCCGTTCGGGCATTGCGCGGCGCCGCGCCCGC‑GCGGGGTTTAGGCCCGGGCCAAGCCGCGCGGCGGGGTTTACTGGCCGCCCACCGCAAACGccc                                                                                  >  2:999667/1‑135 (MQ=255)
                                                                ctgggccGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCCGC‑GCGGGGTTTAGGGCAGGACAACGCCGCGCTGCGTGATTACCAGGGCGCCCACGGCAAACCCCGGGTGCgcggcg                                                                       >  1:109259/1‑136 (MQ=255)
                                                                ctgggccGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCCGA‑GCGGGATTTAGGCCCGGCCATCGCCGAGCGGGGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCgcggcg                                                                       >  1:245219/1‑136 (MQ=255)
                                                                ctgggccGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGC‑GCGGGATTTAGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACCGCAAACGCCGGGTTCgccgcg                                                                       >  2:736766/1‑136 (MQ=255)
                                                                ctgggccGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCCCCCCGCCCGA‑GCGGGATTTAGGGCAGGACATAGCCCCGCGGCGTGATTTACATGCCCCCAACCGCCAAACCCCGGTGCgcgggg                                                                       >  2:160441/1‑134 (MQ=255)
                                                                ctgggccGGGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCgcggcg                                                                       <  2:807649/136‑1 (MQ=255)
                                                                   ggccGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTAc                                                                    <  2:760911/136‑1 (MQ=255)
                                                                    gccGTGCGCCAGGCACGGCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTAcc                                                                   <  1:207339/136‑1 (MQ=255)
                                                                      cGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTAGGCCAGGACATCGCCGCGCTGCGTGATTACCTGGCCGCCAACGGCAAACGCCGGGTGCGCGGCGTACCgg                                                                 >  1:214102/1‑136 (MQ=255)
                                                                        cgcgcCAGCCACGGCCGCCCCGACCGTTCGGGCAATGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGGGACTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGt                                                               <  2:668128/135‑1 (MQ=255)
                                                                          cgcCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTAGGCCAGGCCATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTgc                                                             >  2:934868/1‑136 (MQ=255)
                                                                           gcCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCCGC‑GCGGGGTTTCGGCCAGGACCTCGCCGCGCGGCGGGATTACCTGGCCGCCAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTgcg                                                            >  1:565059/1‑136 (MQ=255)
                                                                           gcCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGCGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTgcg                                                            >  2:21515/1‑136 (MQ=255)
                                                                           gcCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTgcg                                                            >  1:850982/1‑136 (MQ=255)
                                                                           gcCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTgcg                                                            >  1:649154/1‑136 (MQ=255)
                                                                           gcCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTAGGCCAGGCCATCGCCGCGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTTCCGGGTgcg                                                            >  2:915194/1‑136 (MQ=255)
                                                                           gcCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTAGGCCAGGACATCGCCGAGCTGGGTGGTTTCCTGGCCGCCAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTgcg                                                            >  1:655714/1‑136 (MQ=255)
                                                                           gcCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTAGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTgcg                                                            >  2:836933/1‑136 (MQ=255)
                                                                           gcCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCCCGCCAGA‑GCGGGATTTAGGCCAGGACAAAGGCGCGCTGCGGGATTACCATGCCGCCCACCGCAAACGCCGGGTGCGCGGCGTACCGGGTgcc                                                            >  2:103551/1‑135 (MQ=255)
                                                                           gcCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCCCGCCAGA‑GCCGGATTTAGGCCAGGACATAGCCCAGCTGCGTGATTAACAGGCCGACAACCGCAAACGCCCGGTGCGCGGCGTACCGGGTgcc                                                            >  1:117419/1‑135 (MQ=255)
                                                                                     tCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGc                                                  <  2:493369/136‑1 (MQ=255)
                                                                                       cGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGCGCTGCGTGTTTTCCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGCTGTGCCg                                                >  2:741974/1‑136 (MQ=255)
                                                                                       cGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCg                                                >  2:880097/1‑136 (MQ=255)
                                                                                       cGGGCCGACCGTTCGGGCATTGCGCCGCCCCGCGCCAGA‑GCGGGATTTAGGCCACGACAACGCCGAGCAGCGGGATTACCTGGGCGACAACCGCAAACGCCCGGTGCGCGGCGTACCGGGTGCGCACCATGTGCCg                                                >  2:33852/1‑136 (MQ=255)
                                                                                            cGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGtgtg                                           >  2:497078/1‑136 (MQ=255)
                                                                                               ccGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGAt                                        <  2:674516/136‑1 (MQ=255)
                                                                                                 gTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTAGGCCAAGACAAAGCCGAGCAGGGTGATTAACAGGCCCCCAACGGCAAAAGCCGGGTGGGGCGGGTACCGGGTGCGCACGATTTGCCCGTGTgggtgc                                      >  2:963299/1‑132 (MQ=255)
                                                                                                         tttGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTcc                              <  2:238814/135‑1 (MQ=255)
                                                                                                         aTTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTcc                              <  2:780658/136‑1 (MQ=255)
                                                                                                              gcggcgCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCGAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCt                         <  2:625778/136‑1 (MQ=255)
                                                                                                              gcggcgCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCt                         <  1:777003/136‑1 (MQ=255)
                                                                                                              gcggcgCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCt                         <  1:421723/136‑1 (MQ=255)
                                                                                                              gcggcgCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCt                         <  2:565059/136‑1 (MQ=255)
                                                                                                              gcggcgCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCt                         <  1:668230/136‑1 (MQ=255)
                                                                                                              gcggcgCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCt                         <  1:588485/136‑1 (MQ=255)
                                                                                                                       cgccAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCGGGTGTGGATCCTCGGGTCCAGCCTGGTCGGAGc                >  1:619600/1‑136 (MQ=255)
                                                                                                                                cGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCTGTTCGGTGCTGACCTgg        <  1:579616/136‑1 (MQ=255)
                                                                                                                                cGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTACGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCTGTTCGGTGCTGACCTgg        <  2:494930/136‑1 (MQ=255)
                                                                                                                                 gggATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCTGTTCGGTGCTGACCTggc       <  2:473263/136‑1 (MQ=255)
                                                                                                                                      ttCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCTGTTCGGTGCTGACCTGGCGGCGa  <  2:559831/136‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                     
GTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGA‑GCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCTGTTCGGTGCTGACCTGGCGGCGA  >  NC_002947/3595181‑3595449

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: ‑