Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,143,125 C→T D415N (GAC→AAC)  rne ← endoribonuclease; RNA‑binding protein;RNA degradosome binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,143,1250CT90.9% 119.6 / 4.0 44D415N (GAC→AAC) rneendoribonuclease; RNA‑binding protein;RNA degradosome binding protein
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (20/20);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

ATTGGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAA  >  NC_000913/1143028‑1143224
                                                                                                 |                                                                                                   
aTTGGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTtac                                                                                                   <  1:460387/100‑1 (MQ=255)
 ttGGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTtacg                                                                                                  <  1:860944/100‑1 (MQ=255)
 ttGGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTtacg                                                                                                  <  2:1367877/100‑1 (MQ=255)
   ggCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTtacgca                                                                                                <  1:1183076/100‑1 (MQ=255)
   ggCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTtacgca                                                                                                >  1:449535/1‑100 (MQ=255)
    gCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTtacgcac                                                                                               >  1:1576422/1‑100 (MQ=255)
      acaGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACgg                                                                                             <  2:822925/100‑1 (MQ=255)
         ggAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGc                                                                                          >  2:724385/1‑100 (MQ=255)
             cAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGt                                                                                      >  2:341986/1‑100 (MQ=255)
              aaTGGCGTGAACTGCCTGGGTGATCTCTTTCAGCGCTTCTTCTTCGATCAGTGGCAGAAGTGAGAGCGACAGCGATTCGTTGTTACGCACGGTGcc                                                                                         <  2:730414/96‑1 (MQ=255)
                 ggCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCa                                                                                  <  1:200612/100‑1 (MQ=255)
                 ggCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCa                                                                                  >  2:1547493/1‑100 (MQ=255)
                  gCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAg                                                                                 <  2:1519707/100‑1 (MQ=255)
                  gCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGTACCAg                                                                                 <  2:1119508/100‑1 (MQ=255)
                      gAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGCCAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGaaca                                                                             >  1:166528/1‑100 (MQ=255)
                           tCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAAcgcg                                                                        <  1:724387/100‑1 (MQ=255)
                            ccTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCgg                                                                       <  1:578292/100‑1 (MQ=255)
                              tGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGAc                                                                     >  2:1092462/1‑100 (MQ=255)
                                   nttCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAAc                                                                 >  2:439312/2‑99 (MQ=255)
                                   gtTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACg                                                                >  2:909039/1‑100 (MQ=255)
                                       tctTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGtgat                                                            >  1:810102/1‑100 (MQ=255)
                                        ctTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGtgatg                                                           <  2:420256/100‑1 (MQ=255)
                                        ctTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGtgatg                                                           >  2:161338/1‑100 (MQ=255)
                                            cAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTg                                                       >  2:574128/1‑100 (MQ=255)
                                               cgcTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGAt                                                    <  2:224311/100‑1 (MQ=255)
                                                  ttcttcttcGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCa                                                 <  1:14944/100‑1 (MQ=255)
                                                     ttcttcGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCAccc                                              >  1:209454/1‑100 (MQ=255)
                                                        ttcGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGt                                           >  2:1077835/1‑100 (MQ=255)
                                                            aTCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATg                                       >  1:1389287/1‑100 (MQ=255)
                                                            aTCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATg                                       >  1:142134/1‑100 (MQ=255)
                                                              cAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATggg                                     <  2:1346268/100‑1 (MQ=255)
                                                               aGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGc                                    <  2:890887/100‑1 (MQ=255)
                                                                gACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGAGCCAGTAGCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCt                                   >  2:280857/1‑100 (MQ=255)
                                                                 aCGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCt                                   >  1:253734/1‑99 (MQ=255)
                                                                 aCGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTc                                  >  1:915471/1‑100 (MQ=255)
                                                                  cGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCa                                 >  1:457105/1‑100 (MQ=255)
                                                                      gAATAGAGCGCGCCAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGgcg                             >  2:775808/1‑100 (MQ=255)
                                                                      gAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGGCTGGATTCACCCAGTGATGGGCTCa                                 >  2:870952/1‑96 (MQ=255)
                                                                        aTAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGgcgc                            <  2:668264/99‑1 (MQ=255)
                                                                         tAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTg                          <  1:598278/100‑1 (MQ=255)
                                                                         tAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTg                          <  2:1384570/100‑1 (MQ=255)
                                                                          agagagCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGa                         >  2:718713/1‑100 (MQ=255)
                                                                           gagagCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGAc                        >  1:95635/1‑100 (MQ=255)
                                                                               gCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGa                    <  2:239918/100‑1 (MQ=255)
                                                                                cGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGan                   <  1:628877/100‑2 (MQ=255)
                                                                                cGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGAc                   <  1:1550481/100‑1 (MQ=255)
                                                                                cGACAGCGATTCGTTGTCACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGAc                   <  2:457104/100‑1 (MQ=255)
                                                                                 gACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACa                  <  2:1479864/100‑1 (MQ=255)
                                                                                    aGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACAttt               <  2:1094188/100‑1 (MQ=255)
                                                                                           cGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCcagcag        <  1:1561945/100‑1 (MQ=255)
                                                                                            gttgttACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAgg       <  1:1519492/100‑1 (MQ=255)
                                                                                                 cacgcacgGTGCCAGTACCAGAACAACGCGGACAAACATGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCGGGCCaaa  >  2:848500/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
ATTGGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAA  >  NC_000913/1143028‑1143224

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: