Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,556,033 A→G Y273H (TAC→CAC)  rtcA ← RNA 3'‑terminal phosphate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,556,0330AG98.1% 184.4 / ‑6.3 53Y273H (TAC→CAC) rtcARNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (24/28);  minor base T (1/0);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT  >  NC_000913/3555904‑3556164
                                                                                                                                 |                                                                                                                                    
cgaccGGGATAATGGTCAGGAGATGGGATGAGGGGTTGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCAc                                                                                                                            <  2:1180898/137‑1 (MQ=255)
ccaccGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCAc                                                                                                                            <  2:643735/139‑1 (MQ=255)
  accGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCCGCGCCATCGGTAGCATCTGTTTGTCAGCAAGCTATTCCCCCTCCGCCGTCGTGCTTGCCAGGTGGCGTTTCAccg                                                                                                                          <  2:86140/139‑2 (MQ=255)
  accGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                                          <  1:222041/139‑1 (MQ=255)
    cgcgATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACctct                                                                                                                        >  2:1474214/1‑139 (MQ=255)
             ggTCAGCAGATGGCATGAGGGATGGGCAACCGAAAATTCCCCCGCGCCGGCCACCGCCATCGTTAGCACAAGTTGGTAAGAAGGTTATTCCCCCACCGCCGCGTTGTTTCCAGGTTGCGGTTTCACCCCTTTCAcaacc                                                                                                               >  1:1180898/1‑135 (MQ=255)
                        ggCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac                                                                                                    <  1:1061000/139‑1 (MQ=255)
                             gAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACctct                                                                                                                        >  1:133300/1‑114 (MQ=255)
                             gAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACctct                                                                                                                        <  2:133300/114‑1 (MQ=255)
                                  aTGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCa                                                                                          <  1:1474214/139‑1 (MQ=255)
                                   tGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCAc                                                                                         >  1:1314505/1‑139 (MQ=255)
                                   tGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCAc                                                                                         >  1:1088335/1‑139 (MQ=255)
                                        gACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAcg                                                                                    <  2:631489/139‑1 (MQ=255)
                                        gACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAcg                                                                                    <  2:350951/139‑1 (MQ=255)
                                         aCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAcgc                                                                                   <  2:152283/139‑1 (MQ=255)
                                                   cccccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCAcc                                                                         <  2:1314505/139‑1 (MQ=255)
                                                    ccccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAaccac                                                                                                 <  2:1058329/114‑1 (MQ=255)
                                                    ccccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAaccac                                                                                                 >  1:1058329/1‑114 (MQ=255)
                                                      ccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac                                                                                                    <  1:1201608/109‑1 (MQ=255)
                                                      ccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac                                                                                                    >  2:1201608/1‑109 (MQ=255)
                                                           cccgccAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACaa                                                                 >  2:964033/1‑139 (MQ=255)
                                                                    gcCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGtt                                                        <  1:1402002/139‑1 (MQ=255)
                                                                    gcCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGtt                                                        <  1:793178/139‑1 (MQ=255)
                                                                    gcCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGtt                                                        >  2:1402002/1‑139 (MQ=255)
                                                                     ccATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTc                                                       <  2:1088335/139‑1 (MQ=255)
                                                                      cATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAcgcg                                                                                  <  1:451717/111‑1 (MQ=255)
                                                                      cATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAcgcg                                                                                  >  2:451717/1‑111 (MQ=255)
                                                                               cACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATAttt                                             >  1:1263470/1‑139 (MQ=255)
                                                                                 ccAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGtt                                                        <  2:969279/126‑1 (MQ=255)
                                                                                 ccAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGtt                                                        >  1:969279/1‑126 (MQ=255)
                                                                                    gTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCa                                                                                                                             >  1:80729/1‑54 (MQ=255)
                                                                                    gTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCa                                                                                                                             <  2:80729/54‑1 (MQ=255)
                                                                                    gTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACt                                        >  2:1585701/1‑139 (MQ=255)
                                                                                      tGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCa                                                                                                      >  2:802806/1‑75 (MQ=255)
                                                                                      tGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCa                                                                                                      <  1:802806/75‑1 (MQ=255)
                                                                                            gCAAGATATTCCCCCACCACCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACtt                                >  2:542461/1‑139 (MQ=255)
                                                                                             cAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAaccac                                                                                                 <  1:464648/73‑1 (MQ=255)
                                                                                             cAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAaccac                                                                                                 >  2:464648/1‑73 (MQ=255)
                                                                                                     tCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCgag                        >  2:1407270/1‑138 (MQ=255)
                                                                                                     tCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCgag                        <  1:1407270/138‑1 (MQ=255)
                                                                                                          cACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                                          <  1:530829/35‑1 (MQ=255)
                                                                                                          cACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                                          >  2:530829/1‑35 (MQ=255)
                                                                                                            ccgccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCa                                                                                                      >  1:992564/1‑53 (MQ=255)
                                                                                                            ccgccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCa                                                                                                      <  2:992564/53‑1 (MQ=255)
                                                                                                               ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATAtttt                                            >  1:977370/1‑108 (MQ=255)
                                                                                                               ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATAtttt                                            <  2:977370/108‑1 (MQ=255)
                                                                                                               ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc             <  1:1585701/139‑1 (MQ=255)
                                                                                                                     tGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGccc       >  1:633894/1‑139 (MQ=255)
                                                                                                                     tGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTACCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGccc       >  1:1043882/1‑139 (MQ=255)
                                                                                                                          gCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGt  <  2:702339/139‑1 (MQ=255)
                                                                                                                          gCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGt  >  2:1389294/1‑139 (MQ=255)
                                                                                                                            cAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc             <  2:891203/126‑1 (MQ=255)
                                                                                                                            cAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc             >  1:891203/1‑126 (MQ=255)
                                                                                                                                 |                                                                                                                                    
CCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT  >  NC_000913/3555904‑3556164

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: