Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,057,905 A→G 83.2% intergenic (+121/+7) PP_1832 → / ← PP_1833 oxidase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,057,9050AG83.2% 90.7 / 11.0 42intergenic (+121/+7)PP_1832/PP_1833oxidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/5);  new base G (17/17);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTATCGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAGGTCACGACGGATGGCACTGATCGCCCAGGACAGCTGTTCGGCACTGTGCAATTGCCCG  >  NC_002947/2057768‑2058050
                                                                                                                                         |                                                                                                                                                 
gTATGGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTTTTACATTGACTGCGCGATCCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAg                                                                                                                                         <  1:107545/148‑1 (MQ=39)
  aTCGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCgg                                                                                                                                      >  7:241167/1‑149 (MQ=255)
     gacgaGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTa                                                                                                                                    <  3:40924/148‑1 (MQ=255)
           cTGTAATCGCATTGCCGCCCGCTTCGCGCGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCggc                                                                                                                             <  5:192811/149‑1 (MQ=255)
              tAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCggcaag                                                                                                                          <  8:241167/149‑1 (MQ=255)
                 tCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCGACATTGACTGCGCGAGGCCGAAGGGCGGCGCAACCGTCGGGCGTAGGCGCGGGTTCCCCCCCGAAGCGGCCAGTTCAGGCAATTCACCCCGACGCGGTAGGTAgggcaggc                                                                                                                        >  6:228610/1‑148 (MQ=255)
                                             cccGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATc                                                                                               >  6:89548/1‑145 (MQ=255)
                                                    cTACATTGACTGCGCGATGCCTAAGGGCGGCGAAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGAcc                                                                                      <  5:328469/147‑1 (MQ=255)
                                                        tatGACTGCGCGATGCCTAAGGGCGGCTCAACGGGCGGTCGTAGGAGCGGGTTCACCGGCGAAGAGGGCAGTTCAGGCAATACACGCCTCAGCGGTAGGTCGGCAAGTCAAACCGCTGCTGGGTCTGCTTCATCGAGATGACCGGCAAt                                                                                <  2:76471/147‑1 (MQ=255)
                                                              tGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTg                                                                          >  1:17717/1‑149 (MQ=255)
                                                                cgcgATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTg                                                                          <  2:17717/147‑1 (MQ=255)
                                                                cgcgATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTg                                                                          <  7:207209/147‑1 (MQ=255)
                                                                   gATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCctc                                                                     >  3:127875/1‑149 (MQ=255)
                                                                       ccTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCg                                                                 >  5:210853/1‑149 (MQ=255)
                                                                         tAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCgac                                                               >  2:227476/1‑147 (MQ=255)
                                                                           aGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGc                                                               <  1:236361/147‑1 (MQ=255)
                                                                                   gcAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAGGTCacga                                                     >  1:193164/1‑149 (MQ=255)
                                                                                   gcAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAGGTCacga                                                     >  7:155323/1‑149 (MQ=255)
                                                                                        cGTCGGGCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAGGTCACGACGGAt                                                <  7:116818/149‑1 (MQ=255)
                                                                                         gtcggtcgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAAtt                                                                               <  1:236948/117‑1 (MQ=255)
                                                                                         gtcggtcgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAAtt                                                                               >  2:236948/1‑117 (MQ=255)
                                                                                             gtcgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGACCAGGTCACGACGGATGGCAc                                           <  1:227476/149‑1 (MQ=255)
                                                                                              tcgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAGGTCACGACGGATGGCACt                                          >  2:126351/1‑149 (MQ=255)
                                                                                                gTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAAt                                                                                <  7:17020/109‑1 (MQ=255)
                                                                                                gTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAAt                                                                                >  8:17020/1‑109 (MQ=255)
                                                                                                      gCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAg                                                            <  1:134480/123‑1 (MQ=255)
                                                                                                      gCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAg                                                            >  2:134480/1‑123 (MQ=255)
                                                                                                         ggTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAGGTCACGACGGATg                                               >  2:74607/1‑133 (MQ=255)
                                                                                                          gTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCa                                                                             >  7:222074/1‑102 (MQ=255)
                                                                                                          gTTCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCa                                                                             <  8:222074/102‑1 (MQ=255)
                                                                                                           ttCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAGGTCACGACGGATGGCACTGATCGCCCAGGAc                             <  3:235129/149‑1 (MQ=255)
                                                                                                           ttCACCCGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAGGTCACGACGGATGGCACTGATCGCCCAGGAc                             <  8:155323/149‑1 (MQ=255)
                                                                                                               cccGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAGGTCACGCCGGATGGCACTGATCGCCCAGGACAGCt                         >  8:225226/1‑149 (MQ=255)
                                                                                                               cccGCGAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAGGTCACGACGGATGGCACTGATCGCCCAGGACAGCt                         <  6:210853/149‑1 (MQ=255)
                                                                                                                  gcgAAGAGGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAGGTCACGACGGATGGCACGTATCTACCAGGAAAGCTGtt                      >  6:230276/1‑149 (MQ=255)
                                                                                                                       gagGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAGGTCACGACGGATGGCACTGATCGCCCAGGACAg                           >  6:166290/1‑139 (MQ=255)
                                                                                                                       gagGCCAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAGGTCACGACGGATGGCACTGATCGCCCAGGACAg                           <  5:166290/139‑1 (MQ=255)
                                                                                                                            cAGTTCAGGCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAGGTCACGACGGATGGCACTGATCGCCCAGGACAGCTGTTCGGCACtgtg            >  8:75616/1‑149 (MQ=255)
                                                                                                                                    gCAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAGGTCACGACGGATGGCACTGATCGCCCAGGACAGCTGTTCGGCACTGTGCAATTGcc    <  5:116468/149‑1 (MQ=255)
                                                                                                                                     cAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCACTTCACTGGCCTCTTCGGCCAGGTCACGACGGATGGCACTGATCGCCCAGGACAGCTGTTCGGCACTGTGCAATTGccc   <  5:228610/149‑1 (MQ=255)
                                                                                                                                     cAATGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAGGTCACGACGGATGGCACTGATCGCCCAGGACAGCTGTTCGGCACTGTGCAATTGccc   <  4:127875/149‑1 (MQ=255)
                                                                                                                                      aaTGCACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCAAGATGACCGGCAATTCACTGGCCTCTTCGGCCAGGTCACGACGGATGGCACTGATCGCCCAGGACAGCTGTTCGGCACTGTGCAATTGCCCg  >  2:261641/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                 
GTATCGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAGGTCACGACGGATGGCACTGATCGCCCAGGACAGCTGTTCGGCACTGTGCAATTGCCCG  >  NC_002947/2057768‑2058050

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: ‑