Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,826,830 G→C 12.4% W81C (TGG→TGC)  PP_2479 → cytochrome c family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,826,8300GC12.4% 153.8 / 11.5 55W81C (TGG→TGC) PP_2479cytochrome c family protein
Reads supporting (aligned to +/- strand):  ref base G (22/26);  new base C (5/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACAGCCTGCAACCGG  >  NC_002947/2826684‑2826978
                                                                                                                                                  |                                                                                                                                                    
tACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGAc                                                                                                                                                    <  6:261092/149‑1 (MQ=255)
tACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGAc                                                                                                                                                    <  3:260550/149‑1 (MQ=255)
    tGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGAt                                                                                                                                                >  1:68403/1‑149 (MQ=255)
        gcgcgcTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGTCATCGGTCCATGCACCGATGCCg                                                                                                                                            <  6:249049/149‑1 (MQ=255)
         cgcgcTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGa                                                                                                                                           <  2:99962/149‑1 (MQ=255)
          gcgcTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAg                                                                                                                                          >  6:247408/1‑149 (MQ=255)
                ccAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTc                                                                                                                                    <  2:68403/149‑1 (MQ=255)
                       tgtgGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTgcgc                                                                                                                             >  2:125/1‑149 (MQ=255)
                                 cATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGAc                                                                                                                                             <  8:165257/123‑1 (MQ=255)
                                 cATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGAc                                                                                                                                             >  7:165257/1‑123 (MQ=255)
                                 cATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCgg                                                                                                                   >  5:61109/1‑149 (MQ=255)
                                   tACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGAt                                                                                                                                                <  4:261770/118‑1 (MQ=255)
                                       gTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTgg                                                                                                             <  6:61109/149‑1 (MQ=255)
                                           cTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGcc                                                                                                           >  7:20637/1‑147 (MQ=255)
                                           cTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGcc                                                                                                           <  8:20637/147‑1 (MQ=255)
                                             ggCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGa                                                                                                       >  4:209048/1‑149 (MQ=255)
                                                ggTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACgg                                                                                                    >  6:30877/1‑149 (MQ=255)
                                                      cccTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCa                                                                                              >  3:155096/1‑149 (MQ=255)
                                                       ccTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCac                                                                                             <  4:155096/149‑1 (MQ=255)
                                                           cGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCacta                                                                                           <  1:125/147‑1 (MQ=255)
                                                                gtggtgTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCg                                                                                    <  5:247408/149‑1 (MQ=255)
                                                                 tggtgTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCgg                                                                                   >  6:170637/1‑149 (MQ=255)
                                                                    tgtTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGt                                                                                >  3:10346/1‑149 (MQ=255)
                                                                      ttGAGTTCAAGCTGCCGTTCGGCTAGCTGTTTTCGCCCAAGATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCTATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGttt                                                                              <  3:216872/149‑1 (MQ=255)
                                                                        gAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTcc                                                                            >  2:38452/1‑149 (MQ=255)
                                                                          gTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTccc                                                                           <  2:86839/148‑1 (MQ=255)
                                                                           ttCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTcc                                                                            <  2:186710/146‑1 (MQ=255)
                                                                           ttCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTcc                                                                            >  1:186710/1‑146 (MQ=255)
                                                                           ttCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCta                                                                         <  1:38452/149‑1 (MQ=255)
                                                                                     cGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTAc                                                               <  5:30877/149‑1 (MQ=255)
                                                                                           gCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCac                                                                                             >  4:15550/1‑113 (MQ=255)
                                                                                                   ttttCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTcgcg                                                 >  1:57985/1‑149 (MQ=255)
                                                                                                       cGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGa                                              <  1:134183/148‑1 (MQ=255)
                                                                                                         cccAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTcc                                                                            <  1:268062/116‑1 (MQ=255)
                                                                                                         cccAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTcc                                                                            >  2:268062/1‑116 (MQ=255)
                                                                                                          ccAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGa                                           >  4:110312/1‑148 (MQ=255)
                                                                                                             aCATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATc                                       >  5:71962/1‑149 (MQ=255)
                                                                                                                ttACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTc                                    >  5:85915/1‑149 (MQ=255)
                                                                                                                ttACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGCCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTc                                    <  5:170637/149‑1 (MQ=255)
                                                                                                                    cccGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTcgc                                                  <  4:71571/131‑1 (MQ=255)
                                                                                                                    cccGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTcgc                                                  >  3:71571/1‑131 (MQ=255)
                                                                                                                    cccGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGa                                >  8:35659/1‑149 (MQ=255)
                                                                                                                        gACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCaa                            <  3:209048/149‑1 (MQ=255)
                                                                                                                          cccGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAgg                          <  6:85915/149‑1 (MQ=255)
                                                                                                                           ccGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTGCTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAggg                         <  2:56525/149‑1 (MQ=255)
                                                                                                                                aaCCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGGCATCCTCGCGATCAAGGGTTAtc                    <  8:59088/149‑1 (MQ=255)
                                                                                                                                aaCCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTAtc                    <  6:71962/149‑1 (MQ=255)
                                                                                                                                 ctcgtcATCGGTCCATGCACCGATGCCGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTAtct                   >  2:95823/6‑149 (MQ=255)
                                                                                                                                 ctcgtcATCGGTCCATGCACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTAtct                   >  5:196818/6‑149 (MQ=255)
                                                                                                                                  tcgtcATCGGTCCATGCACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGAt                                             <  2:19607/118‑1 (MQ=255)
                                                                                                                                  tcgtcATCGGTCCATGCACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGAt                                             >  1:19607/5‑122 (MQ=255)
                                                                                                                                    gtcATCGGTCCATGCACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGa                >  5:51916/3‑149 (MQ=255)
                                                                                                                                    gtcATCGGTCCATGCACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGa                >  5:12718/3‑149 (MQ=255)
                                                                                                                                        acGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATAGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTCTACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACAGc            >  6:129555/2‑149 (MQ=255)
                                                                                                                                               aTGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACAGCCTGCAAc     <  4:10346/149‑1 (MQ=255)
                                                                                                                                                  gACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACAGCCTGCAACCgg  >  3:152840/1‑149 (MQ=255)
                                                                                                                                                  gACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACAGCCTGCAACCg   <  3:110312/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACAGCCTGCAACCGG  >  NC_002947/2826684‑2826978

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