Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,022,727 G→C 7.9% L293L (CTG→CTC PP_5273 → oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,022,7270GC 7.9% 153.5 / 4.7 52L293L (CTG→CTCPP_5273oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (26/21);  new base C (2/2);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

CAATCGCGAGCTGGGCGGCAAGGTGCTGCCGGCCGGCAGCTACATCATCGCCACCGAGCCGCTGGGCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACATCGCCGCCTACATGCGGCCGAAGATGCTCAAGGTGTTCCCGCAGCTTGCCGACGTGCGCATCGACTACCAGTGGGGTGGCATGA  >  NC_002947/6022583‑6022875
                                                                                                                                                |                                                                                                                                                    
cAATCGCGAGCTGGGCGGCAAGGTGCTGCCGGCCGGCAGCTACATCATCGCCACCGAGCCGCTGGGCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTCTccgc                                                                                                                                                 >  1:142014/1‑150 (MQ=255)
cAATCGCGAGCTGGGCGGCAAGGTGCTGCCGGCCGGCAGCTACATCATCGCCACCGAGCCGCTGGGCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTCTccgc                                                                                                                                                 >  1:142024/1‑150 (MQ=255)
 aaTCGCGAGCTGGGCGGCAAGGTGCTGCCGGCCGGCAGCTACATCATCGCCACCGAGCCGCTGGGCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTCTccgcc                                                                                                                                                <  2:142024/150‑1 (MQ=255)
 aaTCGCGAGCTGGGCGGCAAGGTGCTGCCGGCCGGCAGCTACATCATCGCCACCGAGCCGCTGGGCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTCTccgcc                                                                                                                                                <  2:142014/150‑1 (MQ=255)
   tCGCGAGCTGGGCGGCAAGGTGCTGCCGGCCGGCAGCTACATCATCGCCACCGAGCCGCTGGGCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGa                                                                                                                                              >  1:238335/1‑150 (MQ=255)
                        gctgcCGGCCGGCAGCTACATCATCGCCACCGAGCCGCTGGGCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCgg                                                                                                                         >  2:497542/1‑150 (MQ=255)
                         ctgcCGGCCGGCAGCTACATCATCGCCACCGAGCCGCTGGGCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGt                                                                                                                        >  2:374685/1‑150 (MQ=255)
                         ctgcCGGCCGGCAGCTACATCATCGCCACCGAGCCGCTGGGCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTATCCGCCGACCACCGCCTGCTGTTTGGCGGt                                                                                                                        <  2:586060/150‑1 (MQ=255)
                             cggccggcAGCTACATCATCGCCACCGAGCCGCTGGGCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGtgcct                                                                                                                    <  2:176726/150‑1 (MQ=255)
                                        tACATCATCGCCACCGAGCCGCTGGGCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCa                                                                                                         <  2:185507/150‑1 (MQ=255)
                                           atcatcGCCACCGAGCCGCTGGGCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGc                                                                                                      <  1:497542/150‑1 (MQ=255)
                                              atcGCCACCGAGCCGCTGGGCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCcgc                                                                                                   >  1:118642/1‑150 (MQ=255)
                                              atcGCCACCGAGCCGCTGGGCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCcgc                                                                                                   >  1:226001/1‑150 (MQ=255)
                                                  ccACCGAGCCGCTGGGCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGa                                                                                                 >  2:399895/1‑148 (MQ=255)
                                                  ccACCGAGCCGCTGGGCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGa                                                                                                 <  1:399895/148‑1 (MQ=255)
                                                   cACCGAGCCGCTGGGCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGAccc                                                                                              <  1:681765/150‑1 (MQ=255)
                                                      cGAGCCGCTGGGCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGaa                                                                                           >  1:188677/1‑150 (MQ=255)
                                                            gCTGGGCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACAt                                                                                     >  2:206508/1‑150 (MQ=255)
                                                             cTGGGCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACATc                                                                                    <  1:70879/150‑1 (MQ=255)
                                                                 gCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGtgcctgc                                                                                                                  <  1:465693/116‑1 (MQ=255)
                                                                 gCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGtgcctgc                                                                                                                  >  2:465693/1‑116 (MQ=255)
                                                                 gCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACATcgccg                                                                                >  1:752464/1‑150 (MQ=255)
                                                                    aagaaCGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACATcgccgcaa                                                                             <  2:733167/150‑3 (MQ=255)
                                                                     agaaCGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCg                                                                                                                          <  1:719478/104‑1 (MQ=255)
                                                                     agaaCGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCg                                                                                                                          >  2:719478/1‑104 (MQ=255)
                                                                      gaaCGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGCCCCGCAGGCCATCGCCGCCTAc                                                                           >  2:325674/1‑150 (MQ=255)
                                                                          gcgcACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGc                                                                                                      >  2:277915/1‑119 (MQ=255)
                                                                          gcgcACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGc                                                                                                      <  1:277915/119‑1 (MQ=255)
                                                                             cacgcacATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACATCGCCGCCTACATGCGGc                                                                    <  2:118642/150‑1 (MQ=255)
                                                                              acgcacATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACATc                                                                                    >  1:139920/1‑133 (MQ=255)
                                                                              acgcacATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACATc                                                                                    <  2:139920/133‑1 (MQ=255)
                                                                              acgcacATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACATCGCCGCCTACATGCGGcc                                                                   >  1:710593/1‑150 (MQ=255)
                                                                                    aTTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCACGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACATCGCCGCCTACATGCGGCCGAAGAt                                                             >  1:661968/1‑150 (MQ=255)
                                                                                             gCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAg                                                                                                        >  1:85860/1‑98 (MQ=255)
                                                                                             gCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAg                                                                                                        <  2:85860/98‑1 (MQ=255)
                                                                                                   tatGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACATCGCCGCCTACATGCGGCCGAAGATGCTAACGGTGTTccc                                              >  2:655332/1‑150 (MQ=255)
                                                                                                     tGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACATCGCCGCCTACATGCGGCCGAAGATGCTCAAGGTGTTCCCGc                                            >  1:358553/1‑150 (MQ=255)
                                                                                                         ggTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCAGCAGGCCGTCGAAGTGTAAATGGGCTTGAAGAGGATCGAGGGGGTCCCGCAGCt                                        >  2:485024/1‑150 (MQ=255)
                                                                                                          gtgtgtGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGCCGTACATCGCCGCCTACATGCGGCCGAAGATGCTCAAGGTGTTCCCGCAGCtt                                       <  1:655332/150‑1 (MQ=255)
                                                                                                                    aGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAgg                                                                                                       <  1:131557/76‑1 (MQ=255)
                                                                                                                    aGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAgg                                                                                                       >  2:131557/1‑76 (MQ=255)
                                                                                                                    aGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCAGAAGGACATCGCCGCCTACATGCGGCCGAAGATGCTCAAGGTGTTCCCGCAGCTTGCCGACGTgc                             <  1:325674/150‑1 (MQ=255)
                                                                                                                            gcgcTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACATCGCCGCCTACATGCGGCCGAAGATGCTCAAGGTGTTCCCGCAGCTTGCCGACGTGCGCATCGAc                     <  1:374685/150‑1 (MQ=255)
                                                                                                                                 ggACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACATCGCCGCCTACATGCGGCCGAAGATGCTCAAGGTGTTCCCGCAGCTTGCCGACGTggttga                         >  1:773085/1‑136 (MQ=255)
                                                                                                                                 ggACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACATCGCCGCCTACATGCGGCCGAAGATGCTCAAGGTGTTCCCGCAGCTTGCCGACGTggttga                         <  2:773085/141‑6 (MQ=255)
                                                                                                                                       ctacCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACATCGCCGCCTACATGCGGCCGAAGATGCTCAAGGTGTTCCCGCAGCTTGCCGACGTGCGCATCGACTACCAGTgg            >  2:487001/1‑148 (MQ=255)
                                                                                                                                       ctacCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACATCGCCGCCTACATGCGGCCGAAGATGCTCAAGGTGTTCCCGCAGCTTGCCGACGTGCGCATCGACTACCAGTgg            <  1:487001/148‑1 (MQ=255)
                                                                                                                                        tacCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACATCGCCGCCTACATGCGGCCGAAGATGCTCAAGGTGTTCCCGCAGCTTGCCGACGTGCGCATCGACTACCAGTGGGGt         >  1:140074/1‑150 (MQ=255)
                                                                                                                                          ccGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACATCGCCGCCTACATGCGGCCGAAGATGCTCAAGGTGTTGCCGCAGCTTGCCg                                   >  2:544802/1‑122 (MQ=255)
                                                                                                                                          ccGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACATCGCCGCCTACATGCGGCCGAAGATGCTCAAGGTGTTCCCGCAGCTTGCCg                                   <  1:544802/122‑1 (MQ=255)
                                                                                                                                            gCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACATCGCCGCCTACATGCGGCCGAAGATGCTCAAGGTGTTCCCGCAGCTTGCCGACGTGCGCATCGACTACCAGTGGGGTGGCa     <  2:710593/150‑1 (MQ=255)
                                                                                                                                               tGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGGAGCGACGCGAAGGACATCGCCGCCTACATGCGGCCGAAGATGCTCAAGGTGTTCCCGCAGCTTGCCGACGTGCGCATCGACTACCAGTGGGGTGGCATGa  >  1:454551/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CAATCGCGAGCTGGGCGGCAAGGTGCTGCCGGCCGGCAGCTACATCATCGCCACCGAGCCGCTGGGCGAAGAACGCGCACGCACATTGCTGCCGCAGAATATGGCGGTGTGTGACCAGCGCGTGGCGCTGGACTACTACCGCCTGTCCGCCGACCACCGCCTGCTGTTTGGCGGTGCCTGCCACTATTCAGGCCGCGACCCGAAGGACATCGCCGCCTACATGCGGCCGAAGATGCTCAAGGTGTTCCCGCAGCTTGCCGACGTGCGCATCGACTACCAGTGGGGTGGCATGA  >  NC_002947/6022583‑6022875

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: