Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,200,212 T→A 8.9% R276* (AGA→TGA)  PP_3690 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,200,2120TA 8.9% 116.7 / 3.8 45R276* (AGA→TGA) PP_3690hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (18/23);  new base A (2/2);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

AAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGATTACTGAACTC  >  NC_002947/4200066‑4200356
                                                                                                                                                  |                                                                                                                                                
aaGCTCCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTAt                                                                                                                                                <  7:131890/149‑1 (MQ=255)
 aGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTAtt                                                                                                                                               <  3:52357/149‑1 (MQ=255)
   ctacATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGa                                                                                                                                             <  2:74486/146‑1 (MQ=255)
   cTGCATCACGGCGGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACACACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGa                                                                                                                                             <  2:105710/149‑1 (MQ=255)
    tGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATTCATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTAtt                                                                                                                                               >  2:54262/1‑146 (MQ=255)
        tCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTg                                                                                                                                        >  6:103377/1‑149 (MQ=255)
        tCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTg                                                                                                                                        >  2:316873/1‑149 (MQ=255)
             gCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCa                                                                                                                                   <  1:288049/149‑1 (MQ=255)
                           cttgGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAac                                                                                                                      <  5:103377/148‑1 (MQ=255)
                              ggAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAacac                                                                                                                    <  1:207520/147‑1 (MQ=255)
                                            aaTCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTc                                                                                                     >  2:91569/1‑148 (MQ=255)
                                                gTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCaaaaa                                                                                                <  8:171714/149‑1 (MQ=255)
                                                           cGAGAGCTTCAAAGTTCGCCTTCTGAACGCCCTCGATACATAGGCTAAGCTTCAGTGCTGCCGAATCTATGCCTGCCTGCCATGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGt                                                                                     <  4:151853/149‑1 (MQ=255)
                                                                  ttCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTa                                                                               >  2:62921/1‑148 (MQ=255)
                                                                  ttCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCAATAGAGCtt                                                                                                                                         <  4:158290/90‑1 (MQ=255)
                                                                  ttCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCAATAGAGCtt                                                                                                                                         >  3:158290/1‑90 (MQ=255)
                                                                   tCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTaa                                                                              >  1:6327/1‑148 (MQ=255)
                                                                    cAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTa                                                                                    >  8:96936/1‑141 (MQ=255)
                                                                    cAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTa                                                                                    <  7:96936/141‑1 (MQ=255)
                                                                          aCGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGAc                                                                      >  3:269385/1‑149 (MQ=255)
                                                                                tCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAAAAGACCCg                                                                   >  2:43715/1‑146 (MQ=255)
                                                                                         ccTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGaa                                                       <  1:91569/149‑1 (MQ=255)
                                                                                         ccTCGATACATACGCTAACGTTCAGTGCTGCAGAATCTATGCCTGCCTCCCAAGCTCTATTGAGCTTGATCCAGCTCTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTCACCAGTAGCCTAAACACAGAACCAaagaagaaaga                                                        >  6:252701/1‑143 (MQ=255)
                                                                                          cTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGATTGTTCCAGCTGTCATACAACACGCCTCCGTCACGCTCAAAAAATCTGACGACTAGGCTAAAAACAGACCCGAAGAAGTCAGAAg                                                      >  1:243029/1‑149 (MQ=255)
                                                                                               tacatacGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACaa                                                 <  4:269385/149‑1 (MQ=255)
                                                                                                     cGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCAAGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGaa                                            <  5:266281/148‑1 (MQ=255)
                                                                                                      gCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAgg                                          <  1:54262/149‑1 (MQ=255)
                                                                                                        tAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGaa                                                       <  1:203788/134‑1 (MQ=255)
                                                                                                        tAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGaa                                                       >  2:203788/1‑134 (MQ=255)
                                                                                                          aCGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGTGCTTGTTCCAGCTGTCATAAAACACGCCTCCTTCACGCTCAAAAAATCTTACGAGTAGGCTAAAAACAGAGCCGAAGAAGTCAGAAGCACAATTGAAGTGTAt                                      <  1:43715/149‑1 (MQ=255)
                                                                                                          aCGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGt                                        <  5:28122/147‑1 (MQ=255)
                                                                                                            gTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTAtct                                    <  2:6327/149‑1 (MQ=255)
                                                                                                                  tgctgcCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCgaagaa                                                              <  7:80977/117‑1 (MQ=255)
                                                                                                                  tgctgcCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCgaagaa                                                              >  8:80977/1‑117 (MQ=255)
                                                                                                                     tgcCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAgg                                          <  7:2650/134‑1 (MQ=255)
                                                                                                                     tgcCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAgg                                          >  8:2650/1‑134 (MQ=255)
                                                                                                                         gAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGca                                                    >  5:62200/1‑120 (MQ=255)
                                                                                                                         gAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGca                                                    <  6:62200/120‑1 (MQ=255)
                                                                                                                               aTGCCTGCATTCCAATCTCTTTTGAGCTTGTTCCAGCTGTAATAAAACACGGCTCCGTCACGCTAAAAAAATGTGACGTTTAGGCTAAAAACTGACCCGAAGAAGTGAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGtt                  <  5:252701/148‑1 (MQ=255)
                                                                                                                                          ccAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGTAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGATTACTg       >  8:124707/1‑148 (MQ=255)
                                                                                                                                          acAAGCTCAATAGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAg                       <  3:120141/131‑1 (MQ=255)
                                                                                                                                          acAAGCTCAATAGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGCGTATCTCGACAGGCGGAAg                       >  4:120141/2‑132 (MQ=255)
                                                                                                                                           cAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGATTACTGa      >  8:238658/1‑148 (MQ=255)
                                                                                                                                            aaGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGATTACTGAAc    <  1:316873/149‑1 (MQ=255)
                                                                                                                                             aGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGATTACTGAACt   >  4:26417/1‑149 (MQ=255)
                                                                                                                                               ctctATTGAGCTTTTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGATAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCTACAGGCGGAAGGAGTTGACGATTAATGaaata  >  6:213020/1‑145 (MQ=255)
                                                                                                                                               ctctATTGAGCTAGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGCAGTTGACGATTACTGAAATc  >  6:208431/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
AAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGATTACTGAACTC  >  NC_002947/4200066‑4200356

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑