Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,176,896 C→A 9.9% Y511* (TAC→TAA PP_3676 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,176,8960CA 9.9% 121.4 / 2.0 40Y511* (TAC→TAAPP_3676hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (21/15);  new base A (2/2);  total (23/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 = 9.94e-01

AAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGAAGGTGAACCTTGGCAAGTTGCGGCATCCGACTTCACTACAAAAGAGCAAATGCAGGACTTCA  >  NC_002947/4176747‑4177038
                                                                                                                                                     |                                                                                                                                               
aaaTTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTAc                                                                                                                                                 <  1:930910/150‑1 (MQ=255)
  aTTTGAGGCAGCTGGCGTCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTACCCGACTTATTGGCAAATACATGTCAAGTTGTAGAGATTTGTTACAt                                                                                                                                               >  2:345971/1‑150 (MQ=255)
         gCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAg                                                                                                                                        >  1:836911/1‑150 (MQ=255)
            gCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTgg                                                                                                                                     <  1:19289/150‑1 (MQ=255)
                         gcgGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCt                                                                                                                        >  1:646426/1‑150 (MQ=255)
                           ggAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTgg                                                                                                                      <  2:635390/150‑1 (MQ=255)
                                 cTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACt                                                                                                                >  1:664341/1‑150 (MQ=255)
                                  tCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTa                                                                                                               >  2:688342/1‑150 (MQ=255)
                                            ccGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAg                                                                                                     <  2:79254/150‑1 (MQ=255)
                                                aTTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGAATTCTTAGGAta                                                                                                 >  2:922646/1‑150 (MQ=255)
                                                                     cTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTCTCAGACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATtctc                                                                            <  1:922646/150‑1 (MQ=255)
                                                                                        gcagcaTTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGAAGGTGAAc                                                        <  1:453244/150‑1 (MQ=255)
                                                                                            caTTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGAAGGTGAACCTTg                                                    <  2:342755/150‑1 (MQ=255)
                                                                                            caTTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGAAGGTGAACCTTg                                                    <  1:349773/150‑1 (MQ=255)
                                                                                             aTTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTg                                                                                                                       >  1:738908/1‑83 (MQ=255)
                                                                                             aTTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTg                                                                                                                       <  2:738908/83‑1 (MQ=255)
                                                                                              ttCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGCTTTGTTACATCCCTAAGTGGTAAAAATACCGCTGGGGCACTAGACTTCTTAGCATACTACAACAGAAATTTATTCTCGTAGCAACTGAAGGTGACCCTTGGc                                                  >  2:571190/1‑150 (MQ=255)
                                                                                              ttCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGCAGGTGaaacctcga                                                  >  1:661714/1‑143 (MQ=255)
                                                                                                 aaaGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGAAGGTGAACCTTGGCAAg                                               >  2:838420/1‑150 (MQ=255)
                                                                                                 aaaGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCACTGGTGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGAAGGTGAACCTTGGCAAg                                               >  1:980001/1‑150 (MQ=255)
                                                                                                   agagTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTAtt                                                                               <  2:101784/117‑1 (MQ=255)
                                                                                                   agagTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTAtt                                                                               >  1:101784/1‑117 (MQ=255)
                                                                                                           tGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGAAgg                                                             <  1:505143/126‑1 (MQ=255)
                                                                                                           tGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGAAgg                                                             >  2:505143/1‑126 (MQ=255)
                                                                                                               ttACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGAAGGTGAACCTTGGCAAGTTGCGGCATCCGAc                                 >  1:559552/1‑150 (MQ=255)
                                                                                                                   tGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGaaa                                                                                     <  1:909145/95‑1 (MQ=255)
                                                                                                                   tGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGaaa                                                                                     >  2:909145/1‑95 (MQ=255)
                                                                                                                   tGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGAAGGTGAACCTTGGCAAGTTGCGGCATCCGACTTCa                             >  1:1122872/1‑150 (MQ=255)
                                                                                                                     ggAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTAt                                                                                <  2:33513/98‑1 (MQ=255)
                                                                                                                     ggAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTAt                                                                                >  1:33513/1‑98 (MQ=255)
                                                                                                                      gAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGAAGGTGAACCTTGGCAAGTTGCGGCATCCGACTTCACTa                          <  2:3119/150‑1 (MQ=255)
                                                                                                                          tACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGAAGGTGAACCTTGGCAAGTTGCGGCATCCGACTTCACTACaaa                      >  2:959633/1‑150 (MQ=255)
                                                                                                                             aTGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGAAGGTGAACCTTGGCAAGTTGCGGCATCCGACTTCACTACaaaa                     >  2:725894/1‑148 (MQ=255)
                                                                                                                             aTGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGAAGGTGAACCTTGGCAAGTTGCGGCATCCGACTTCACTACaaaa                     <  1:725894/148‑1 (MQ=255)
                                                                                                                                ccAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGAAGGTGAACCTTGGCAAGTTGCGGCATCCGACTTCACTACAAAAGAGCa                >  1:254512/1‑150 (MQ=255)
                                                                                                                                 cAAGTTGTAGAGATTTGTTAAATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGAAGGTGAACCTTGGCAAGTTGCGGCATCCGACTTCACTACAAAAGAGCaa               >  1:454763/1‑150 (MQ=255)
                                                                                                                                 cAAGTTGTAGAGATTTGTTAAATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGAAGGTGAACCTTGGCAAGTTGCGGCATCCGACTTCACTACAAAAGAGCaa               >  1:1028067/1‑150 (MQ=255)
                                                                                                                                           agaTTTGTTACATCCCTAAGTGGTTAAAATAACTCATGGTCTCTAGACTTCTTATGATCCTACAACAGAAATTCATTTATTCTGCTGGGCGTCGGTTACCTTGGCAAGTTGCGGCATCCGACTTCACTACAAAAGAGCAAATTCAGGACt     >  1:1102277/1‑150 (MQ=255)
                                                                                                                                             aTTTGTTAAATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGAAGGTGAACCTTGGCAAGTTGCGGCATCCGACTTCACTACAAAAGAGCAAATGCAGGACTTc   <  2:454763/150‑1 (MQ=255)
                                                                                                                                             aTTTGTTAAATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGAAGGTGAACCTTGGCAAGTTGCGGCATCCGACTTCACTACAAAAGAGCAAATGCAGGACTTc   <  2:1028067/150‑1 (MQ=255)
                                                                                                                                              tttGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGAAGGTGAACCTTGGCAAGTTGCGGCATCCGACTTCACTACAAAAGAGCAAATGCAGGACTTCa  <  2:836911/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
AAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGTTACATCCCTAAGTGGTTAAAATACCGCTGGGGCACTAGACTTCTTAGGATACTACAACAGAAATTTATTCTCGTAGCAACTGAAGGTGAACCTTGGCAAGTTGCGGCATCCGACTTCACTACAAAAGAGCAAATGCAGGACTTCA  >  NC_002947/4176747‑4177038

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: