Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,706,587 T→A 8.8% L202Q (CTG→CAG)  PP_5009 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,706,5870TA 8.8% 165.0 / 5.5 57L202Q (CTG→CAG) PP_5009hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/24);  new base A (3/2);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

GGGCTCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCATGTC  >  NC_002947/5706439‑5706732
                                                                                                                                                    |                                                                                                                                                 
gggCTCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCt                                                                                                                                                   <  5:52877/149‑1 (MQ=255)
  gCTCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTgg                                                                                                                                                 <  3:315059/149‑1 (MQ=255)
    tCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCt                                                                                                                                               >  6:67259/1‑149 (MQ=255)
       acgaAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGTACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTc                                                                                                                                            >  3:249798/1‑149 (MQ=255)
       acgaAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTc                                                                                                                                            >  1:170823/1‑149 (MQ=255)
        cgaAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCAGGCTGTCg                                                                                                                                           <  3:88806/149‑1 (MQ=255)
         gaAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCAGGCTGTCgt                                                                                                                                          <  7:265467/149‑2 (MQ=255)
         gaAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAACGCTGCATTGCTCCAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCACTGCTCGACGGGGAACACGACACCCAGCCTCTCgt                                                                                                                                          >  8:265467/1‑148 (MQ=255)
                   aaGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACaa                                                                                                                                <  6:102446/149‑1 (MQ=255)
                      tctcGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGt                                                                                                                              >  6:98414/1‑148 (MQ=255)
                        tcGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAgg                                                                                                                                        <  8:266384/136‑1 (MQ=255)
                        tcGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAgg                                                                                                                                        >  7:266384/1‑136 (MQ=255)
                            gcATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTc                                                                                                                       <  5:131461/149‑1 (MQ=255)
                                          gAGCACACCATCACGGCTGCATTGCTCGAGCCCCTCCCTCGCCACGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGCAACGGATGGCgg                                                                                                         <  2:291004/149‑1 (MQ=255)
                                            gCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCgggg                                                                                                       <  3:273809/149‑1 (MQ=255)
                                              acacCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGca                                                                                                     >  8:7251/1‑149 (MQ=255)
                                               cacCATCAAGGCTGCATTGCTCGAGTACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGcac                                                                                                    <  4:249798/149‑1 (MQ=255)
                                                           tGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTAccc                                                                                        >  6:281661/1‑149 (MQ=255)
                                                                tGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCg                                                                                   >  3:144264/1‑149 (MQ=255)
                                                                 gCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCg                                                                                   <  3:321403/148‑1 (MQ=255)
                                                                  cTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGa                                                                                 <  1:129174/149‑1 (MQ=255)
                                                                   tCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGa                                                                                 >  3:67964/1‑148 (MQ=255)
                                                                     gagaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGAt                                                                                <  5:67259/147‑1 (MQ=255)
                                                                     gagaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGAt                                                                                <  7:292399/147‑1 (MQ=255)
                                                                         aCCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTg                                                                          <  5:281661/149‑1 (MQ=255)
                                                                             ccctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCtca                                                                      >  5:92101/1‑149 (MQ=255)
                                                                              cctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCtcat                                                                     <  2:237235/149‑1 (MQ=255)
                                                                                   ccAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATcacca                                                                >  3:153681/1‑149 (MQ=255)
                                                                                    cAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATcacca                                                                >  4:72013/1‑148 (MQ=255)
                                                                                          cTGGCAGATGATGCAAAACCTGCAATGTACGCCCTGCTCGACCGGGTACACGACAGCCTCGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACAGGATCTATTACTAATCACCACAAGTgg                                                         >  2:225849/1‑149 (MQ=255)
                                                                                           tGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGTTGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTgg                                                         >  4:318506/1‑148 (MQ=255)
                                                                                             gCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTAGACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGgcaa                                                      >  1:241691/1‑149 (MQ=255)
                                                                                                gatgatGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGgcaaccg                                                   <  6:92101/149‑1 (MQ=255)
                                                                                                   gatgCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGgcaaccgcaa                                                >  3:95231/1‑149 (MQ=255)
                                                                                                    atgCAAAAGCTGCATTGTACGCCCTTCTCGACGGTGTACACGACAGCCTGGCTGTCGATGCGATACAAGTCACGCTCGGAACTGATGGCGGGGCACACACCGTTACCCTAGCGGATCTATTGCTCATCACCACAAGTGgcaaccgcaa                                                <  1:225849/148‑1 (MQ=255)
                                                                                                          aaGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGtt                                         <  4:95231/149‑1 (MQ=255)
                                                                                                          aaGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGtt                                         <  2:241691/149‑1 (MQ=255)
                                                                                                            gCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCt                                       >  1:37478/1‑149 (MQ=255)
                                                                                                              tGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAACTCACGCTCGGAATGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCtg                                      >  4:11552/1‑148 (MQ=255)
                                                                                                                cAATGTACGCCCTGCTCGAGGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGTTGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCtgcaa                                   <  3:318506/149‑1 (MQ=255)
                                                                                                                 aaTGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCtgcaat                                  <  4:67964/149‑1 (MQ=255)
                                                                                                                 aaTGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCtgcaat                                  <  7:7251/149‑1 (MQ=255)
                                                                                                                 aaTGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCtgcaat                                  >  6:93962/1‑149 (MQ=255)
                                                                                                                  aTGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTg                                                          <  6:272638/124‑1 (MQ=255)
                                                                                                                  aTGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTg                                                          >  5:272638/1‑124 (MQ=255)
                                                                                                                   tGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCtgcaatg                                 >  4:269024/1‑148 (MQ=255)
                                                                                                                      aCGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCTGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCtgcaatgcaa                              <  3:72013/148‑1 (MQ=255)
                                                                                                                             gCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGcg                      >  7:147144/1‑149 (MQ=255)
                                                                                                                                       gTACACGACAGCCTGGCTGTCGAGGCGATACAACTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGgttcggttcg            >  6:11215/1‑149 (MQ=255)
                                                                                                                                        tACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTtcgt           >  2:5046/1‑149 (MQ=255)
                                                                                                                                         acacGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACAACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTtcgtc          >  3:140390/1‑149 (MQ=255)
                                                                                                                                          ctcGACAGCCAGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTtcgtcg         >  1:85979/3‑149 (MQ=255)
                                                                                                                                          ctcGACAGCCAGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTtcgtcg         >  8:293621/3‑149 (MQ=255)
                                                                                                                                          cacGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTtcgtc          <  2:88038/148‑1 (MQ=255)
                                                                                                                                           acGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGc        >  5:193817/1‑149 (MQ=255)
                                                                                                                                                 gCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCATGTc  <  5:11215/149‑1 (MQ=255)
                                                                                                                                                  ccTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCATGTc  >  3:321439/1‑148 (MQ=255)
                                                                                                                                                  ccTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCATGTc  >  2:22332/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGGCTCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCATGTC  >  NC_002947/5706439‑5706732

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