Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,188,263 G→T 7.5% V53V (GTG→GTT)  PP_3684 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,188,2630GT 7.5% 172.5 / 3.3 59V53V (GTG→GTT) PP_3684hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/29);  new base T (2/2);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

CGCGTGGTTTAACATCCACCATGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCT  >  NC_002947/4188118‑4188402
                                                                                                                                                 |                                                                                                                                           
cgcgTGGTTTAACATCCACCATGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGaa                                                                                                                                           <  2:265116/148‑1 (MQ=255)
   gTGGTTTAACATCCACCATGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGc                                                                                                                                       <  8:249849/149‑1 (MQ=255)
      gTTTAACATCCACCATGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCTGTTCACCGCCTc                                                                                                                                    <  6:191877/149‑1 (MQ=255)
       tttAACATCCACCATGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATACGGCATGAGAGCCGATTCATCGAGGCGGTCAACAGCCTc                                                                                                                                    >  4:191651/1‑148 (MQ=255)
       tttAACATCCACCATGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTc                                                                                                                                    >  2:161186/1‑148 (MQ=255)
                 accaTGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCAATATACAAACTTTGGGATTACGTCTTCCCATACAGAAGACCAGAAGCAGAGATAATGCATGAGCGCCAATTCAACGAGGCTGTGAACAGCAACAAAACCTTGAg                                                                                                                         >  3:92375/1‑149 (MQ=255)
                  ccaTGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAgg                                                                                                                        <  8:43904/149‑1 (MQ=255)
                   caTGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAggg                                                                                                                       >  6:324982/1‑149 (MQ=255)
                         cAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACt                                                                                                                  >  8:112728/1‑148 (MQ=255)
                         cAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGAAAGAGCGCCAATTAATAGAGGCGGTAAAAAGCCTCAAAACCTTGAGGGTCaat                                                                                                                  >  7:284387/1‑146 (MQ=255)
                            aCCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCaaaa                                                                                                                                >  6:18228/1‑131 (MQ=255)
                            aCCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCaaaa                                                                                                                                <  5:18228/131‑1 (MQ=255)
                               atccatACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAg                                                                                                           >  3:46599/1‑149 (MQ=255)
                                tccatACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAgg                                                                                                          <  3:94481/149‑1 (MQ=255)
                                  catACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGc                                                                                                         <  7:112728/148‑1 (MQ=255)
                                  catACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCg                                                                                                        >  7:41699/1‑149 (MQ=255)
                                   atACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCgg                                                                                                       >  2:122539/1‑149 (MQ=255)
                                   atACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCgg                                                                                                       <  7:116034/149‑1 (MQ=255)
                                     aCTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGCCCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGa                                                                                                     >  7:237820/1‑149 (MQ=255)
                                       tCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCgg                                                                                                       >  1:59124/1‑145 (MQ=255)
                                       tCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCgg                                                                                                       <  2:59124/145‑1 (MQ=255)
                                         ggCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATg                                                                                                   <  8:41699/147‑1 (MQ=255)
                                         ggCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGtc                                                                                                 <  8:122138/149‑1 (MQ=255)
                                         ggCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTAACCGAGGCGGGATGtc                                                                                                 >  4:125483/1‑149 (MQ=255)
                                         ggCAGGTGGAAGATGAAATCACTATAACAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAATGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATg                                                                                                   <  7:213880/147‑1 (MQ=255)
                                                  aaGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGAccc                                                                                        <  1:122539/149‑1 (MQ=255)
                                                      tGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGaa                                                                                    <  2:144605/149‑1 (MQ=255)
                                                            cACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGtcg                                                                                                <  8:225360/131‑1 (MQ=255)
                                                            cACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGtcg                                                                                                >  7:225360/1‑131 (MQ=255)
                                                                   ccAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCGATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGGCGATCGACCCTGAAGAGATACgtgat                                                                        <  3:191651/148‑1 (MQ=255)
                                                                   ccAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACgtgatg                                                                       <  8:38265/149‑1 (MQ=255)
                                                                   ccAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACgtgatg                                                                       <  8:43235/149‑1 (MQ=255)
                                                                             gATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGaa                                                                                    >  6:261994/1‑126 (MQ=255)
                                                                             gATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGaa                                                                                    <  5:261994/126‑1 (MQ=255)
                                                                                  cggcTTCCCAGACCGCAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGATGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGcc                                                        <  8:237820/146‑1 (MQ=255)
                                                                                        cccAGACAGAAGAACATAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAg                                                   <  8:284387/148‑1 (MQ=255)
                                                                                          cAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTTAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTgaga                                                           <  7:134530/138‑1 (MQ=255)
                                                                                          cAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTgaga                                                           >  8:134530/1‑138 (MQ=255)
                                                                                               cGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCaa                                           <  5:324982/149‑1 (MQ=255)
                                                                                                      cAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGATCAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTg                                    >  1:223444/1‑149 (MQ=255)
                                                                                                      cAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATc                                                                  >  6:11868/1‑119 (MQ=255)
                                                                                                      cAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATc                                                                  <  5:11868/119‑1 (MQ=255)
                                                                                                       aGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTAACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTgg                                   <  3:125483/149‑1 (MQ=255)
                                                                                                       aGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCAATGACCGAGGCCGCATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACttcg                                   >  1:215291/1‑147 (MQ=255)
                                                                                                         aaGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTc                                 >  7:289873/1‑149 (MQ=255)
                                                                                                          aGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTc                                 <  8:289873/148‑1 (MQ=255)
                                                                                                                      gCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATc                    >  3:307638/1‑149 (MQ=255)
                                                                                                                      gCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATc                    >  1:35328/1‑149 (MQ=255)
                                                                                                                             cgcTAATTCATCTAGGCGGTGTACAGCCTCAAATCCTTGAGGGTCAGTGACCGAGGCGGGATGTCGATCGACCCTGATGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGGTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGc             <  2:215291/149‑1 (MQ=255)
                                                                                                                                    tCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACTTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACaa       >  6:311779/1‑148 (MQ=255)
                                                                                                                                      tttgAGGCTGTTCACCGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAAc    >  8:287902/4‑149 (MQ=255)
                                                                                                                                      aTCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGCCCGAGGCGGGATGTCGCTCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAACGCGAACaaacc    >  1:246810/1‑147 (MQ=255)
                                                                                                                                         gAGGCTGTTCACCGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCt  >  6:201859/1‑148 (MQ=255)
                                                                                                                                         gAGGCTGTGAACAGCCTCAAAACCTTGAGGTTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCt  <  7:104232/148‑1 (MQ=255)
                                                                                                                                         gAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCt  <  7:108685/148‑1 (MQ=255)
                                                                                                                                              ggTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCaa                                           >  5:75676/1‑102 (MQ=255)
                                                                                                                                              ggTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCaa                                           <  6:75676/102‑1 (MQ=255)
                                                                                                                                              ggTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGt                                  >  8:231157/1‑111 (MQ=255)
                                                                                                                                              ggTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGt                                  <  7:231157/111‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
CGCGTGGTTTAACATCCACCATGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCT  >  NC_002947/4188118‑4188402

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