Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 4,537,135 C→G L96L (CTG→CTC)  RPPX_20270 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099744,537,1350CG97.9% 116.6 / NA 47L96L (CTG→CTC) RPPX_20270hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (25/21);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGTTTGTCCGGGCCGTGGTCCATGAAGTGGCGGAGGTATTCCCATACCTCGGCGACATGCTTGACGCCCGGTTCGGTGGGTTCGTTGCT  >  CP009974/4536987‑4537278
                                                                                                                                                    |                                                                                                                                               
ggcAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGTGGTGAg                                                                                                                                               >  2:538085/1‑151 (MQ=255)
ggcAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAg                                                                                                                                               <  1:6298794/151‑1 (MQ=255)
ggcAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAg                                                                                                                                               >  2:3858373/1‑151 (MQ=255)
 gcAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAgg                                                                                                                                              <  2:643669/151‑1 (MQ=255)
 gcAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAgg                                                                                                                                              <  1:3079513/151‑1 (MQ=255)
  cAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGc                                                                                                                                             <  1:3307879/151‑1 (MQ=255)
  cAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGc                                                                                                                                             >  2:4181657/1‑151 (MQ=255)
  cAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGc                                                                                                                                             >  2:2377390/1‑151 (MQ=255)
  cAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGc                                                                                                                                             <  1:5139082/151‑1 (MQ=255)
   aGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCa                                                                                                                                            <  2:5519146/151‑1 (MQ=255)
   aGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCa                                                                                                                                            <  2:390495/151‑1 (MQ=255)
   aGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCa                                                                                                                                            <  1:2553206/151‑1 (MQ=255)
    gATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAg                                                                                                                                           >  1:3810490/1‑151 (MQ=255)
    gATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAg                                                                                                                                           >  2:1647394/1‑151 (MQ=255)
    gATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAg                                                                                                                                           >  2:1042791/1‑151 (MQ=255)
     aTGCGATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGa                                                                                                                                          >  2:5382319/1‑151 (MQ=255)
     aTGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGa                                                                                                                                          >  2:169197/1‑151 (MQ=255)
     aTGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGa                                                                                                                                          >  1:306132/1‑151 (MQ=255)
      tGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGCACAGCGCGGCCCAGGACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGAt                                                                                                                                         <  2:1785304/151‑1 (MQ=255)
      tGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGAt                                                                                                                                         <  1:1770845/151‑1 (MQ=255)
      tGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGAt                                                                                                                                         <  2:1786368/151‑1 (MQ=255)
       gCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATg                                                                                                                                        >  1:4896614/1‑151 (MQ=255)
       gCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATg                                                                                                                                        <  2:2857561/151‑1 (MQ=255)
       gCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATg                                                                                                                                        >  1:4896404/1‑151 (MQ=255)
       gCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATg                                                                                                                                        >  2:2039982/1‑151 (MQ=255)
         cATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGc                                                                                                                                      >  1:6074279/1‑151 (MQ=255)
         cATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGc                                                                                                                                      <  1:3642935/151‑1 (MQ=25)
           tACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGt                                                                                                                                    >  1:6050660/1‑151 (MQ=255)
           tACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGACGt                                                                                                                                    >  2:1765477/1‑151 (MQ=255)
            aCCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTg                                                                                                                                   >  1:3205579/1‑151 (MQ=21)
            aCCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTg                                                                                                                                   <  1:4139962/151‑1 (MQ=19)
             ccAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGc                                                                                                                                  <  1:78710/151‑1 (MQ=255)
              cagcagAGGGCGACGGTCAGGTTGCAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCCTCTGCCGGGGTGAGGCAGATCGCGTGcc                                                                                                                                 >  2:1829986/1‑150 (MQ=11)
               agcagAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTGACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGGGCCAGCCTTCTGCCGGGGTGGGGCAGAGCGCGTGccc                                                                                                                                >  2:5113493/1‑149 (MQ=11)
               agcagAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAg                                                                                                                                <  1:1748696/151‑1 (MQ=16)
                 cagAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTg                                                                                                                              <  2:2048085/151‑1 (MQ=12)
                   gagGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGc                                                                                                                            >  2:2250714/1‑151 (MQ=9)
                   gagGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGc                                                                                                                            >  2:2251095/1‑151 (MQ=14)
                    agGGCGCCGGTCAGGTCGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTGACCGGGCCCGCCTGGCCCCCATGGGGCGTAGCGGCCCCAGTCTTCCGCCGGGAGGCGGCAGCTGCCGCGCAGTgccg                                                                                                                           >  2:3745956/1‑148 (MQ=2)
                    agGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGCt                                                                                                                           >  2:2233030/1‑151 (MQ=14)
                    agGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGCt                                                                                                                           >  1:5713699/1‑151 (MQ=14)
                     gggCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGCtt                                                                                                                          >  1:6469902/1‑151 (MQ=14)
                     gggCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGCtt                                                                                                                          <  1:381167/151‑1 (MQ=9)
                      ggCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGCTtg                                                                                                                         <  1:3475672/151‑1 (MQ=2)
                                                                                                                         gtggTGGCGCCAGGCTTCTGGCGGGGTCGGGCGGATGGCGTGCAGTGGCTTGTGGCACCAGGTGGGTTCGCGGGGCGGGGGGAGTTTGTCCGGGCCGTGGTCCATGAAGTGGGGGGGGTATTCCCATACCTCGGCGACATGCTTGACGccc                      <  2:2321991/148‑1 (MQ=1)
                                                                                                                                  cggggtttCGGCCGAGGTGAGGGACATGGGGTGCAGGGGCTTCTGCCGCGCGTTGGGTTCACCGGGGGGCGGGAGTTTGTCCGCGCCGTGGTCTATGAAGTGGCGGAGGTATTCCGATACCTCGGCGACATGCTTGACGCCCGGTTCGGTg             <  2:3756579/145‑1 (MQ=1)
                                                                                                                                             cctgggTGAGGCAGATGGCGTGCAGTGGCTTGTGCCACTAGTTGGGTTCACCGGGCGGGGGGAGTTTGTCCGGGCCGTGGTCCATGAAGTGGCGGAGGTATTCCCATACCTCGGCGACATGCTTGACGCCCGGTTCGGTGGGTTCGTTGCt  <  2:3813827/148‑1 (MQ=0)
                                                                                                                                                    |                                                                                                                                               
GGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGTTTGTCCGGGCCGTGGTCCATGAAGTGGCGGAGGTATTCCCATACCTCGGCGACATGCTTGACGCCCGGTTCGGTGGGTTCGTTGCT  >  CP009974/4536987‑4537278

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