Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,015,846 A→G 9.1% Y183H (TAT→CAT)  PP_1795 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,015,8460AG 9.1% 119.8 / 2.9 44Y183H (TAT→CAT) PP_1795hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (20/20);  new base G (2/2);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01

ACCAGGCCTGTCAGTTCATAGTGTTTGCGCAATGGCCCTACCAGGCTAAAACCTACGCACAGCAGCGGCAGGCCTTGAGCACGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGTGCCCGGGAGGCCGACCCAGCCCTCGCTGGACATCACTTCACCCGCCTTTTTCAGGTGCGCTGTAAGCGTGATGGGCAGGGGC  >  NC_002947/2015698‑2015994
                                                                                                                                                    |                                                                                                                                                    
accaGGCCTGTCAGTTCATAGTGTTTGCGCAATGGCCCTACCAGGCTAAAACCTACGCACAGCAGCGGCAGGCCTTGAGCACGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGtata                                                                                                                                                      <  5:223088/149‑1 (MQ=255)
   aGGCCTGTCAGTTCATAGTGTTTGCGCAATGGCCCTACCAGGCTAAAACCTACGCACAGCAGCGGCAGGCCTTGAGCACGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATagag                                                                                                                                                   <  5:134940/149‑1 (MQ=255)
                  tAGTGTTTGCGCAATGGCCCTACCAGGCTAAAACCTACGCACAGCAGCGGCAGGCCTTGAGCACGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGAc                                                                                                                                    >  5:206406/1‑149 (MQ=255)
                     tgtTTGCGCAATGGCCCTACCAGGCTAAAACCTACGCACAGCAGCGGCAGGCCTTGAGCACGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACggg                                                                                                                                 >  5:122227/1‑149 (MQ=255)
                              aaTGGCCCTACCAGGCTAAAACCTACGCACAGCAGCGGCAGGCCTTGAGCACGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGc                                                                                                                        <  2:287489/149‑1 (MQ=255)
                                  gCCCTACCAGGCTAAAACCTACGCACAGCAGCGGCAGGCCTTGAGCACGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATagag                                                                                                                                                   <  5:114916/118‑1 (MQ=255)
                                  gCCCTACCAGGCTAAAACCTACGCACAGCAGCGGCAGGCCTTGAGCACGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATagag                                                                                                                                                   >  6:114916/1‑118 (MQ=255)
                                    ccTACCAGGCTAAAACCTACGCACAGCAGCGGCAGGCCTTGAGCACGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATc                                                                                                                  >  1:210686/1‑149 (MQ=255)
                                     cTACCAGGCTAAAACCTACGCACAGCAGCGGCAGGCCTTGAGCACGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCAc                                                                                                                                          >  6:305212/1‑124 (MQ=255)
                                        ccAGGCTAAAACCTACGCACAGCAGCGGCAGGCCTTGAGCACGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGAc                                                                                                               <  1:158534/148‑1 (MQ=255)
                                             cTAAAACCTACGCACAGCAGCGGCAGGCCTTGAGCACGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCaaa                                                                                                         <  6:122227/149‑1 (MQ=255)
                                             cTAAAACCTACGCACAGCAGCGGCAGGCCTTGAGCACGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCaaa                                                                                                         <  2:210686/149‑1 (MQ=255)
                                                aaaCCTACGCACAGCAGCGGCAGGGCTTGAGCAGGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCTGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAgg                                                                                                       <  4:210614/148‑1 (MQ=255)
                                                 aaCCTACGCACAGCAGCGGCAGGCCTTGAGCACGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGAt                                                                                                     >  4:237277/1‑149 (MQ=255)
                                                  aCCTACGCACAGCAGCGGCAGGCCTTGAGCACGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGAtt                                                                                                    <  2:279958/149‑1 (MQ=255)
                                                    cTACGCACAGCAGCGGCAGGCCTTGAGCACGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTc                                                                                                   >  4:357225/1‑148 (MQ=255)
                                                     tACGCACAGCAGCGGCAGGCCTTGAGCACGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGa                                                                                                 >  4:174909/1‑149 (MQ=255)
                                                         cacaGCAGCGGCAGGCCTTGAGCACGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGaa                                                                                             <  8:242402/149‑1 (MQ=255)
                                                           cagcagCGGCAGGCCTTGAGCACGAGCGCCGACAAAAGTGCCTGTCGAGACATTGCGGCAGGGGCAGGAGCCCGCCACCATGCAGGTATAGAGCACATCCACCCCGACGCGATTTTCCGCCCAATCGACGGCAAAGGATTCGATGAAgt                                                                                           >  2:38044/1‑149 (MQ=255)
                                                                        ccTTGAGCACGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGTGCCCgg                                                                              <  6:206406/149‑1 (MQ=255)
                                                                            gAGCACGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGTGCCCGGGa                                                                            <  3:237277/147‑1 (MQ=255)
                                                                                  gAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGTGCCCGGGAGGCCGAcc                                                                    >  5:173986/1‑149 (MQ=255)
                                                                                   aGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGTGCCCGGGAGGCCGAcc                                                                    >  6:167097/1‑148 (MQ=255)
                                                                                       ccGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATGGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGt                                                                                    >  1:181937/1‑128 (MQ=255)
                                                                                       ccGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATGGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGt                                                                                    <  2:181937/128‑1 (MQ=255)
                                                                                       cagACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATGGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGt                                                                                    <  2:181940/126‑1 (MQ=255)
                                                                                       cagACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATGGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCCATGAAGTCGTCAGt                                                                                    >  1:181940/3‑128 (MQ=255)
                                                                                         gACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGaa                                                                                             <  1:170799/117‑1 (MQ=255)
                                                                                         gACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGaa                                                                                             >  2:170799/1‑117 (MQ=255)
                                                                                            aaaaGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGTGCCCGGGAGGCCGACCCAGCCCTCGc                                                          <  7:207831/149‑1 (MQ=255)
                                                                                              aaGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGTGCCCGGGAGGCCGACCCAGCCCTCGCTg                                                        <  8:130293/149‑1 (MQ=255)
                                                                                                gTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGa                                                                                                 <  5:329924/106‑1 (MQ=255)
                                                                                                gTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGa                                                                                                 >  6:329924/1‑106 (MQ=255)
                                                                                                         gagaCATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGTGCCCGGGAGGCCGACCCAGCCCTCGCTGGACATCACtt                                              >  8:207672/1‑148 (MQ=255)
                                                                                                         gagaCATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGTGCCCGGGAGGCCGACCCAGCCCTCGCTGGACATCACTTc                                             <  3:357225/149‑1 (MQ=255)
                                                                                                              aTTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGATGTCGTCAGTGCCCGGGAGGCCGACCCAGCCCTCGCTGGACATCACTTCACCCg                                        >  1:295497/1‑149 (MQ=255)
                                                                                                              aTTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGTGCCCGGGAGGCCGACCCAGCCCTCGCTGGACATCACTTCACCCg                                        >  2:204659/1‑149 (MQ=255)
                                                                                                                   ggCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGTGCCCGGGAGGCCGACCCAg                                                                 >  6:292989/1‑119 (MQ=255)
                                                                                                                   ggCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGTGCCCGGGAGGCCGACCCAGCCCTCGCTGGACATCACTTCACCCGCCttt                                   >  7:312208/1‑149 (MQ=255)
                                                                                                                      aGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGGCAGTGCCCGGGAGGCCGACCCAGCCCTCGCTGGACATCACTTCACCCGCCTTTTTc                                <  5:167097/149‑1 (MQ=255)
                                                                                                                              gAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGTGCCCGGGAGGCCGACCCAGCCCTCGCTGGACATCACTTCACCCGCCTTTTTCAGGTgcgc                        <  8:312208/149‑1 (MQ=255)
                                                                                                                               aGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGTGCCCGGGAGGCCGACCCAGCCCTCGCTGGACATCACTTCACCCGCCTTTTTCAGGTGCGCt                       >  3:292331/1‑149 (MQ=255)
                                                                                                                                    gCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGTGCCCGGGAGGCCGACCCAGCCCTCGCTGGACATCACTTCACCCGCCTTTTTCAGGTGCGCTGTAAg                  >  2:314559/1‑149 (MQ=255)
                                                                                                                                     ccaccaTGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCTATGGACGGCAAAGGATTCGATGAAGTCTTCAGTGCCGGGGAGGCCGACCCAGCCCTCGCTGGACATCTCTTCACCCGCCTTTTTCAGGTGCGCTGTAAGc                 <  1:38044/149‑1 (MQ=255)
                                                                                                                                         caTGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGTGCCCGGGAGGCCGACCCAGCCCTCGCTGGACATCACTTCACCCGCCTTTTTCAGGTGCGCTGTAAGCGTGa             >  7:249755/1‑149 (MQ=255)
                                                                                                                                          aTGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGTGCCCGGGAGGCCGACCCAGCCCTCGCTGGACATCACTTCACCCGCCTTTTTCAGGTGCGCTGTAAGCGTGAt            <  3:174909/149‑1 (MQ=255)
                                                                                                                                             cAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGTGCCCGGGAGGCCGACCCAGCCCTCGCTGGACATCACTTCACCCGCCTTTTTCAGGTGCGCTGTAAGCGTGATgg          <  8:249755/148‑1 (MQ=255)
                                                                                                                                             cAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGTGCCCGGGAGGCCGACCCAGCCCTCGCTGGACATCACTTCACCCGCCTTTTTCAGGTGCACTGTAAGCGTGATggg         <  2:295497/149‑1 (MQ=255)
                                                                                                                                                    agagCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGTGCCCGGGAGGCCGACCCAGCCCTCGCTGGACATCACTTCACCCGCCTTTTTCAGGTGCGCTGTAAGCGTGATGGGCAGGGGc  >  1:258339/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACCAGGCCTGTCAGTTCATAGTGTTTGCGCAATGGCCCTACCAGGCTAAAACCTACGCACAGCAGCGGCAGGCCTTGAGCACGAGCGCCGACAAAAGTGCCTGTGGAGACATTGCGGCAGGGGCCGGAGCCGGCCACCATGCAGGTATAGAGCAGATCCACGCCGACGGGCTTTTCAGCCCAATCGACGGCAAAGGATTCGATGAAGTCGTCAGTGCCCGGGAGGCCGACCCAGCCCTCGCTGGACATCACTTCACCCGCCTTTTTCAGGTGCGCTGTAAGCGTGATGGGCAGGGGC  >  NC_002947/2015698‑2015994

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