Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 733,852 C→T 8.0% intergenic (+19/+172) PP_t11 → / ← ndh tRNA‑Lys/NADH dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947733,8520CT 8.0% 170.0 / 2.6 59intergenic (+19/+172)PP_t11/ndhtRNA‑Lys/NADH dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (29/24);  new base T (2/3);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCGGGTTCACCCGCGAATGCGGCAGTGCAGCCACCATCGCATTCACAGGTAAACCCG  >  NC_002947/733712‑733998
                                                                                                                                            |                                                                                                                                                  
ggAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGc                                                                                                                                           <  1:324041/150‑1 (MQ=255)
       gATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCttttt                                                                                                                                    >  2:618366/1‑150 (MQ=255)
          gTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCtt                                                                                                                                 >  2:94910/1‑150 (MQ=255)
          gTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCtt                                                                                                                                 >  1:433646/1‑150 (MQ=255)
            cAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCtttt                                                                                                                               <  1:94910/150‑1 (MQ=255)
              ggAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGc                                                                                                                             <  1:772483/150‑1 (MQ=255)
               gAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGcc                                                                                                                            >  2:857235/1‑150 (MQ=255)
                  ttCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCg                                                                                                                         >  2:556735/1‑150 (MQ=255)
                  ttCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCATGCGAGATGCCCtttg                                                                                                                                     >  2:1189263/1‑137 (MQ=14)
                  ttCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCATGCGAGATGCCCtttg                                                                                                                                     <  1:1189263/138‑2 (MQ=14)
                     gAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGaaa                                                                                                                      >  1:395019/1‑150 (MQ=255)
                      aaTCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGaaaa                                                                                                                     >  1:525576/1‑150 (MQ=255)
                          cccTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATgc                                                                                                                 <  1:1010821/150‑1 (MQ=255)
                           cttaaCCCCGATCATTTGTGGGTCGTTAGCTCATTTGTTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAACCCTACACTACCCACCATGAACAAAGGGCATCTCGTATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATgcg                                                                                                                <  2:351366/146‑1 (MQ=21)
                           ccTCACCCCGACCATTTTTTTTTTGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATgcg                                                                                                                >  2:933773/1‑150 (MQ=255)
                            cTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGcc                                                                                                                                          >  2:940516/1‑123 (MQ=255)
                            cTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGcc                                                                                                                                          <  1:940516/123‑1 (MQ=255)
                                  ccGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCg                                                                                                         <  2:433646/150‑1 (MQ=255)
                                     aCCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCtttt                                                                                                                               <  2:250220/125‑1 (MQ=255)
                                     aCCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCtttt                                                                                                                               >  1:250220/1‑125 (MQ=255)
                                        aTTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCaa                                                                                                   <  1:857235/150‑1 (MQ=255)
                                                gTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCt                                                                                           >  2:635625/1‑150 (MQ=255)
                                                           aGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAACAATGCGCCCCTCGCCGCAATCCCACCTGCATCCATTTg                                                                                >  2:27365/1‑150 (MQ=255)
                                                            gTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGGCCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGc                                                                               >  1:91623/1‑150 (MQ=255)
                                                            gTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGc                                                                               <  2:525576/150‑1 (MQ=255)
                                                              tGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCca                                                                             <  2:127795/150‑1 (MQ=255)
                                                                gTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCttt                                                                                                                                <  1:522474/97‑1 (MQ=255)
                                                                gTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCttt                                                                                                                                >  2:522474/1‑97 (MQ=255)
                                                                       aGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGAATGAGATGCCCTTTTTCTTTTGCCTCGAAACATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGccgc                                                                    >  1:820453/1‑148 (MQ=255)
                                                                              ttttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGc                                                                               <  2:574630/132‑1 (MQ=255)
                                                                              ttttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGc                                                                               >  1:574630/1‑132 (MQ=255)
                                                                              ttttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAg                                                             >  1:1056642/1‑150 (MQ=255)
                                                                               tttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCAGCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAgg                                                            <  2:310300/150‑1 (MQ=255)
                                                                                    ccAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTAttt                                                                                 <  1:1011836/124‑1 (MQ=255)
                                                                                    ccAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTAttt                                                                                 >  2:1011836/1‑124 (MQ=255)
                                                                                    ccAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCgg                                                       >  1:1151173/1‑150 (MQ=255)
                                                                                     cAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCggg                                                      >  2:238290/1‑150 (MQ=255)
                                                                                     cAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCggg                                                      <  2:1151173/150‑1 (MQ=255)
                                                                                     cAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCggg                                                      <  1:933773/150‑1 (MQ=255)
                                                                                         tGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCGGGTTCa                                                  <  1:618366/150‑1 (MQ=255)
                                                                                           gTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCAcc                                                                                            >  2:1063056/1‑106 (MQ=255)
                                                                                           gTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCAcc                                                                                            <  1:1063056/106‑1 (MQ=255)
                                                                                            tCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCGGGTTCAccc                                               >  2:501304/1‑150 (MQ=255)
                                                                                            tCGTAGGTGCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCAGCCTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCGGGTTCAccc                                               <  1:27365/150‑1 (MQ=255)
                                                                                                   ttCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCGGGTTCACCCGCGAATg                                        >  1:86821/1‑150 (MQ=255)
                                                                                                    tCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCGGGTTCACCCGCGAATGc                                       >  2:361800/1‑150 (MQ=255)
                                                                                                     cGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCGGGTTCACCCGCGAATGCg                                      <  2:86821/150‑1 (MQ=255)
                                                                                                                    ccaccaTGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCCTCGGCGCAATCCCAACTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCGGGTTCACCCGCGAATGc                                       <  1:349481/134‑1 (MQ=255)
                                                                                                                    ccaccaTGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCGGGTTCACCCGCGAATGc                                       >  2:349481/1‑134 (MQ=255)
                                                                                                                    ccaccaTGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCGGGTTCACCCGCGAATGCGGCAGTGCAGccacca                       >  1:387128/1‑150 (MQ=255)
                                                                                                                      accaTGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCGGGTTCACCCGCGAATGCGGCAGTGCAGCCACCATc                     <  2:387128/150‑1 (MQ=255)
                                                                                                                      accaTGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCGGGTTCACCCGCGAATGCGGCAGTGCAGCCACCATc                     >  1:814559/1‑150 (MQ=255)
                                                                                                                      accaTGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCGGGTTCACCCGCGAATGCGGCAGTGCAGCCACCATc                     >  1:493339/1‑150 (MQ=255)
                                                                                                                      accaTGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCGGGTTCACCCGCGAATGCGGCAGTGCAGCCACCATc                     <  2:1056642/150‑1 (MQ=255)
                                                                                                                          tGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCGGGTTCACCCGCGAATGCGGCAGTGCAGCCACCATCGCAt                 <  2:493339/150‑1 (MQ=255)
                                                                                                                            gaaaaaGGGCATCTCATGCGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCCATGTAGGAGCGGGTTCACCCGCGAATGCGGCAGTGCAGCCACCATc                     >  2:95284/4‑144 (MQ=255)
                                                                                                                            gaaaaaGGGCATCTCATGCGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCCATGTAGGAGCGGGTTCACCCGCGAATGCGGCAGTGCAGCCACCATc                     <  1:95284/141‑1 (MQ=255)
                                                                                                                                   ggCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCGGGTTCACCCGCGAATGCGGCAGTGCAGCCACCATCGCATTCACAGGTa        >  2:631274/1‑150 (MQ=255)
                                                                                                                                         tcGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCGGGTTGGCCCGCGAATGCGGCAGTGCAGCCACCGTCGCGTTGGCAGGTAAACCcg  >  1:89054/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
GGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAGGGCATCTCGCATGAGATGCCCTTTTTCTTTTGCCTCGAAAAATGCGCCCATCGGCGCAATCCCACCTGCATCTATTTGCCACACCGCCCCTATGTAGGAGCGGGTTCACCCGCGAATGCGGCAGTGCAGCCACCATCGCATTCACAGGTAAACCCG  >  NC_002947/733712‑733998

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