Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,966,720 C→T R102C (CGC→TGC)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,966,7200CT95.2% 203.4 / 0.9 62R102C (CGC→TGC) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (33/26);  total (35/27)
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.85e-01

CAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTG  >  NC_000913/3966623‑3966816
                                                                                                 |                                                                                                
cAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgc                                                                                                >  2:347745/1‑100 (MQ=255)
  gCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgccc                                                                                              >  2:3071/1‑100 (MQ=255)
   cTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgcccg                                                                                             <  2:1269784/100‑1 (MQ=255)
   cTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgcccg                                                                                             >  1:1616214/1‑100 (MQ=255)
      cTACCTCGCCGGTACTGATGACATCTACGTTTCCCCTACCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgcccgcc                                                                                           >  1:1371479/1‑99 (MQ=255)
         ccTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGaaa                                                                                       >  2:1142862/1‑100 (MQ=255)
          cTCGCCGGTCCTGATGACATCTCCGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAaag                                                                                      >  1:1293781/1‑100 (MQ=255)
          cTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAaag                                                                                      >  1:886656/1‑100 (MQ=255)
          cTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAaag                                                                                      <  2:1738127/100‑1 (MQ=255)
                   ccTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAAc                                                                             <  2:1802462/100‑1 (MQ=255)
                       atgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTa                                                                         >  2:1645828/1‑100 (MQ=255)
                          aCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTAttt                                                                      <  1:277460/100‑1 (MQ=255)
                          aCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTAttt                                                                      <  1:2149654/100‑1 (MQ=255)
                           cATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTAtttt                                                                     >  2:1248574/1‑100 (MQ=255)
                            aTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTg                                                                    <  2:1919360/100‑1 (MQ=255)
                             tCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgc                                                                   <  1:1046034/100‑1 (MQ=255)
                             tCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgc                                                                   <  1:2059536/100‑1 (MQ=255)
                               tACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgcgc                                                                 >  2:1251130/1‑100 (MQ=255)
                                 cGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctg                                                               <  2:1999054/100‑1 (MQ=255)
                                        ccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAg                                                        >  1:1917747/1‑100 (MQ=255)
                                         cTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGt                                                       <  2:308767/100‑1 (MQ=255)
                                            gCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTaa                                                    >  1:2213519/1‑100 (MQ=255)
                                            gCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTaa                                                    >  2:208142/1‑100 (MQ=255)
                                            gCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTaa                                                    >  1:490327/1‑100 (MQ=255)
                                              cAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg                                                  >  2:1881437/1‑100 (MQ=255)
                                              cAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg                                                  <  2:1558830/100‑1 (MQ=255)
                                                aaTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGaa                                                >  1:1118452/1‑100 (MQ=255)
                                                 aTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAg                                               >  1:1877514/1‑100 (MQ=255)
                                                        gTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACtt                                        <  1:2127112/100‑1 (MQ=255)
                                                        gTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACtt                                        <  2:1321980/100‑1 (MQ=255)
                                                          ttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCg                                      >  2:401337/1‑100 (MQ=255)
                                                          ttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCg                                      >  1:1092689/1‑100 (MQ=255)
                                                          ttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCg                                      <  1:956878/100‑1 (MQ=255)
                                                            cAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                    >  2:708085/1‑100 (MQ=255)
                                                            cAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                    <  2:1683838/100‑1 (MQ=255)
                                                                cTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAc                                <  1:347744/100‑1 (MQ=255)
                                                                 tCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAcc                               >  1:894276/1‑100 (MQ=255)
                                                                 tCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAcc                               <  2:886657/100‑1 (MQ=255)
                                                                 tCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAcc                               <  1:141055/100‑1 (MQ=255)
                                                                 tCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAcc                               <  1:1804884/100‑1 (MQ=255)
                                                                 tCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTAAACGATATTTTGCGCTGATGAACGTTACCGAAATTAACTTCGACAAccc                               >  2:1858479/1‑100 (MQ=255)
                                                                         ggTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACg                       >  2:459559/1‑100 (MQ=255)
                                                                         ggTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACg                       >  1:634595/1‑100 (MQ=255)
                                                                         ggTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCAGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACg                       >  2:2215128/1‑100 (MQ=255)
                                                                              tACCATCTCTGGTAAGCTTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCg                   >  1:1805423/1‑99 (MQ=255)
                                                                              tACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                  >  2:719749/1‑100 (MQ=255)
                                                                               aCCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCAGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGca                 <  1:2215114/100‑1 (MQ=255)
                                                                                 cATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCAcaac               >  2:79420/1‑100 (MQ=255)
                                                                                      ctGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt          >  1:352125/1‑100 (MQ=255)
                                                                                      ctGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt          >  2:1376497/1‑100 (MQ=255)
                                                                                      ctGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt          <  1:1142861/100‑1 (MQ=255)
                                                                                        ggaaaGCTTTGCCCGCCGCAAGCACCTGAACGCCATTATGCGCTGCTGAAAGTTACCGACGTTAACTTCGACAAACCTGAAAACGCCCGCAACACAATcc        <  2:657978/97‑1 (MQ=255)
                                                                                        ggTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc        <  1:1442673/100‑1 (MQ=255)
                                                                                        ggTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc        <  1:690224/100‑1 (MQ=255)
                                                                                        ggTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc        >  2:347469/1‑100 (MQ=255)
                                                                                        ggTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc        >  2:2210971/1‑100 (MQ=255)
                                                                                        ggTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc        >  1:1643330/1‑100 (MQ=255)
                                                                                         gTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCct       >  1:1120909/1‑100 (MQ=255)
                                                                                         gTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCct       >  2:667980/1‑100 (MQ=255)
                                                                                         gTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCct       >  2:2038657/1‑100 (MQ=255)
                                                                                          tAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctc      <  2:276814/100‑1 (MQ=255)
                                                                                          tAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctc      >  2:846504/1‑100 (MQ=255)
                                                                                          tAAGATTCGCCCGCCGAAAGAAGGTGAACGCTCTGTTGCGCTTCTGAAAGTTAACGAAGTTAACTCCCACAAACCTGAAAACACCCGCAACAAAAccc        >  1:591326/1‑98 (MQ=255)
                                                                                           aaGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctct     <  2:490328/100‑1 (MQ=255)
                                                                                           aaGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctct     <  2:1371481/100‑1 (MQ=255)
                                                                                            aGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCtt    <  2:1293783/100‑1 (MQ=255)
                                                                                              aTTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg  >  1:41743/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTG  >  NC_000913/3966623‑3966816

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