Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,563,208 A→G 9.9% L220S (TTG→TCG)  PP_3146 ← oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,563,2080AG 9.9% 128.5 / 2.6 48L220S (TTG→TCG) PP_3146oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (27/16);  new base G (2/3);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTCGGGCAAGGGTTCGGTGGCCACCACATAGCTGCCGGCCGGCAGTACCTTGCCGCTCAGGCGCGGTTCCAGTTCATCCAGGTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCAGCGGGTGCAGGTGGCCCGAGCCCATGTCGACCAGGCCGCCAGCGTACTGGTCG  >  NC_002947/3563061‑3563357
                                                                                                                                                   |                                                                                                                                                     
gcTTCGGGCAAGGGTTCGGTGGCCACCACATAGCTGCCGGCCGGCAGTACCTTGCCGCTCAGGCGCGGTTCCAGTTCATCCAGGTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAAc                                                                                                                                                     <  2:216664/150‑1 (MQ=255)
      ggCAAGGGTTCGGTGGCCACCACATAGCTGCCGGCCGGCAGTACCTTGCCGCTCAGGCGCGGTTCCAGTTCATCCAGGTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGc                                                                                                                                               >  2:54099/1‑150 (MQ=255)
      ggCAAGGGTTCGGTGGCCACCACATAGCTGCCGGCCGGCAGTACCTTGCCGCTCAGGCGCGGTTCCAGTTCATCCAGGTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGc                                                                                                                                               >  2:54197/1‑150 (MQ=255)
       gCAAGGGTTCGGTGGCCACCACATAGCTGCCGGCCGGCAGTACCTTGCCGCTCAGGCGCGGTTCCAGTTCATCCAGGTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGcc                                                                                                                                              >  1:606145/1‑150 (MQ=255)
       gCAAGGGTTCGGTGGCCACCACATAGCTGCCGGCCGGCAGTACCTTGCCGCTCAGGCGCGGTTCCAGTTCATCCAGGTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGcc                                                                                                                                              >  1:525812/1‑150 (MQ=255)
                     gCCACCACATAGCTGCCGGCCGGCAGTACCTTGCCGCTCAGGCGCGGTTCCAGTTCATCCAGGTGAGAGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATgcgc                                                                                                                                <  2:606145/150‑1 (MQ=255)
                          cacaTAGCTGCCGGCCGGCAGTACCTTGCCGCTCAGGCGCGGTTCCAGTTCATCCAGGTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCAc                                                                                                                           >  1:1001688/1‑150 (MQ=255)
                                      ggccggcAGTACCTTGCCGCTCAGGCGCGGTTCCAGTTCATCCATGTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTc                                                                                                               <  2:421300/150‑1 (MQ=255)
                                           gcAGTACCTTGCCGCTCAGGCGCGGTTCCAGTTCATCCAGGTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGaaaa                                                                                                          >  2:400325/1‑150 (MQ=255)
                                                aCCTTGCCGCTCAGGCGCGGTTCCAGTTCATCCAGGTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCGACCt                                                                                                                                                   <  1:108590/104‑1 (MQ=255)
                                                aCCTTGCCGCTCAGGCGCGGTTCCAGTTCATCCAGGTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCGACCt                                                                                                                                                   >  2:108599/1‑104 (MQ=255)
                                                aCCTTGCCGCTCAGGCGCGGTTCCAGTTCATCCAGGTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCGACCt                                                                                                                                                   >  2:108590/1‑104 (MQ=255)
                                                aCCTTGCCGCTCAGGCGCGGTTCCAGTTCATCCAGGTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCGACCt                                                                                                                                                   <  1:108599/104‑1 (MQ=255)
                                                 ccTTGCCGCTCAGGCGCGGTTCCAGTTCATCCAGGTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCa                                                                                                    >  1:268252/1‑150 (MQ=255)
                                                                gcgGTTCCAGTTCATCCAGGTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTggg                                                                                     >  2:861789/1‑150 (MQ=255)
                                                                       cAGTTCATCCAGGTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGGGCAGCGTGAAAATGCGCACGCCCAGGCTGTGGGCcgcgcg                                                                              >  1:607910/1‑150 (MQ=255)
                                                                       cAGTTCATCCAGGTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCcgcgcg                                                                              >  2:1107476/1‑150 (MQ=255)
                                                                         gTTCATCCAGGTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGTACGCCCAGGCTGTGGGCCGCGCggg                                                                            >  2:258239/1‑150 (MQ=255)
                                                                             aTCCAGGTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTc                                                                        >  2:356676/1‑150 (MQ=255)
                                                                               ccAGGTGAGCGTGGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCGACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGAGCTTCGc                                                                      <  2:442050/150‑1 (MQ=255)
                                                                                cAGGTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGcc                                                                     >  1:810857/1‑150 (MQ=255)
                                                                                  ggTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCt                                                                   >  1:1058488/1‑150 (MQ=255)
                                                                                   gTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTT‑CCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGc                                                                 >  2:650894/1‑150 (MQ=255)
                                                                                          ttGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGa                                                           <  2:787235/150‑1 (MQ=255)
                                                                                          ttGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGa                                                           <  2:1001688/150‑1 (MQ=255)
                                                                                              aGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCa                                                       >  1:969102/1‑150 (MQ=255)
                                                                                                cccAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCAGc                                                     >  1:511815/1‑150 (MQ=255)
                                                                                                cccAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTAGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCAGc                                                     >  1:786863/1‑150 (MQ=255)
                                                                                                  cAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCAGCgg                                                   >  2:1032869/1‑150 (MQ=255)
                                                                                                        cAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTc                                                                        <  2:909888/123‑1 (MQ=255)
                                                                                                        cAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTc                                                                        >  1:909888/1‑123 (MQ=255)
                                                                                                        cAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCAGCGGGTGCAg                                             <  1:632423/150‑1 (MQ=255)
                                                                                                         aGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCAGCGGGTGCAgg                                            <  1:1032869/150‑1 (MQ=255)
                                                                                                          ggCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCAGCGGGTGCAGGt                                           >  2:979362/1‑150 (MQ=255)
                                                                                                           gctgctTGCGCGCACCTTGCCGCGCGCCGTGAGCAGCGTCAACCGGGGCCCTTGCTAAATAGGCAGCACGTGGCTTTATCCATTGAAGAACAATTTTAGTCGGTGGGCCGCGCGGGCTTCGCCCTGCCCCAGATCCAGCGGGTGCAGGTg                                          >  1:252600/1‑150 (MQ=255)
                                                                                                            ctgctTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCAGCGGGTGCAGGTgg                                         <  1:400325/150‑1 (MQ=255)
                                                                                                                ttGCGCGCACCTT‑CCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCAGCGGGTGCAGGTGGCCCGa                                    <  1:650894/150‑1 (MQ=255)
                                                                                                                    gcgcACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCAGCGGGTGCAGGTGGCCCGAGcc                                 >  1:553566/1‑150 (MQ=255)
                                                                                                                     cgcACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCAGCGGGTGCAGGTGGCCCGAGccc                                <  2:1058488/150‑1 (MQ=255)
                                                                                                                      gcACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCAGCGGGTGCAGGTGGCCCGAGCCCa                               >  2:935989/1‑150 (MQ=255)
                                                                                                                          cTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCAGCGGGTGCAGGTGGCCCGAGCCCATGTc                           <  2:305812/150‑1 (MQ=255)
                                                                                                                           ttGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCAGCGGGTGCAGGTGGCCCGAGCCCATGTCg                          >  2:740262/1‑150 (MQ=255)
                                                                                                                                  gcgcCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCAGCGGGTGCAGGTGGCCCGAGCCCATGTCGACCAGgc                   <  2:717474/150‑1 (MQ=255)
                                                                                                                                     ccGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCAGCGGGTGCAGGTGGCCCGGGCCCATGTCGACCAGgccgc                >  2:373921/1‑150 (MQ=255)
                                                                                                                                     ccGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGAGCGGGCTTCGCCCTGCACCAGATCCAGCGGGTGCAGGTGGCCCGAGCCCATGTCGACCAGgccgc                <  2:268252/150‑1 (MQ=255)
                                                                                                                                          atctgCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCAGCGGGTGCAGGTGGCCCGAGCCCATGTCGACCAGGCCGCCAGCg           <  2:992947/146‑1 (MQ=255)
                                                                                                                                            cAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCAGCGGGTGCAGGTGGCCCGAGCCCATGTCGACCAGGCCGCCAGCGTa         <  1:97690/150‑1 (MQ=255)
                                                                                                                                               cGNCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTNGTTCGANAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCAGCGGGTGCAGGTGGCCCGAGCCCATGTCGACCAGGCCGCCAGCGTagat      >  1:992947/1‑147 (MQ=255)
                                                                                                                                                   aaCCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCAGCGGGTGCAGGTGGCCCGAGCCCATGTCGACCAGGCCGCCAGCGTACTGGTCg  >  2:399257/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCTTCGGGCAAGGGTTCGGTGGCCACCACATAGCTGCCGGCCGGCAGTACCTTGCCGCTCAGGCGCGGTTCCAGTTCATCCAGGTGAGCGTTGCAGCCCAGCACCAGGCTGCTTGCGCGCACCTTGCCGCGCGCCGTGTGCAGCGTCAACCTGGGCCCGTGCTCAATGCGCAGCACCGGGCTTTGTTCGAAAATGCGCACGCCCAGGCTGTGGGCCGCGCGGGCTTCGCCCTGCACCAGATCCAGCGGGTGCAGGTGGCCCGAGCCCATGTCGACCAGGCCGCCAGCGTACTGGTCG  >  NC_002947/3563061‑3563357

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