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

CCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT  >  NC_000913/3555904‑3556164
                                                                                                                                 |                                                                                                                                    
cgaccGGGATAATGGTCAGGAGATGGGATGAGGGGTTGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCAc                                                                                                                            <  2:1180898/137‑1 (MQ=255)
ccaccGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCAc                                                                                                                            <  2:643735/139‑1 (MQ=255)
  accGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCCGCGCCATCGGTAGCATCTGTTTGTCAGCAAGCTATTCCCCCTCCGCCGTCGTGCTTGCCAGGTGGC‑GTTTCAccg                                                                                                                          <  2:86140/139‑2 (MQ=255)
  accGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACct                                                                                                                          <  1:222041/139‑1 (MQ=255)
    cgcgATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACctct                                                                                                                        >  2:1474214/1‑139 (MQ=255)
             ggTCAGCAGATGGCATGAGGGATGGGCAACCGAAAATTCCCCCGCGCCGGCCACCGCCATCGTTAGCACAAGTTGGTAAGAAGGTTATTCCCCCACCGCCG‑CGTTGTTTCCAGGTTGCGGTTTCACCCCTTTCAcaacc                                                                                                               >  1:1180898/1‑135 (MQ=255)
                        ggCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAac                                                                                                    <  1:1061000/139‑1 (MQ=255)
                             gAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACctct                                                                                                                        >  1:133300/1‑114 (MQ=255)
                             gAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACctct                                                                                                                        <  2:133300/114‑1 (MQ=255)
                                  aTGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCa                                                                                          <  1:1474214/139‑1 (MQ=255)
                                   tGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCAc                                                                                         >  1:1314505/1‑139 (MQ=255)
                                   tGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCAc                                                                                         >  1:1088335/1‑139 (MQ=255)
                                        gACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAcg                                                                                    <  2:631489/139‑1 (MQ=255)
                                        gACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAcg                                                                                    <  2:350951/139‑1 (MQ=255)
                                         aCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAcgc                                                                                   <  2:152283/139‑1 (MQ=255)
                                                   cccccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCAcc                                                                         <  2:1314505/139‑1 (MQ=255)
                                                    ccccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAaccac                                                                                                 <  2:1058329/114‑1 (MQ=255)
                                                    ccccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAaccac                                                                                                 >  1:1058329/1‑114 (MQ=255)
                                                      ccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAac                                                                                                    <  1:1201608/109‑1 (MQ=255)
                                                      ccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAac                                                                                                    >  2:1201608/1‑109 (MQ=255)
                                                           cccgccAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACaa                                                                 >  2:964033/1‑139 (MQ=255)
                                                                    gcCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGtt                                                        <  1:1402002/139‑1 (MQ=255)
                                                                    gcCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGtt                                                        <  1:793178/139‑1 (MQ=255)
                                                                    gcCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGtt                                                        >  2:1402002/1‑139 (MQ=255)
                                                                     ccATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTc                                                       <  2:1088335/139‑1 (MQ=255)
                                                                      cATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAcgcg                                                                                  <  1:451717/111‑1 (MQ=255)
                                                                      cATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAcgcg                                                                                  >  2:451717/1‑111 (MQ=255)
                                                                               cACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATAttt                                             >  1:1263470/1‑139 (MQ=255)
                                                                                 ccAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGtt                                                        <  2:969279/126‑1 (MQ=255)
                                                                                 ccAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGtt                                                        >  1:969279/1‑126 (MQ=255)
                                                                                    gTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCa                                                                                                                             >  1:80729/1‑54 (MQ=255)
                                                                                    gTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCa                                                                                                                             <  2:80729/54‑1 (MQ=255)
                                                                                    gTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACt                                        >  2:1585701/1‑139 (MQ=255)
                                                                                      tGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCa                                                                                                      >  2:802806/1‑75 (MQ=255)
                                                                                      tGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCa                                                                                                      <  1:802806/75‑1 (MQ=255)
                                                                                            gCAAGATATTCCCCCACCACCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACtt                                >  2:542461/1‑139 (MQ=255)
                                                                                             cAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAaccac                                                                                                 <  1:464648/73‑1 (MQ=255)
                                                                                             cAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAaccac                                                                                                 >  2:464648/1‑73 (MQ=255)
                                                                                                     tCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCgag                        >  2:1407270/1‑138 (MQ=255)
                                                                                                     tCCCCCACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCgag                        <  1:1407270/138‑1 (MQ=255)
                                                                                                          cACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACct                                                                                                                          <  1:530829/35‑1 (MQ=255)
                                                                                                          cACCGCCGCCGTGCTTGCCAGGTGGC‑GTTTCACct                                                                                                                          >  2:530829/1‑35 (MQ=255)
                                                                                                            ccgccgccgTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCa                                                                                                      >  1:992564/1‑53 (MQ=255)
                                                                                                            ccgccgccgTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCa                                                                                                      <  2:992564/53‑1 (MQ=255)
                                                                                                               ccgccgTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATAtttt                                            >  1:977370/1‑108 (MQ=255)
                                                                                                               ccgccgTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATAtttt                                            <  2:977370/108‑1 (MQ=255)
                                                                                                               ccgccgTGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc             <  1:1585701/139‑1 (MQ=255)
                                                                                                                     tGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGccc       >  1:633894/1‑139 (MQ=255)
                                                                                                                     tGCTTGCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTACCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGccc       >  1:1043882/1‑139 (MQ=255)
                                                                                                                          gCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGt  <  2:702339/139‑1 (MQ=255)
                                                                                                                          gCCAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGt  >  2:1389294/1‑139 (MQ=255)
                                                                                                                            cAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc             <  2:891203/126‑1 (MQ=255)
                                                                                                                            cAGGTGGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc             >  1:891203/1‑126 (MQ=255)
                                                                                                                                 |                                                                                                                                    
CCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGC‑GTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT  >  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: ‑