Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,621,0210A.45.1% 90.0 / 87.4 51coding (31/417 nt)PP_5654phosphotriesterase
Reads supporting (aligned to +/- strand):  ref base A (12/16);  new base . (11/12);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 7.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AACTCTCTCGACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATGGTAAAGAAA  >  NC_002947/4620875‑4621167
                                                                                                                                                  |                                                                                                                                                  
aaCTCTCTCGACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCaa                                                                                                                                                   >  4:363050/1‑148 (MQ=255)
     tctcGACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCaaaaaaa                                                                                                                                              >  3:48735/1‑148 (MQ=255)
           ccTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAATCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCaaaaaaaaacaaa                                                                                                                                        <  7:285132/148‑5 (MQ=255)
           ccTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCaaaaaaaaaaacaa                                                                                                                                       <  4:48735/149‑5 (MQ=255)
            cTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAg                                                                                                                                      >  4:31289/1‑148 (MQ=255)
              tCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAGCATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGAc                                                                                                                                    >  7:149384/1‑148 (MQ=255)
               ccAACAATCCAAGTTACACCCAATTATTTAAAAAAATCAAACAATAAAGCCCAACCAAAACTCACAAAAACCCCCGTAAAACTCATAATTTGAACATTCAAATGAACAAACCTAAAATCTGGATATGGTCCAAAAAAAAAACAAAGAc                                                                                                                                    <  8:112009/148‑1 (MQ=255)
                  agaagCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACa                                                                                                                                <  6:154659/149‑1 (MQ=255)
                      ccgtaGTTAGAGGGAATTATTTAAAAAAATCAAACAATAAAGCCCTATCAAAACTCAGACAAACTCCACTAAAACTCATATGTTGTACATGCATATGAACAAAGCTAACTTCTTGAGTTGTTCCAAAAAAAAAACAAAGACTACATTaa                                                                                                                            <  4:67612/145‑1 (MQ=255)
                       cGAAGTTAGAGCTAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGTACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTaaa                                                                                                                           <  8:33076/149‑1 (MQ=255)
                        gAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAATTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACCAAAGACTACATTaaa                                                                                                                           <  7:224097/148‑1 (MQ=255)
                        gAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAg                                                                                                                          >  4:296992/1‑148 (MQ=255)
                         aaGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAGCATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGaa                                                                                                                        <  8:149384/149‑1 (MQ=255)
                                 gCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATaaaa                                                                                                                 <  6:128994/149‑1 (MQ=255)
                                     attattTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAgg                                                                                                            <  3:363050/149‑1 (MQ=255)
                                     attattTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAg                                                                                                             <  3:31289/148‑1 (MQ=255)
                                      ttattTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGga                                                                                                           <  3:58513/149‑1 (MQ=255)
                                         tttAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACAAAAAAGCAg                                                                                                             >  6:129178/1‑145 (MQ=255)
                                             aaaTAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCa                                                                                                              >  7:207200/1‑139 (MQ=255)
                                             aaaTAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCa                                                                                                              <  8:207200/139‑1 (MQ=255)
                                             aaaTAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCAAAAAAAAAAACAAAGACTAAATTAAAGAACATAAAAGAAGGAGATAtt                                                                                                     >  8:305609/1‑149 (MQ=255)
                                                   ccAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGc                                                                                               <  4:136720/148‑1 (MQ=255)
                                                     aaaCAATAAAGCCCAAGCAAAACTCAAACAAAATACAGTAAAACTCAGAAGTGTAACATGCAAATGAACAAAGCTAACATCTTGAGATTGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCtt                                                                                             <  7:346748/148‑1 (MQ=255)
                                                       aCAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCa                                                                                           >  3:395446/1‑149 (MQ=255)
                                                          aTAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCa                                                                                       <  3:296992/149‑1 (MQ=255)
                                                           tAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGACATATTTTCTGCTTCACCCAg                                                                                      >  7:16176/1‑149 (MQ=255)
                                                                cccAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCAAAAAAAAAAACAAAGACTACATTAAACAACATAAAAGAAGGATGAATTTTGTCCTTCACCCAgg                                                                                     >  4:113640/1‑146 (MQ=255)
                                                                            tCATACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAAAAAAGAATAAATTAAAGAACAAAAAAGCAGAAGATATTTAGTGCTTCACCCACGAAGACGGAACaaat                                                                       >  3:266485/1‑148 (MQ=255)
                                                                                  aaaCTCCCGTAAAACTCATAAGTGTAACATGAAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTg                                                                <  7:305609/149‑1 (MQ=255)
                                                                                      tCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGc                                                             >  3:356415/1‑148 (MQ=255)
                                                                                          gTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAaata                                                        >  8:109929/1‑148 (MQ=255)
                                                                                                cTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATc                                                   >  1:72474/1‑148 (MQ=255)
                                                                                                cTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTg                                                                                                >  5:418927/1‑102 (MQ=255)
                                                                                                cTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTg                                                                                                <  6:418927/102‑1 (MQ=255)
                                                                                                 tCAGAAGTGGAACATGCAAATGAACAAAGCTAAAATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGACCATAAACGCAGGAGATAATTTCTGCTAAACCCATGAACACGGAACATATTATTTTGGGCAAATAATATCa                                                  >  6:158681/1‑148 (MQ=255)
                                                                                                  cAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCaa                                                 >  5:227166/1‑148 (MQ=255)
                                                                                                  cAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGAATACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCAAGGAAGACGGAACATATTATTTTGGGGAAATAATATCaa                                                 >  8:180195/1‑148 (MQ=255)
                                                                                                   aGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCaaaa                                               <  2:72474/149‑1 (MQ=255)
                                                                                                    taatttgAAAATTCAAATTAAAAAATCTAACATATGTATCTGTTACAAAAAAAAAACAAATACTACATTAAAGAACATAAAATCAGTATATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCaaaa                                               <  7:260962/142‑1 (MQ=255)
                                                                                                     aaGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAaat                                             >  8:307921/1‑148 (MQ=255)
                                                                                                     aaGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAaat                                             >  3:186521/1‑148 (MQ=255)
                                                                                                          taaCATGCAAATTAACAAAGCTAAAATCTGTAGATGGTCAAAAAAAAAAACAAAGACTACATTAAAGAACTTAAAAGCAGTAGATATTTTCTGCTTCACCCAGGAAGACGGAACATTTTTTTTTGCGCAAATAATATCAAAAATAATc                                         <  5:158681/147‑1 (MQ=255)
                                                                                                          gAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGAAATTTTCTGCTTCACCCACGAAGACGGAACATATTATTTTTGGAGAATAATATCAAAaaaaata                                         >  4:83564/1‑147 (MQ=255)
                                                                                                            aCATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACt                                      <  4:395446/149‑1 (MQ=255)
                                                                                                            aCATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTg                                     <  7:307921/149‑1 (MQ=255)
                                                                                                                 cAAATGAAGAAAGGTAACATCTGGAGATGGTCC‑AAAATATAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTGACCCAGGAAGACGGAACATATTATTTTGGGCAAATATTATCATAAATTATCTCTGGcga                                 <  8:16176/148‑3 (MQ=255)
                                                                                                                 cAAATGAACATATCTAACTTCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGTACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGTcc                                  <  2:129187/148‑1 (MQ=255)
                                                                                                                 cAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGcc                                  <  3:146098/148‑1 (MQ=255)
                                                                                                                   aaTGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGAAAATAATATCAAAAATAATCACTGGGCATg                               >  5:79311/1‑149 (MQ=255)
                                                                                                                    aTGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATATTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATg                               >  3:209631/1‑147 (MQ=255)
                                                                                                                             gCTAAATTCTGGATATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAat                     <  3:113640/149‑1 (MQ=255)
                                                                                                                                        gagaTGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATgg          <  4:186521/148‑1 (MQ=255)
                                                                                                                                          gaTGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATGGTa        >  7:48506/1‑148 (MQ=255)
                                                                                                                                            tGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAGATAATCACTGGCCATGTAGCTGAAATATTCAGAATGGTaaa      <  7:109929/148‑1 (MQ=255)
                                                                                                                                                 cAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATGGTaaagaaa  <  5:129178/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AACTCTCTCGACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATGGTAAAGAAA  >  NC_002947/4620875‑4621167

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