Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,914,888 T→G Q216P (CAG→CCG)  proQ ← RNA chaperone, putative ProP translation regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,914,8880TG96.6% 204.7 / ‑1.8 59Q216P (CAG→CCG) proQRNA chaperone, putative ProP translation regulator
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (28/29);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

AAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCT  >  NC_000913/1914796‑1914983
                                                                                            |                                                                                               
aaGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAACGACATACCCGAATTCAGCGggacgcg                                                                                          >  1:174153/1‑100 (MQ=255)
 aGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGggacgcgg                                                                                         <  2:1395462/100‑1 (MQ=255)
  gCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGggacgcgga                                                                                        >  2:1154489/1‑100 (MQ=255)
   ccTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGggacgcggac                                                                                       >  2:368621/1‑100 (MQ=255)
   ccTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGggacgcggac                                                                                       <  1:1491828/100‑1 (MQ=255)
        aaaaCATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGt                                                                                  >  1:704320/1‑100 (MQ=255)
        aaaaCATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGt                                                                                  <  2:1189746/100‑1 (MQ=255)
        aaaaCATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGt                                                                                  <  1:389418/100‑1 (MQ=255)
             aTGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTg                                                                             <  2:687699/100‑1 (MQ=255)
             aTGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTg                                                                             <  2:1071646/100‑1 (MQ=255)
              tGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTgg                                                                            >  1:950789/1‑100 (MQ=255)
                   aggCCTGGCCCGGCGTCCGTTTCAGCACACCAGGTGTTATGCGCGCACGCGCACAGACATACCCGAAATCAGCGGGACGCGGACGCCGTCTTTGGTGagt                                                                       >  1:987716/3‑98 (MQ=255)
                     gCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTc                                                                     >  2:1401442/1‑100 (MQ=255)
                      ccTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCt                                                                    <  2:499398/100‑1 (MQ=255)
                       cTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTa                                                                   >  2:517424/1‑100 (MQ=255)
                            ccGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACg                                                              <  2:236451/100‑1 (MQ=255)
                                   ccGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGcat                                                       <  1:1536016/100‑1 (MQ=255)
                                     gTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTgg                                                          <  1:799790/95‑1 (MQ=255)
                                     gTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGcatcc                                                     >  1:1589936/1‑100 (MQ=255)
                                       cgcAGAACACCTGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGcatccat                                                   >  1:1305962/3‑100 (MQ=255)
                                        gcAGAACACCTGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGcatccatc                                                  >  2:146772/2‑100 (MQ=255)
                                         cAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATcg                                                 >  1:1321784/1‑100 (MQ=255)
                                         cAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATcg                                                 >  1:450768/1‑100 (MQ=255)
                                          aGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATcgc                                                >  2:44308/1‑100 (MQ=255)
                                           gAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATcgcg                                               >  1:1530171/1‑100 (MQ=255)
                                               acCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGtttt                                           <  2:1094930/100‑1 (MQ=255)
                                               acCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGtttt                                           >  2:1475741/1‑100 (MQ=255)
                                                tctggTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTg                                          <  2:594066/97‑1 (MQ=255)
                                                ccAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTg                                          <  1:154781/100‑1 (MQ=255)
                                                 cAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGa                                         >  1:363899/1‑100 (MQ=255)
                                                   ggTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGaga                                       >  1:594062/1‑98 (MQ=255)
                                                         cTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCttt                                 >  2:1299480/1‑100 (MQ=255)
                                                         cTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCttt                                 >  1:266940/1‑100 (MQ=255)
                                                          tGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGTCATCCATCGCGTTTTGACCCGCTTTc                                <  2:974355/100‑1 (MQ=255)
                                                          tGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTc                                <  1:490184/100‑1 (MQ=255)
                                                               gcACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCtt                           >  2:212016/1‑100 (MQ=255)
                                                                cacaATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTc                          <  1:1340835/100‑1 (MQ=255)
                                                                cacaATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTc                          <  1:623517/100‑1 (MQ=255)
                                                                  caATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAg                        >  1:1160396/1‑100 (MQ=255)
                                                                   aaTCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAgg                       <  1:262426/100‑1 (MQ=255)
                                                                     tCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGc                     <  1:264025/100‑1 (MQ=255)
                                                                     tCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGc                     >  1:1509796/1‑100 (MQ=255)
                                                                      cAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGcn                    <  1:797619/100‑2 (MQ=255)
                                                                      cAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCt                    <  2:1202266/100‑1 (MQ=255)
                                                                       aaaGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCtt                   >  1:1037932/1‑100 (MQ=255)
                                                                        aaGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTg                  <  1:1154475/100‑1 (MQ=255)
                                                                          gACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTc                <  1:517420/100‑1 (MQ=255)
                                                                           aCATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTcc               <  2:25364/100‑1 (MQ=255)
                                                                            cATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCg              <  2:878104/100‑1 (MQ=255)
                                                                             aTACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGa             >  2:927332/1‑100 (MQ=255)
                                                                              tACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGAc            <  1:1299466/100‑1 (MQ=255)
                                                                              tACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGAc            >  1:804512/1‑100 (MQ=255)
                                                                               aCCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACa           >  2:1095792/1‑100 (MQ=255)
                                                                                cccGAATTAAGCGGGACGCGGCCGCCGTCTTTGGTGATTTCTAATACGGTGGCGTCCATCGCGTTTTGACCTGCTTTCACCTTCAGGGCTTGTCCGACAg          >  1:1109759/1‑100 (MQ=255)
                                                                                   gAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCa       <  2:266943/100‑1 (MQ=255)
                                                                                   gAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCa       <  1:1061241/100‑1 (MQ=255)
                                                                                     aTTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCaga     <  2:32052/100‑1 (MQ=255)
                                                                                        cAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCt  <  1:1090989/100‑1 (MQ=255)
                                                                                        cAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCt  >  2:1012661/1‑100 (MQ=255)
                                                                                            |                                                                                               
AAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCT  >  NC_000913/1914796‑1914983

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