Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,165,021 A→G 10.8% V187A (GTA→GCA) ‡ PP_1920 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,165,0210AG10.8% 97.3 / 3.3 37V187A (GTA→GCA) ‡PP_1920hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/12);  new base G (2/2);  total (23/14)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

TCGCAGTGGTAAAGACTTCATAGAACGTCGTGAAGTCTTTGTTGTACTCGATCATATCTTCGAAGGAGTTGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGCCGGTATCGTGAACTGATGCTTCCTCAGTAAACGAGCCATGTCCAAACCAATGATGCAAAAGCGGGTCGGTGTAGTTTTGCTCGCAATCACTCCAGCGCCAGTTGAC  >  NC_002947/2164875‑2165160
                                                                                                                                                  |                                                                                                                                           
tCGCAGTGGTAAAGACTTCATAGAACGTCGTGAAGTCTTTGTTGTACTCGATCATATCTTCGAAGGAGTTGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAg                                                                                                                                          >  1:530624/1‑150 (MQ=255)
               cTTCATAGAACGTCGTGAAGTCTTTGTTGTACTCGATCATATCTTCGAAGGAGTTGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTagag                                                                                                                           >  1:238811/1‑150 (MQ=255)
                  cATAGAACGTCGTGAAGTCTTTGTTGTACTCGATCATATCTTCGAAGGAGTTGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGAtt                                                                                                                        >  2:393277/1‑150 (MQ=255)
                       aaCGTCGTGAAGTCTTTGTTGTACTCGATCATATCTTCGAAGGAGTTGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCaa                                                                                                                   <  2:219039/150‑1 (MQ=255)
                          gtcgtGAAGTCTTTGTTGTACTCGATCATATCTTCGAAGGAGTTGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGAtt                                                                                                                        >  2:229535/1‑142 (MQ=255)
                          gtcgtGAAGTCTTTGTTGTACTCGATCATATCTTCGAAGGAGTTGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGAtt                                                                                                                        <  1:229535/142‑1 (MQ=255)
                             gtGAAGTCTTTGTTGTACTCGATCATATCTTCGAAGGAGTTGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAg                                                                                                             >  1:540881/1‑150 (MQ=255)
                              tGAAGTCTTTGTTGTACTCGATCATATCTTCGAAGGAGTTGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGc                                                                                                            >  2:301547/1‑150 (MQ=255)
                               gAAGTCTTTGTTGTACTCGATCATATCTTCGAAGGAGTTGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGca                                                                                                           >  2:536802/1‑149 (MQ=255)
                                aaGTCTTTGTTGTACTCGATCATATCTTCGAAGGAGTTGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGCCg                                                                                                          <  2:530624/150‑1 (MQ=255)
                                           gtACTCGATCATATCTTCGAAGGAGTTGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGACAAGGAAAGCCGGTATCGTGAAc                                                                                               >  1:480152/1‑150 (MQ=255)
                                               tCGATCATATCTTCGAAGGAGTTGGTCCACCCGCCACCAATAAGATCAGTACTATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGCCGGTATCGTGAACTGAt                                                                                           >  2:304793/1‑150 (MQ=255)
                                                 gATCATATCTTCGAAGGAGTTGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAg                                                                                                                                          >  1:363325/1‑101 (MQ=255)
                                                 gATCATATCTTCGAAGGAGTTGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAg                                                                                                                                          <  2:363325/101‑1 (MQ=255)
                                                   tCATATCTTCGAAGGAGTTGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGCCGGTATCGTGAACTGATGCtt                                                                                       >  2:544578/1‑150 (MQ=255)
                                                         cTTCGAAGGAGTTGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGCCGGTATCGTGAACTGATGCTTCCTCAg                                                                                 <  1:393277/150‑1 (MQ=255)
                                                          ttCGAAGGAGTTGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTGTAGCCATgt                                                                                                                                    >  1:6838/1‑97 (MQ=255)
                                                          ttCGAAGGAGTTGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTGTAGCCATgt                                                                                                                                    <  2:6838/98‑2 (MQ=255)
                                                          ttCGAAGGAGTTGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGCCGGTATCGTGAACTGATGCTTCCTCAGt                                                                                >  1:413101/1‑150 (MQ=255)
                                                                    ttGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACa                                                                                                                                           <  1:571056/81‑1 (MQ=255)
                                                                    ttGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACa                                                                                                                                           >  2:571056/1‑81 (MQ=255)
                                                                          cACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTGTAGCCATgttg                                                                                                                                  <  2:136622/84‑4 (MQ=255)
                                                                          cACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTGTAGCCATggcg                                                                                                                                  >  1:136622/1‑81 (MQ=255)
                                                                             ccGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGCCGGTATCGTGAACTGATGCTTCCTCAGTAAACGAGCCATGCCCAAAc                                                             >  1:104955/1‑150 (MQ=255)
                                                                                ccaccaATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGCCGGTATCGTGAACTGATGCTTCCTCAGTAAACGAGCCATGTCCAAACCaa                                                          <  1:544578/150‑1 (MQ=255)
                                                                                  accaATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGCCGGTATCGTGAACTGATGCTTCCTCAGTAAACGAGCCATGTCCAAACCAatg                                                        >  1:369343/1‑150 (MQ=255)
                                                                                    caATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGCCGGTATCGTGAACTGATGCTTCCTCAGTAAACGAGCCATGTCCAAACCAatgat                                                      <  1:301547/150‑1 (MQ=255)
                                                                                     aaTAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGCCg                                                                                                          <  1:472250/97‑1 (MQ=255)
                                                                                     aaTAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGCCg                                                                                                          >  2:472250/1‑97 (MQ=255)
                                                                                          gATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGCCGGTATCGTGAACTGAt                                                                                           >  2:280735/1‑107 (MQ=255)
                                                                                          gATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGCCGGTATCGTGAACTGAt                                                                                           <  1:280735/107‑1 (MQ=255)
                                                                                                 aaTATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGCCGGTATCGTGAACTGATGCTTCCTCAGTAAACGAGCCATGTCCAAACCAATGATGCAAAAGCGGGTc                                         <  2:369343/150‑1 (MQ=255)
                                                                                                                   cGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGCCGGTATCGTGAACTGATGCTTCCTCAGTAAACGAGCCATGTCCAAACCAATGATGCAAAAGCGGGTCGGTGTAGTTTTGCTCGCa                       >  1:470308/1‑150 (MQ=255)
                                                                                                                              gCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGCCGGTATCGTGAACTGATGCTTCCTCAGTAAACGAGCCATGTCCAAACCAATGATGCAAAAGCGGGTCGGTGTAGTTTTGCTCGCAATCACTCCAgc            >  1:323858/1‑150 (MQ=255)
                                                                                                                                     ctcCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGCAGGTATCGTGAACTGAAGCTTCCTCAGTATACGATCCATGTCCAAACCAATGATGCAAAAGCGGGTCGGTGTAGTTTTGCTCGCAATCACTCCAGCGCCAGtt     <  1:536802/150‑1 (MQ=255)
                                                                                                                                       ccAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGCCGGTATCGTGAACTGATGCTTCCTCAGTAAACGAGCCATGTCCAAACCAATGATGCAAAAGCGGGTCGGTGTAGTTTTGCTCGCAATCACTCCAGCGCCAGTTGa   >  1:291671/1‑150 (MQ=255)
                                                                                                                                        cAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGCCGGTATCGTGAACTGATGCTTCCTCAGTAAACGAGCCATGTCCAAACCAATGATGCAAAAGCGGGTCGGTGTAGTTTTGCTCGCAATCACTCCAGCGCCAGTTGAc  >  1:606049/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
TCGCAGTGGTAAAGACTTCATAGAACGTCGTGAAGTCTTTGTTGTACTCGATCATATCTTCGAAGGAGTTGGTCCACCCGCCACCAATAAGATCAGTAATATCGAGGGTGCACTCCGCATCAACGGGCAGAATCTCCAACATGGCTACAGCCATGCTTTCTAGAGATTGGCAAGGAAAGCCGGTATCGTGAACTGATGCTTCCTCAGTAAACGAGCCATGTCCAAACCAATGATGCAAAAGCGGGTCGGTGTAGTTTTGCTCGCAATCACTCCAGCGCCAGTTGAC  >  NC_002947/2164875‑2165160

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