Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,159,126 A→G 16.9% intergenic (+30/‑179) PP_3661 → / → PP_3662 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,159,1260AG16.9% 63.4 / 4.3 29intergenic (+30/‑179)PP_3661/PP_3662membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (10/14);  new base G (2/3);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCT  >  NC_002947/4158981‑4159253
                                                                                                                                                 |                                                                                                                               
tCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAAcc                                                                                                                             <  2:163880/150‑1 (MQ=255)
  aGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCg                                                                                                                           >  1:405229/1‑150 (MQ=255)
   gCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGGCTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGc                                                                                                                          <  1:459788/150‑1 (MQ=255)
          tCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACa                                                                                                                   <  2:19294/150‑1 (MQ=255)
            gCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAgg                                                                                                                 >  2:239984/1‑150 (MQ=255)
                 gCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGATCCTACAggtgggt                                                                                                            <  2:85685/150‑1 (MQ=255)
                                cTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACg                                                                                             <  1:464937/150‑1 (MQ=255)
                                              gATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGccc                                                                                    >  2:41378/1‑145 (MQ=255)
                                              gATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGccc                                                                                    <  1:41378/145‑1 (MQ=255)
                                                       ccTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa                                                                            <  2:111298/144‑1 (MQ=255)
                                                       ccTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa                                                                            >  1:111298/1‑144 (MQ=255)
                                                             gCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCAccc                                                                >  2:330951/1‑150 (MQ=255)
                                                               tGGTTGGCGGGGGGCTGGAGATTAACATCGGGGTGCAGAAGCTGTTTTGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCcgc                                                              <  2:387123/150‑1 (MQ=38)
                                                                                gCGATTCCCCTGGGGGTGCAGAAGCTGTTGTGCTGATAGGAGAGCGAGTCGGCTGTTTCGCGGGTGAACCCGCTCCTACAGGTGGGTGCTGCCCTCAAGACGGCAAGGTCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGt                                             <  2:163103/150‑1 (MQ=255)
                                                                                   aTTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCggg                                                                       >  2:161134/1‑121 (MQ=255)
                                                                                   aTTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCggg                                                                       <  1:161134/121‑1 (MQ=255)
                                                                                    ttACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCa                                         <  2:508081/150‑1 (MQ=255)
                                                                                                    gAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGccc                         >  1:158675/1‑150 (MQ=255)
                                                                                                            tGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGtt                 <  2:391819/150‑1 (MQ=255)
                                                                                                              gtcTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCAGCTCATACCGGCGGGTGCTGCCCTGAAGACGTCAAGGCCACTGTCGAAACGGATTCATTCTTTAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTga               >  1:391819/3‑149 (MQ=255)
                                                                                                                     ttgagaCCGAGCTGCCTGCTTCGCGGGTGAACCCGCTCCTACAGGTGGGCGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa                                                                            >  2:191310/3‑82 (MQ=17)
                                                                                                                     ttgagaCCGAGCTGCCTGCTTCGCGGGTGAACCCGCTCCTACAGGTGGGCGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa                                                                            >  2:191311/3‑82 (MQ=17)
                                                                                                                     ttgagaCCGAGCTGCCTGCTTCGCGGGTGAACCCGCTCCTACAGGTGGGCGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa                                                                            <  1:191311/80‑1 (MQ=18)
                                                                                                                     ttgagaCCGAGCTGCCTGCTTCGCGGGTGAACCCGCTCCTACAGGTGGGCGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa                                                                            <  1:191310/80‑1 (MQ=17)
                                                                                                                           ccGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTtct  <  2:405229/150‑1 (MQ=255)
                                                                                                                              aGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTg                <  1:76112/133‑1 (MQ=255)
                                                                                                                              aGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTg                >  2:76112/1‑133 (MQ=255)
                                                                                                                                        ttCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGt                                             >  1:502473/1‑94 (MQ=255)
                                                                                                                                        ttCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGt                                             <  2:502473/94‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                               
TCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCT  >  NC_002947/4158981‑4159253

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