Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,265,163 G→T 8.2% P86Q (CCG→CAG)  PP_1106 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,265,1630GT 8.2% 152.3 / 4.1 49P86Q (CCG→CAG) PP_1106hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/17);  new base T (2/2);  total (30/19)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.95e-01

GTCGACTCGGTGATGAAGGTGTGGCACGCCACCGGTTCGAAGGGCGTACAGGTGCCATCGGGTTCTACCTTGTAGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCGGGGCTGGCGGTCAGATAATGGCCATTGCCGGTGCGTGCGT  >  NC_002947/1265016‑1265307
                                                                                                                                                   |                                                                                                                                                
gTCGACTCGGTGATGAAGGTGTGGCACGCCACCTGTTCGAAGGGCGTACAGGTGCCATCGGGTTCTACCTTGTAGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCccggccgg                                                                                                                                                 <  2:218083/149‑1 (MQ=255)
gTCGACTCGGTGATGAAGGTGTGGCACGCCACCGGTTCGAAGGGCGTACAGGTGCCATCGGGTTCTACCTTGTAGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCccggccgg                                                                                                                                                 <  7:21362/149‑1 (MQ=255)
   gACTCGGTGATGAAGGTGTGGCACGCCACCGGTTCGAAGGGCGTACAGGTGCCATCGGGTTCTACCTTGTAGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGt                                                                                                                                               <  7:277441/148‑1 (MQ=255)
    aCTCGGTGATGAAGGTGTGGCACGCCACCGGTTCGAAGGGCGTACAGGTGCCATCGGGTTCTACCTTGTAGTCACCGGATGCAACCCAGACTTCACCTCGCTACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGc                                                                                                                                             >  4:216539/1‑149 (MQ=255)
               aaGGTGTGGCACGCCACCGGTTCGAAGGGCGTACAGGTGCCATCGGGTTCTACCTTGTAGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACg                                                                                                                                   <  4:265572/148‑1 (MQ=255)
                           gCCACCGGTTCGAAGGGCGTACAGGTGCCATCGGGTTCTACCTTGTGGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCgtggtg                                                                                                                      <  1:5639/149‑1 (MQ=255)
                              aCCGGTTCGAAGGGCGTACAGGTGCCATCGGGTTCTACCTCGTAGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGTGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTgcag                                                                                                                   <  5:78074/149‑1 (MQ=255)
                                   ttCGAAGGGCGTACAGGTGCCATCGGGTTCTACCTTGTAGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGgc                                                                                                              >  6:134242/1‑149 (MQ=255)
                                         gggCGTACAGGTGCCATCGGGTTCTACCTTGTAGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCAc                                                                                                        <  7:93194/149‑1 (MQ=255)
                                              tACAGGTGCCATCGGGTTCTACCTTGTAGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCgtggtg                                                                                                                      <  7:13184/130‑1 (MQ=255)
                                              tACAGGTGCCATCGGGTTCTACCTTGTAGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCgtggtg                                                                                                                      >  8:13184/1‑130 (MQ=255)
                                                    tGCCATCGGGTTCTACCTTGTAGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCTGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAg                                                                                                   <  5:3445/143‑1 (MQ=255)
                                                    tGCCATCGGGTTCTACCTTGTAGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCTGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAg                                                                                                   >  4:44774/1‑143 (MQ=255)
                                                    tGCCATCGGGTTCTACCTTGTAGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCTGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAg                                                                                                   >  6:3445/1‑143 (MQ=255)
                                                    tGCCATCGGGTTCTACCTTGTAGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCTGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAg                                                                                                   <  3:44774/143‑1 (MQ=255)
                                                         tCGGGTTCTACCTTGTAGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGGAGGGTCTg                                                                                        >  8:266462/1‑149 (MQ=255)
                                                         tCGGGTTCTACCTTGTAGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTg                                                                                        >  1:68890/1‑149 (MQ=255)
                                                         tCGGGTTCTACCTTGTAGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTg                                                                                        <  2:240881/149‑1 (MQ=255)
                                                              ttCTACCTTGTAGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGAt                                                                                   >  4:30673/1‑149 (MQ=255)
                                                                 tACCTTGTAGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGCGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGa                                                                                >  8:174909/1‑149 (MQ=255)
                                                                     ttGTAGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTc                                                                            >  4:161929/1‑149 (MQ=255)
                                                                           tCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGcc                                                                      >  6:75992/1‑149 (MQ=255)
                                                                           tCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCAGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGcc                                                                      >  5:5143/1‑149 (MQ=255)
                                                                                      aaCCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCGGCAGCTTGAAGTCCTGGCCCAGGCGGCTgc                                                           >  2:221987/1‑149 (MQ=255)
                                                                                              cTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATg                                                   >  7:172046/1‑149 (MQ=255)
                                                                                                    cTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCa                                             >  6:273086/1‑149 (MQ=255)
                                                                                                      cGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCAcc                                           >  5:76888/1‑149 (MQ=255)
                                                                                                       gATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCg                                                >  8:8683/1‑143 (MQ=255)
                                                                                                       gATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCg                                                <  7:8683/143‑1 (MQ=255)
                                                                                                       gATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCAc                                            <  1:221987/147‑1 (MQ=255)
                                                                                                       gATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCg                                          >  7:289639/1‑149 (MQ=255)
                                                                                                         tACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCggg                                        <  3:285683/149‑1 (MQ=255)
                                                                                                           cTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCgggg                                       <  8:172046/148‑1 (MQ=255)
                                                                                                            tCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCGGGGCt                                     >  8:75880/1‑149 (MQ=255)
                                                                                                              cAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCGGGGCTgg                                   >  2:86722/1‑149 (MQ=255)
                                                                                                               aGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCGGGGCTGGc                                  <  5:134242/149‑1 (MQ=255)
                                                                                                                  atctccTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCGGGGCTGGCGGt                               <  3:176448/145‑1 (MQ=255)
                                                                                                                    cACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCGCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCGGGGCTGGCGGTCa                             <  3:216539/149‑1 (MQ=255)
                                                                                                                      ccTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCGGGGCTGGCGGTCAGa                           <  3:30673/149‑1 (MQ=255)
                                                                                                                       cTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCGGGGCTGGCGGt                               >  4:176448/1‑144 (MQ=255)
                                                                                                                       cTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCGGGGCTGGCGGTCAGa                           >  7:35961/1‑148 (MQ=255)
                                                                                                                       cTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCGGGGCTGGCGGTCAGa                           >  2:136721/1‑148 (MQ=255)
                                                                                                                            ccGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCGGGGCTGGCGGTCAGATAATgg                     >  2:52682/1‑149 (MQ=255)
                                                                                                                             cGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCGGGGCTGGCGGTCAGATAATg                      <  3:254953/147‑1 (MQ=255)
                                                                                                                             cGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCGGGGCTGGCGGTCAGATAATg                      >  4:254953/1‑147 (MQ=255)
                                                                                                                             cGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCGGGGCTGGCGGTCAGATAATGGc                    <  2:68890/149‑1 (MQ=255)
                                                                                                                                  ccAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCGGGGCTGGCGGTCAGATAATGGCCATTg               >  2:204096/1‑149 (MQ=255)
                                                                                                                                  ccAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCGGGGCTGGCGGTCAGATAATGGCCATTg               >  4:292254/1‑149 (MQ=255)
                                                                                                                                       aCGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCGGGGCTGGCGGTCAGATAATGGCCATTGCCGgt          >  4:277864/1‑149 (MQ=255)
                                                                                                                                         gTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCTCCGGGGCTGGCGGTCAGATAATGGCCATTGCCGgtgc        >  7:288136/1‑149 (MQ=255)
                                                                                                                                         gTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCGGGGCTGGCGGTCAGATAATGGCCATTGCCGgtgc        >  8:125843/1‑149 (MQ=255)
                                                                                                                                               ggccggGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCGGGGCTGGCGGTCAGATAATGGCCATTGCCGGTgcgtgcgt  >  5:288168/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GTCGACTCGGTGATGAAGGTGTGGCACGCCACCGGTTCGAAGGGCGTACAGGTGCCATCGGGTTCTACCTTGTAGTCACCGGATGCAACCCAGACTTCACCTCGATACTCCAGGCGCACCTGTGCCGAACCCAGCACGTGCCCGGCCGGGTGCAGGCTCAACGTCACGCCGTGGTGCAGCAGGCGCTCACCGTAGGGCAGGGTCTGCAGATTGATGTCCTGGCCCAGGCGGCTGCGCAGAATGCCGGCACCGGGGCTGGCGGTCAGATAATGGCCATTGCCGGTGCGTGCGT  >  NC_002947/1265016‑1265307

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