Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,200,215 T→A 9.1% N275Y (AAT→TAT)  PP_3690 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,200,2150TA 9.1% 121.1 / 3.2 44N275Y (AAT→TAT) PP_3690hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (19/21);  new base A (2/2);  total (21/23)
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.95e-01

AGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGATTACTGAACTC  >  NC_002947/4200067‑4200356
                                                                                                                                                    |                                                                                                                                             
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                                                                                                                                        >  2:316873/1‑149 (MQ=255)
       tCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTg                                                                                                                                        >  6:103377/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                                                                                    <  7:96936/141‑1 (MQ=255)
                                                                   cAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTa                                                                                    >  8:96936/1‑141 (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                                                       >  2:203788/1‑134 (MQ=255)
                                                                                                       tAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGaa                                                       <  1:203788/134‑1 (MQ=255)
                                                                                                         aCGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGTGCTTGTTCCAGCTGTCATAAAACACGCCTCCTTCACGCTCAAAAAATCTTACGAGTAGGCTAAAAACAGAGCCGAAGAAGTCAGAAGCACAATTGAAGTGTAt                                      <  1:43715/149‑1 (MQ=255)
                                                                                                         aCGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGt                                        <  5:28122/147‑1 (MQ=255)
                                                                                                           gTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTAtct                                    <  2:6327/149‑1 (MQ=255)
                                                                                                                 tgctgcCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCgaagaa                                                              >  8:80977/1‑117 (MQ=255)
                                                                                                                 tgctgcCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCgaagaa                                                              <  7:80977/117‑1 (MQ=255)
                                                                                                                    tgcCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAgg                                          >  8:2650/1‑134 (MQ=255)
                                                                                                                    tgcCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAgg                                          <  7:2650/134‑1 (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)
                                                                                                                                                    tGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGATTACt        <  1:46863/136‑1 (MQ=255)
                                                                                                                                                    tGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGATTACt        >  2:46863/1‑136 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
AGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGATTACTGAACTC  >  NC_002947/4200067‑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: ‑