Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,127,349 T→A 9.3% V907E (GTG→GAG)  PP_1887 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,127,3490TA 9.3% 114.6 / 3.0 43V907E (GTG→GAG) PP_1887hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (18/21);  new base A (2/2);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

CAGCATCAAGGTGGGGAACACCGAAAGTGCCGCCGTGTTGTTGGGCGAACAGGTGTTCGATGGCCTTGAGCGGCGAACCGTAGCCATTACCGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAGGAAGCGAAATACGAGTATGACCTGAAAAATGCCAGGCTTACCCATTGCGAAGAGCCTGGGGAT  >  NC_002947/2127203‑2127497
                                                                                                                                                  |                                                                                                                                                    
cAGCATCAAGGTGGGGAACACCGAAAGTGCCGCCGTGTTGTTGGGCGAACAGGTGTTCGATGGCCTTGAGCGGCGAACCGTAGCCATTACCGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTg                                                                                                                                                     >  7:218491/1‑148 (MQ=255)
  gCATCAAGGTGGGGAACACCGAAAGTGCCGCCGTGTTGTTGGGCGAACAGGTGTTCGATGGCCTTGAGCGGCGAACCGTAGCCATTACCGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAg                                                                                                                                                  <  2:19867/149‑1 (MQ=255)
        aGGTGGGGAACACCGAAAGTGCCGCCGTGTTGTTGGGCGAACAGGTGTTCGATGGCCTTGAGCGGCGAACCGTAGCCATTACCGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCg                                                                                                                                            >  6:64653/1‑149 (MQ=255)
          gTGGGGAACACCGAAAGTGCCGCCGTGTTGTTGGGCGAACAGGTGTTCGATGGCCTTGAGCGGCGAACCGTAGCCATTACCGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGa                                                                                                                                          >  4:295758/1‑149 (MQ=255)
           tGGGGAACACCGAAAGTGCCGCCGTGTTGTTGGGCGAACAGGTGTTCGATGGCCTTGAGCGGCGAACCGTAGCCATTACCGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGAt                                                                                                                                         <  7:44459/149‑1 (MQ=255)
                     cGAAAGTGCCGCCGTGTTGTTGGGCGAACAGGTGTTCGATGGCCTTGAGCGGCGAACCGTAGCCATTACCGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGTCGCCGGATGGTACCCGGa                                                                                                                               >  4:242128/1‑149 (MQ=255)
                      gAAAGGGCCGCCGGGTTGCTGGGCGAACAGGGGGTCGATGGCCTTGAGCGGCGAACCGTCGCCATTACCGGAGGGCGAGAGCAGACATTCGAATATGACCCACGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGAt                                                                                                                              <  8:218491/149‑1 (MQ=255)
                             ccgccgTGTTGTTGGGCGAACAGGTGTTCGATGGCCTTGAGCGGCGAACCGTAGCCATTACCGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGagtac                                                                                                                       >  5:128934/1‑149 (MQ=255)
                                                aCAGGTGTTCGATGGCCTTGAGCGGCGAACCGTAGCCATTACCGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTgg                                                                                                    <  1:227311/149‑1 (MQ=255)
                                                aCAGGTGTTCGATGGCCTTGAGCGGCGAACCGTAGCCATTACCGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTgg                                                                                                    <  2:237781/149‑1 (MQ=255)
                                                    gtgtTCGATGGCCTTGAGCGGCGAACCGTAGCCATTACCGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAgg                                                                                                                                                 <  1:83211/100‑1 (MQ=255)
                                                    gtgtTCGATGGCCTTGAGCGGCGAACCGTAGCCATTACCGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAgg                                                                                                                                                 >  2:83211/1‑100 (MQ=255)
                                                     tgtTCGATGGCCTTGAGCGGCGAACCGTAGCCATTACCGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCg                                                                                                        <  3:200050/140‑1 (MQ=255)
                                                     tgtTCGATGGCCTTGAGCGGCGAACCGTAGCCATTACCGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCg                                                                                                        >  4:200050/1‑140 (MQ=255)
                                                     tgtTCGATGGCCTTGAGCGGCGAACCGTAGCCATTACCGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGaa                                                                                               <  3:295758/149‑1 (MQ=255)
                                                     tgtTCGATGGCCTTGAGCGGCGAACCGTAGCCATTACCGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGaa                                                                                               <  5:64653/149‑1 (MQ=255)
                                                       ttCGATGGCCTTGAGCGGCGAACCGTAGCCATAACCGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGaa                                                                                               >  1:37065/1‑147 (MQ=255)
                                                           aTGGCCTTGAGCGGCGAACCGTAGCCATTACCGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGccc                                                                                         <  6:128934/149‑1 (MQ=255)
                                                                cTTGAGCGGCGAACCGTAGCCATTACCGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCt                                                                                    >  7:145388/1‑149 (MQ=255)
                                                                     gcggcgAACCGTAGCCATTACCGGAGGGCGAGAGCAGACATTCGAATATGACGCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCTTTGGGTGAAGAGCCCGTCCTTCgg                                                                                <  3:242128/148‑1 (MQ=255)
                                                                                    cATTACCGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAg                                                                 >  5:196951/1‑148 (MQ=255)
                                                                                          cGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAGGAAGCg                                                           >  8:20285/1‑148 (MQ=255)
                                                                                               ggCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCgg                                                                                >  8:152696/1‑122 (MQ=255)
                                                                                               ggCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCgg                                                                                <  7:152696/122‑1 (MQ=255)
                                                                                                 cGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAGGAAGCGAAATACGa                                                   <  8:145388/149‑1 (MQ=255)
                                                                                                 cGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAGGAAGCGAAATACGa                                                   <  6:196951/149‑1 (MQ=255)
                                                                                                 cGAGAGCAGACATTCGAATATGACCCAGGGGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAGGAAGCGAAATACGa                                                   <  8:152901/149‑1 (MQ=255)
                                                                                                        aGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGAGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAGGaa                                                              >  1:220147/1‑131 (MQ=255)
                                                                                                        aGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGAGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAGGaa                                                              <  2:220147/131‑1 (MQ=255)
                                                                                                                gAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAGGAAGCGAAATACGAGTATGACCTGaaaaa                                    >  2:86690/1‑149 (MQ=255)
                                                                                                                   tatGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAGGAAGCGAAATACGAGTATGACCTGAAAAATGc                                 >  7:217160/1‑149 (MQ=255)
                                                                                                                      gACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAGGAAGCGAAATACGa                                                   >  7:224346/1‑128 (MQ=255)
                                                                                                                      gACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAGGAAGCGAAATACGa                                                   <  8:224346/128‑1 (MQ=255)
                                                                                                                       aCCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAGGAAGCGAAATACGAGTATGACCTGAAAAATGCCAg                              <  8:217160/148‑1 (MQ=255)
                                                                                                                               gAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTg                                                                                                     <  7:56827/69‑1 (MQ=255)
                                                                                                                               gAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTg                                                                                                     >  8:56827/1‑69 (MQ=255)
                                                                                                                                aGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAGGAAGCGAAATACGAGTATGACCTGAAAAATGCCAGGCTTACCCAt                    <  2:37065/149‑1 (MQ=255)
                                                                                                                                aGCGGCAACCGCATTGGGAGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAGGAAGCGAAATACGAGTATGACCTGAAAAATGCCAGGCTTACCCAt                    >  8:124546/1‑149 (MQ=255)
                                                                                                                                    gCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAGGAAGCGAAATACGAGTATGACCTGAAAAATGCCAGGc                            >  2:200874/1‑137 (MQ=255)
                                                                                                                                    gCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAGGAAGCGAAATACGAGTATGACCTGAAAAATGCCAGGc                            <  1:200874/137‑1 (MQ=255)
                                                                                                                                      aaCCGCATTGGGAGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAGGAAGCGAAATACGAGTATGACCTGAAAAATGCCAGGCTTACCCATTGCGaa              <  7:124546/149‑1 (MQ=255)
                                                                                                                                                ggTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAGGAAGCGAAATACGAGTATGACCTGAAAAATGCCAGGCTTACCCATTGCGAAGAGCCTg       <  6:37865/146‑1 (MQ=255)
                                                                                                                                                ggTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAGGAAGCGAAATACGAGTATGACCTGAAAAATGCCAGGCTTACCCATTGCGAAGAGCCTg       >  5:37865/1‑146 (MQ=255)
                                                                                                                                                  tGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAGGAAGCGAAATACGAGTATGACCTGAAAAATGCCAGGCTTACCCATTGCGAAGAGCCTGGGGAt  >  5:299648/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CAGCATCAAGGTGGGGAACACCGAAAGTGCCGCCGTGTTGTTGGGCGAACAGGTGTTCGATGGCCTTGAGCGGCGAACCGTAGCCATTACCGGAGGGCGAGAGCAGACATTCGAATATGACCCAGGTGAGCGGCAACCGCATTGGGTGAAGGCGCCGGATGGTACCCGGATCGAGTACCAGTACCGGCCAGCGTTGGGTGAAGAGCCCGTCCTTCGGGTCCTCTCCGGCAAGGAAGCGAAATACGAGTATGACCTGAAAAATGCCAGGCTTACCCATTGCGAAGAGCCTGGGGAT  >  NC_002947/2127203‑2127497

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