Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,192,436:1 +C 100% intergenic (+45/‑159) PP_1937 → / → PP_1938 helicase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,192,4361.C91.2% 102.1 / 3.9 34intergenic (+45/‑159)PP_1937/PP_1938helicase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (14/17);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

GATGTTATTCTGAAGAAAATCCAAGGCGTTGACCTTCGACGGCGTTGGGAAGAGCGAATCGTTGCTTGTATTCTTGGGTATGGGCAAAGCGCCGAGGAGTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGTTATGGCAAGCTGCCGTTCGGGCTGAGCACGTACGGCCCGCTCGTTCCTGGGATAGAGCACGGGACACTCGGGGAAAGCCTTACACTGCCGCTAAGAATGCCAGCTC  >  NC_002947/2192290‑2192579
                                                                                                                                                    |                                                                                                                                                
gATGTTATTCTGAAGAAAATCCAAGGCGTTGACCTTCGACGGCGTTGGGAAGAGCGAATCGTTGCTTGTATTCTTGGGTATGGGCAAAGCGCCGAGGAGTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATgccgcc                                                                                                                                               >  1:159323/1‑149 (MQ=255)
          tGAAGAAAATCCAAGGCGTTGACCTTCGACGGCGTTGGGAAGAGCGAATCGTTGCTTGTATTCTTGGGTATGGGCAAAGCGCCGAGGAGTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGt                                                                                                                                      >  7:6745/1‑149 (MQ=255)
              gaaAATCCAAGGCGTTGACCTTCGACGGCGTTGGGAAGAGCGAATCGTTGCTTGTATTCTTGGGTATGGGCAAAGCGCCGAGGAGTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGttttt                                                                                                                                  >  4:221575/1‑149 (MQ=255)
                aaaTCCAAGGCGTTGACCTTCGACGGCGTTGGGAAGAGCGAATCGTTGCTTGTATTCTTGGGTATGGGCAAAGCGCCGAGGAGTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTg                                                                                                                                <  2:159323/149‑1 (MQ=255)
                            ttGACCTTCGACGGCGTTGGGAAGAGCGAATCGTTGCTTGTATTCTTGGGTATGGGCAAAGCGCCGAGGAGTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTGGTGATACGtttt                                                                                                                    <  8:6745/149‑1 (MQ=255)
                                 cTTCGACGGCGTTGGGAAGAGCGAATCGTTGCTTGTATTCTTGGGTATGGGCAAAGCGCCGAGGAGTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTGGTg                                                                                                                             <  5:211656/135‑1 (MQ=255)
                                 cTTCGACGGCGTTGGGAAGAGCGAATCGTTGCTTGTATTCTTGGGTATGGGCAAAGCGCCGAGGAGTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTGGTg                                                                                                                             >  6:211656/1‑135 (MQ=255)
                                    cGACGGCGTTGGGAAGAGCGAAGCGTTGCTTGTATTCTTGGGTATGGGCAAAGCGCCGAGGAGTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTGGTGATACGTTTTTAgtgt                                                                                                              <  7:244249/147‑1 (MQ=255)
                                     gACGGCGTTGGGAAGAGCGAATCGTTGCTTGTATTCTTGGGTATTGGGCAAAGCGCCGAGGAGTGATGGATCTGTACCAAGGTCAAACCGCTACCTTTCATGTGAATGCCGCCCCCCCGGCGTTTTTTGTGGATACGTTTTttgtttgt                                                                                                            >  1:68666/1‑149 (MQ=255)
                                                 aaGAGCGAATCGTTGCTTGTATTCTTGGGTATGGGCAAAGCGCCGAGGAGTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGTTATAGCAAGCTGc                                                                                               <  8:126369/149‑1 (MQ=255)
                                                                 ttGTATTCTTGGGTATGGGCAAAGCGCCGAGGAGTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTTGTGATACGGTTTTTGTGTGTTATTGCAAGCTGCCGTTGGGGCTGAGCAc                                                                               >  4:12114/1‑149 (MQ=255)
                                                                 ttGTATTCTTGGGTATGGGCAAAGCGCCGAGGAGTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGTTATg                                                                                                        <  7:285343/124‑1 (MQ=255)
                                                                 ttGTATTCTTGGGTATGGGCAAAGCGCCGAGGAGTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGTTATg                                                                                                        >  8:285343/1‑124 (MQ=255)
                                                                  tGTATTCTTGGGTAAGGGCAAAGCGCCGAGGAGTGATGGATCTCTACCCAGGTCCAACCGCTACATGTCATGTGAATGCCGCACCCCCGCCGATTATTGGTTAATCGTTTATAGTGGGATATGGCAAGCTGCCCTAGGGGCtacacac                                                                               >  2:241326/1‑142 (MQ=255)
                                                                       tCTTGGGTATGGGCAAAGCGCCGAGGAGTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGTTATGGc                                                                                                      <  3:291449/119‑1 (MQ=255)
                                                                       tCTTGGGTATGGGCAAAGCGCCGAGGAGTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGTTATGGc                                                                                                      >  4:291449/1‑119 (MQ=255)
                                                                                  gggCAAAGCGCCGAGGAGTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGTGTTTTGGTGATACGTTTTTAGTGTGTTATGGCAAGCTGCCGTTGGGGCTGAGCCCGTACGGCCCGGTCGtt                                                               >  8:268707/1‑149 (MQ=255)
                                                                                      aaaGCGCCGAGGAGTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGGTATGGCAAGCTGCCGTTCGGGCTGAGCACGTACGGCCCGCTCGTTCCTg                                                           <  2:68666/149‑1 (MQ=255)
                                                                                       aaGCGCCGAGGAGTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGTTATGGCAAGCTGCCGTTCGGGCTGAGCACGTACGGCCCGCTCGTTCCTgg                                                          >  1:183493/1‑149 (MQ=255)
                                                                                        aGCGCCGAGGAGTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGTTATGGCAAGCTGCCGTTCGGGCTGAGCACGTACGGCCCGCTCGTTCCTggg                                                         >  5:127217/1‑149 (MQ=255)
                                                                                             cGAGGAGTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGTTATGGCAAGCTGCCGTTCGGGCTGAGCACGTACGGCCCGCTCGTTCCTGGGATaga                                                    <  4:140298/149‑1 (MQ=255)
                                                                                                  agTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTTGTGATACGGTTTTAGTTTGTTATGGCAAGCTGCCGTTCGGGCTGAGCACGGACGGCCCGCTACTTCCTGGGATAGAGccac                                               >  8:287732/1‑146 (MQ=255)
                                                                                                  agTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGTTATGGCAAGCTGCCGTTCGGGCTGAGCACGTACGGCCCGCTCGTTCCTGGGATAGAGCACg                                               <  3:12114/149‑1 (MQ=255)
                                                                                                   gTGAGGGAGCTGTACCCAGGTCCAACCGCGACCTGTAATGTGAATGCCGCCCCCCCGGCGTTTTTTGGTGAGACGTTTTTAGTGTGTTATGGCAAGCTGCCGTTCGGGCTGAGCACGTACGGCCCGCTCGTTCCTGGGATAGAGCACgg                                              <  4:185804/149‑1 (MQ=255)
                                                                                                               tACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGTTATGGCAAGCTGCCGTTCGGGCTGAGCACGTACGGCCCGCTCGTTCCTGGGATAGAGCACGGGACACTCg                                      <  6:127217/145‑1 (MQ=255)
                                                                                                                   cAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGTTATGGCAAGCTGCCGTTCGGGCTGAGCACGTACGGCCCGCTCGTTCCTGGGATAGAGCACGGGACACTCGGGGAAAGc                              <  2:183493/149‑1 (MQ=255)
                                                                                                                   cAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGTTATGGCAAGCTGCCGTTCGGGCTGAGCACGTACGGCCCGCTCGTTCCTGGGATAGAGCACGGGACACTCGGGGAAAGc                              <  3:221575/149‑1 (MQ=255)
                                                                                                                          aaCCGCTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGTTATGGCAAGCTGCCGTTCGGGCTGAGCACGTACGGCCCGCTCGTTCCTGGGATAGAGCACGGGACACTCGGGGAAAGCCTTacac                       >  4:202778/1‑149 (MQ=255)
                                                                                                                              gCTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTTGTGATTCGTTTTTAGTGTGTTATGGCAAGCTGCCGTTCGGGCTGAGCACGGACGGCCCCCTGGTTCCGGGGATAGAGCACGGGACACTCGGGGAAAGCCTTAAACAGcc                   >  6:212840/1‑149 (MQ=255)
                                                                                                                               cTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGTTATGGCAAGCTGCCGTTCGGGCTGAGCACGTACGGCCCGCTCGTTCCTggg                                                         <  6:264492/110‑1 (MQ=255)
                                                                                                                               cTACCTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGTTATGGCAAGCTGCCGTTCGGGCTGAGCACGTACGGCCCGCTCGTTCCTggg                                                         >  5:264492/1‑110 (MQ=255)
                                                                                                                                  ccTGTCATGTGAATGCCGCCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGTTATGGCAAGCTGCCGTTCGGGCTGAGCACGTACGGCCCGCTCGTTCCTGGGATAGAGCACGGGACACTCGGGGAAAGCCTTACACTGCCGCt                <  5:212840/148‑1 (MQ=255)
                                                                                                                                        aTGTGAATGCCGCCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGTTATGGCAAGCTGCCGTTCGGGCTGAGCACGTACGGCCCGCTCGTTCCTGGGATAGAGCACGGGACACTCGGGGAAAGCCTTACACTGCCGCTAAGAATg         >  6:298296/1‑149 (MQ=255)
                                                                                                                                         tgtgAATGCCGCCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGTTATGGCAAGCTGCCGTTCGGGCTGAGCACGTACGGCCCGCTCGTTCCTGGGATAGAGCACGGGACACTCGGGGAAAGCCTTACACTGCCGCTAAGAATGc        <  7:268707/149‑1 (MQ=255)
                                                                                                                                               tGCCGCCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGTTATGGCAAGCTGCCGTTCGGGCTGAGCACGTACGGCCCGCTCGTTCCTGGGATAGAGCACGGGACACTCGGGGAAAGCCTTACACTGCCGCTAAGAATGCCAGCTc  <  5:298296/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GATGTTATTCTGAAGAAAATCCAAGGCGTTGACCTTCGACGGCGTTGGGAAGAGCGAATCGTTGCTTGTATTCTTGGGTATGGGCAAAGCGCCGAGGAGTGATGGATCTGTACCCAGGTCCAACCGCTACCTGTCATGTGAATGCCGCCCCCCGGCGTTTTTTGGTGATACGTTTTTAGTGTGTTATGGCAAGCTGCCGTTCGGGCTGAGCACGTACGGCCCGCTCGTTCCTGGGATAGAGCACGGGACACTCGGGGAAAGCCTTACACTGCCGCTAAGAATGCCAGCTC  >  NC_002947/2192290‑2192579

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