Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,610,980 G→A 10.8% P902S (CCC→TCC)  mcoA ← Mn(II) copper oxidase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,610,9800GA10.8% 95.7 / 3.1 37P902S (CCC→TCC) mcoAMn(II) copper oxidase A
Reads supporting (aligned to +/- strand):  ref base G (18/15);  new base A (2/2);  total (20/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 = 7.25e-01

TGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGCGAACTCGCGGAAGCGGATGGCCATTTCCACCTCTTCGGACGAA  >  NC_002947/3610837‑3611122
                                                                                                                                               |                                                                                                                                              
tGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTa                                                                                                                                           >  4:274783/1‑149 (MQ=255)
    gTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTg                                                                                                                                       <  3:274783/149‑1 (MQ=255)
     tCGGTGCGGCACTGTCGTTGGCCACGACAGGTGGGTTGTTGGCCGGGGCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTgt                                                                                                                                      <  7:164145/149‑2 (MQ=255)
          gaggCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCTTCGTCGTCATCGGTTATGGTGCGGAACGTTGGCAGGCCACCCGAGGCCATGTAAAGCACGACATCACATCCGGGCAGCGGGGTGGGCATTACCTCGAACTg                                                                                                                                 >  8:294120/3‑149 (MQ=255)
            ggCACTGTCGTTGGCCACGACAGGTGGTTTGTTAGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGcc                                                                                                                               >  4:129532/1‑149 (MQ=255)
                   gcGTTGGCCACGTCAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGGTTTGGTGCGGTACGTTGGCAGGCCAACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCggg                                                                                                                          <  7:294120/146‑1 (MQ=255)
                          ccACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAgg                                                                                                                 <  6:21937/149‑1 (MQ=255)
                              gACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTccc                                                                                                             <  7:196109/149‑1 (MQ=255)
                                       tttGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGTCAGGACTACCGAGGCCATGTAACGCACGCCATCACATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTAAAGGTCCCAGCgcag                                                                                                      >  5:285278/1‑147 (MQ=255)
                                               ccGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGaa                                                                                            <  3:129532/149‑1 (MQ=255)
                                               ccGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCATACCGAGGACATGTAACGCACGACATCCCATCCTGGCAGCGTGTTGGGCATTACCTGGAACTGCACCGGGTGTTCAAGGTCCCAGCGCAGAAGCATCGaa                                                                                            <  6:285278/149‑1 (MQ=255)
                                                                tcatcgGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTCCAAGGTCCCAGCGCAGCAGCATCGAAGAGccc                                                                                      <  5:15077/138‑1 (MQ=255)
                                                                tcatcgGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTCCAAGGTCCCAGCGCAGCAGCATCGAAGAGccc                                                                                      >  6:15077/1‑138 (MQ=255)
                                                                  atcgGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGaa                                                                                                                                    <  8:223313/90‑1 (MQ=255)
                                                                  atcgGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCgcagca                                                                                                    >  5:132272/1‑122 (MQ=255)
                                                                  atcgGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCgcagca                                                                                                    <  6:132272/122‑1 (MQ=255)
                                                                  atcgGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGgtgt                                                                         >  5:323868/1‑149 (MQ=255)
                                                                       ttttGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTgg                                                                    >  7:52445/1‑149 (MQ=255)
                                                                              gCGGAACTTTGGCAGGCCTACCGAGGCCATGAAACGCACCCCATACCATCCTGACACCCGCCTCGGCATTACCTGGAACTGCCCCCGCTGTACAAGGTCCCAGCGCAGCAGCATCCAAGAGACCTCGTCCTCCCTGTTGAGGCAGTGCt                                                             >  2:118548/1‑149 (MQ=255)
                                                                                    cGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGACCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCATCGTGCTGGGTGTTGTGGCAGTGCTGCATGt                                                       <  6:323868/149‑1 (MQ=255)
                                                                                      ttGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAg                                                     >  6:28439/1‑149 (MQ=255)
                                                                                               cTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGc                                             >  5:136022/1‑148 (MQ=255)
                                                                                                    gAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGCGAACTc                                       <  8:52445/149‑1 (MQ=255)
                                                                                                       gCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGCGAACTcgcg                                    >  5:18915/1‑149 (MQ=255)
                                                                                                                      ccatcccatcCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTcc                                                           <  2:162031/111‑1 (MQ=255)
                                                                                                                      ccatcccatcCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTcc                                                           >  1:162031/1‑111 (MQ=255)
                                                                                                                        atcccatccTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAg                                                     <  6:12068/115‑1 (MQ=255)
                                                                                                                        atcccatccTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAg                                                     >  5:12068/1‑115 (MQ=255)
                                                                                                                        atcccatccTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGCGAACTCGCGGAAGCGGATGGCCAttt                   >  2:132200/1‑149 (MQ=255)
                                                                                                                         tcccatccTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGCGAACTCGCGGAAGCGGATGGCCATTTc                  <  6:18915/149‑1 (MQ=255)
                                                                                                                             atccTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGCGAACTCGCGGAAGCGGATGGCCATTTCCAcc              >  7:98251/1‑149 (MQ=255)
                                                                                                                              tccTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATc                                                                                               >  6:22848/1‑67 (MQ=255)
                                                                                                                                   gCAGCGGGGTGGACATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGCGAACTCGCGGaa                                 <  7:241158/124‑1 (MQ=255)
                                                                                                                                   gCAGCGGGGTGGACATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGCGAACTCGCGGaa                                 <  7:246155/124‑1 (MQ=255)
                                                                                                                                   gCAGCGGGGTGGACATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGCGAACTCGCGGaa                                 >  8:241158/1‑124 (MQ=255)
                                                                                                                                   gCAGCGGGGTGGACATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGCGAACTCGCGGaa                                 >  8:246155/1‑124 (MQ=255)
                                                                                                                                         gggTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGCGAACTCGCGGAAGCGGATGGCCATTTCCACCTCTTCGGACGaa  >  1:314325/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                              
TGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGCGAACTCGCGGAAGCGGATGGCCATTTCCACCTCTTCGGACGAA  >  NC_002947/3610837‑3611122

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