Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,497,350 T→C 10.5% E215G (GAA→GGA)  PP_5651 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,497,3500TC10.5% 123.9 / 3.8 46E215G (GAA→GGA) PP_5651hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (18/23);  new base C (3/2);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TAAGAGTTTCTTGACTTAACTTCATTGGATGAACGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGTATAAAGACTATTCCATGTCCCGGGCTTATATCTTGCAAAGGAAGCACACAGGTCAAT  >  NC_002947/4497204‑4497492
                                                                                                                                                  |                                                                                                                                              
tAAGAGTTTCTTGACTTAACTTCATTGGATGAACGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCt                                                                                                                                             >  2:1396411/1‑150 (MQ=255)
 aaGAGTTTCTTGACTTAACTTCATTGGATGAACGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCtt                                                                                                                                            <  1:1396411/150‑1 (MQ=255)
   gagTTTCTTGACTTAACTTCATTGGATGAACGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTAt                                                                                                                                          <  1:1073227/150‑1 (MQ=255)
    agTTTCTTGACTTAACTTCATTGGATGAACGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTAtt                                                                                                                                         <  2:1071245/150‑1 (MQ=255)
      tttCTTGACTTAACTTCATTGGATGAACGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTgg                                                                                                                                       >  1:891554/1‑150 (MQ=255)
             aCTTAACTTCATTGGATGAACGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTCCCTTATTGGTCCGGGt                                                                                                                                >  1:1049163/1‑150 (MQ=255)
             aCTTAACTTCATTGGATGAACGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTCCCTTATTGGTCCGGGt                                                                                                                                >  1:339272/1‑150 (MQ=255)
                       aTTGGATGAACGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGaaa                                                                                                                      >  2:772625/1‑150 (MQ=255)
                               aaCGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCa                                                                                                              >  1:468558/1‑150 (MQ=255)
                                aCGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCTGGTCgat                                                                                                                            >  1:1354326/1‑133 (MQ=255)
                                aCGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCTGGTCgat                                                                                                                            <  2:1354326/135‑3 (MQ=255)
                                 cGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGCCTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAg                                                                                                            <  1:549943/150‑1 (MQ=255)
                                            tCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTCCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGa                                                                                                 <  2:1049163/150‑1 (MQ=255)
                                            tCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTCCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGa                                                                                                 <  2:339272/150‑1 (MQ=255)
                                               acacCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGaaaa                                                                                                                     <  1:1412285/127‑1 (MQ=255)
                                               acacCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGaaaa                                                                                                                     >  2:1412285/1‑127 (MQ=255)
                                                   ccATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTAt                                                                                                                                          <  1:895050/102‑1 (MQ=255)
                                                   ccATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTAt                                                                                                                                          >  2:895050/1‑102 (MQ=255)
                                                   ccATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTCCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAgg                                                                                          >  2:853469/1‑150 (MQ=255)
                                                       aCGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTg                                                                                      <  1:227278/150‑1 (MQ=255)
                                                           ccAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCggg                                                                                                                                 <  2:367418/103‑1 (MQ=255)
                                                           ccAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCggg                                                                                                                                 >  1:367418/1‑103 (MQ=255)
                                                               tCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGca                                                                                 <  1:1208749/147‑1 (MQ=255)
                                                               tCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGca                                                                                 >  2:1208749/1‑147 (MQ=255)
                                                               tCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACaa                                                                              >  1:23896/1‑150 (MQ=255)
                                                                      ttCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCt                                                                       <  2:468558/150‑1 (MQ=255)
                                                                           tCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACa                                                                           <  2:1226052/141‑1 (MQ=255)
                                                                           tCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACa                                                                           >  1:1226052/1‑141 (MQ=255)
                                                                             aacaaATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAAtt                                                                >  2:431784/1‑150 (MQ=255)
                                                                                  aTCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGt                                                           <  2:891554/150‑1 (MQ=255)
                                                                                                      cAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTcc                                                                                                                                              <  2:1322547/47‑1 (MQ=255)
                                                                                                      cAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTcc                                                                                                                                              >  1:1322547/1‑47 (MQ=255)
                                                                                                            aTTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGTATAAAGACTATTCCATGTCCCGGGCt                                 <  1:908539/150‑1 (MQ=255)
                                                                                                              ttCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGTATAAAGGCTATTCCATGTCCCGGGCTta                               <  1:853469/150‑1 (MQ=255)
                                                                                                               tCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGTATAAAGACTATTCCATGTCCCGGGCTtat                              <  1:431784/150‑1 (MQ=255)
                                                                                                                 aTTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGTATAAAGACTATTCCATGTCCCGGGCTTatat                            >  2:1384114/1‑150 (MQ=255)
                                                                                                                 aTTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGTATAAAGACTATTCCATGTCCCGGGCTTatat                            >  2:752313/1‑150 (MQ=255)
                                                                                                                    aaCTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGTATAAAGACTATTCCATGTCCCGGGCTTATATCtt                         <  1:1384114/150‑1 (MQ=255)
                                                                                                                       tGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGt                                                                                    >  2:598964/1‑88 (MQ=255)
                                                                                                                       tGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGt                                                                                    <  1:598964/88‑1 (MQ=255)
                                                                                                                       tGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGtata                                                        >  1:1009303/1‑116 (MQ=255)
                                                                                                                       tGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGtata                                                        <  2:1009303/116‑1 (MQ=255)
                                                                                                                          aaaaaaGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGTATAAAGACTATTCCATGTCCCGGGCTTATATCTTGCAAAg                   <  1:772625/150‑1 (MQ=255)
                                                                                                                           aaaaaGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCa                                                             >  2:289390/1‑107 (MQ=255)
                                                                                                                           aaaaaGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCa                                                             <  1:289390/107‑1 (MQ=255)
                                                                                                                                           gACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGTATAAAGACTATTCCATGTCCCGGGCTTATATCTTGCAAAGGAAGCACACAGGTCAAt  <  2:23896/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
TAAGAGTTTCTTGACTTAACTTCATTGGATGAACGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGTATAAAGACTATTCCATGTCCCGGGCTTATATCTTGCAAAGGAAGCACACAGGTCAAT  >  NC_002947/4497204‑4497492

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