Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,301,115 C→T 10.9% A79V (GCC→GTC)  PP_3774 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,301,1150CT10.9% 132.5 / 2.7 51A79V (GCC→GTC) PP_3774hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/20);  new base T (3/3);  total (28/23)
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

TGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGTTTC  >  NC_002947/4300975‑4301261
                                                                                                                                            |                                                                                                                                                  
tGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCa                                                                                                                                           <  1:444800/150‑1 (MQ=255)
           ccAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGATGCGAATGCGTGTAAACGCACAACCCGCTGGTGGCGAGGGTCCAGGATCACTGGAATCGcc                                                                                                                                >  2:599772/1‑150 (MQ=255)
            cAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCa                                                                                                                               <  2:664099/150‑1 (MQ=255)
                             tCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCa                                                                                                              >  2:1244337/1‑150 (MQ=255)
                              cAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCaa                                                                                                             <  2:747011/150‑1 (MQ=255)
                                              cAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAg                                                                                             >  2:479189/1‑150 (MQ=255)
                                                    tcAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAggtg                                                                                       <  1:479189/150‑1 (MQ=255)
                                                      aGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAggtggt                                                                                     >  2:448937/1‑150 (MQ=255)
                                                         cccAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCAt                                                                                  >  2:672053/1‑150 (MQ=255)
                                                          ccaaccaaAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATc                                                                                 >  2:956426/1‑150 (MQ=255)
                                                           caaccaaAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCg                                                                                >  1:1173828/1‑150 (MQ=255)
                                                               caaAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCgg                                                                                                  >  1:1105456/1‑128 (MQ=255)
                                                               caaAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCgg                                                                                                  <  2:1105456/128‑1 (MQ=255)
                                                                aaaGGGGCTTCGAGCAAAAGCACCCACAGGACACGGTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAg                                                                           >  1:103254/1‑150 (MQ=255)
                                                                 aaGGGGCTTCGAGCAAAAGCCACCACAGGAAACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGa                                                                          <  2:1173828/150‑1 (MQ=255)
                                                                   ggggCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACt                                                                                                     <  1:548503/121‑1 (MQ=255)
                                                                   ggggCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACt                                                                                                     >  2:548503/1‑121 (MQ=255)
                                                                              cAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCAAGATGCACGTACGCgg                                                             <  1:956426/150‑1 (MQ=255)
                                                                               aaaaGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCggg                                                            >  2:1010069/1‑150 (MQ=255)
                                                                                 aaGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGga                                                          <  1:1010069/150‑1 (MQ=255)
                                                                                   gCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGTACGTACGCGGGgacg                                                        >  1:1096681/1‑150 (MQ=255)
                                                                                   gCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGgacg                                                        >  2:864096/1‑150 (MQ=255)
                                                                                     aCCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCCGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGgacgac                                                      <  1:672053/150‑1 (MQ=255)
                                                                                      cccACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGATAATGACGCTGACCAAGGTCATGTCGGCACATAAGTTTGTCATCGCCCAGATGCACGTACGCGGGGacgacg                                                     <  1:864096/150‑1 (MQ=255)
                                                                                       ccACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATc                                                                                 <  1:165825/121‑1 (MQ=255)
                                                                                       ccACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATc                                                                                 >  2:165825/1‑121 (MQ=255)
                                                                                       ccACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACgg                                                    >  1:631327/1‑150 (MQ=255)
                                                                                         acaGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGcc                                                  >  1:778305/1‑150 (MQ=255)
                                                                                          caGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGccc                                                 >  1:1153051/1‑150 (MQ=255)
                                                                                           aGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCccg                                                >  1:352131/1‑150 (MQ=255)
                                                                                            ggACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCccgc                                               >  1:307135/1‑150 (MQ=255)
                                                                                             gACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAggtgg                                                                                      <  1:45194/110‑1 (MQ=255)
                                                                                             gACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAggtgg                                                                                      >  2:45194/1‑110 (MQ=255)
                                                                                               cacGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTTATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCgg                                            <  1:1224005/150‑1 (MQ=255)
                                                                                                 cGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCATACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGATCATCGCCCAGATGCACGTACGCGGGGACGACGGCccgc                                               >  2:924172/1‑145 (MQ=255)
                                                                                                 cGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCccgc                                               <  1:924172/145‑1 (MQ=255)
                                                                                                   cTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTaaa                                        <  2:760019/150‑1 (MQ=255)
                                                                                                    tGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGaa                                                                                           >  1:24613/1‑98 (MQ=255)
                                                                                                    tGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGaa                                                                                           <  2:24613/98‑1 (MQ=255)
                                                                                                    tGCGAATGGACGGAAACGCAAGACTGGAAGCCGGCCTCCGTCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAg                                       <  2:103254/150‑1 (MQ=255)
                                                                                                       gAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGgtg                                    <  2:1153051/150‑1 (MQ=255)
                                                                                                               ggAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGaa                            >  2:508696/1‑150 (MQ=255)
                                                                                                               ggAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGaa                            >  1:577117/1‑150 (MQ=255)
                                                                                                                    cGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAgg                       >  1:61481/1‑150 (MQ=255)
                                                                                                                    agtgatccTGGGCGCTGGCCAGCGTCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAgg                       >  1:1080127/8‑150 (MQ=255)
                                                                                                                      cAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGTACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGc                     <  2:1096681/150‑1 (MQ=255)
                                                                                                                         tccTGGGCGCTGGCCAGCGTCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAAc                  <  2:1080127/148‑1 (MQ=255)
                                                                                                                         gACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAAc                  >  1:1195306/1‑150 (MQ=255)
                                                                                                                            tGGGCGCTGGCCAGCGTCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAgg                       >  2:373892/1‑142 (MQ=255)
                                                                                                                            tGGGCGCTGGCCAGCGTCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAgg                       <  1:373892/142‑1 (MQ=255)
                                                                                                                               aCGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATCAcc            <  1:469982/150‑1 (MQ=255)
                                                                                                                                         gcgcCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGTTTc  <  1:508696/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
TGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGTTTC  >  NC_002947/4300975‑4301261

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