Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,041,127 T→A 7.8% Q17H (CAA→CAT)  PP_5293 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,041,1270TA 7.8% 146.0 / 3.2 51Q17H (CAA→CAT) PP_5293hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/22);  new base A (2/2);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

TGAAAGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATA‑GCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGC  >  NC_002947/6040984‑6041272
                                                                                                                                                |                                                                                                                                                 
tGAAAGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATA‑GCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCg                                                                                                                                              >  8:183938/1‑149 (MQ=255)
 gAAAGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATA‑GCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCg                                                                                                                                              >  4:226103/1‑148 (MQ=255)
     gTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGGGACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGAtt                                                                                                                                          >  1:245989/1‑149 (MQ=255)
               ccattGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATA‑GCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTggg                                                                                                                               <  8:43812/146‑1 (MQ=255)
                 tttGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATA‑GCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTggg                                                                                                                               <  3:225492/147‑1 (MQ=255)
                 tttGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATA‑GCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCAAACAGCGCGGAGAGGTGGCAGAACATCGCCCAATGCCGAAATCCGGCAATggg                                                                                                                               <  8:222817/147‑1 (MQ=255)
                    gCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATA‑GCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCggg                                                                                                                           <  4:281118/148‑1 (MQ=255)
                             gTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATA‑GCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACa                                                                                                                 <  2:187110/149‑1 (MQ=255)
                                gATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATA‑GCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAgg                                                                                                               <  2:111615/148‑1 (MQ=255)
                                  tAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATA‑GCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTc                                                                                                            >  1:20521/1‑149 (MQ=255)
                                        ggTCTTGCTCGCGCTTGAGGTGCCATA‑GCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCg                                                                                                      <  4:191464/149‑1 (MQ=255)
                                          tCTTGCTCGCGCTTGAGGTGCCATA‑GCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTCGGCGGGGTGATCGACAGGTTCGAATCGCt                                                                                                    >  4:281714/1‑149 (MQ=255)
                                                tCGCGCTTGAGGTGCCATA‑GCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGcc                                                                                              >  8:53054/1‑149 (MQ=255)
                                                tCGCGCTTGAGGTGCCATA‑GCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGcc                                                                                              >  7:151161/1‑149 (MQ=255)
                                                tCGCGCTTGAGGTGCCATA‑GCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGcc                                                                                              >  3:216054/1‑149 (MQ=255)
                                                  gcgcTTGAGGTGCCATA‑GCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGcccc                                                                                            <  7:53054/149‑1 (MQ=255)
                                                      ttGAGGTGCCATA‑GCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGCCAGGTTCGAATCGCTCATGCCCCGCAc                                                                                        >  5:29712/1‑149 (MQ=255)
                                                             gCCATA‑GCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGc                                                                                                                                     >  8:263757/1‑97 (MQ=255)
                                                             gCCATA‑GCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGc                                                                                                                                     <  7:263757/97‑1 (MQ=255)
                                                              ccATA‑GCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTgg                                                                                >  2:181939/1‑149 (MQ=255)
                                                                        cAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTtca                                                                       <  1:181939/149‑1 (MQ=255)
                                                                          gCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCg                                                                                                      >  6:266953/1‑116 (MQ=255)
                                                                          gCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCg                                                                                                      <  5:266953/116‑1 (MQ=255)
                                                                          gCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTtca                                                                       <  8:153978/147‑1 (MQ=255)
                                                                          gCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTtca                                                                       >  7:153978/1‑147 (MQ=255)
                                                                                ccAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGtt                                                                                                             >  4:178482/1‑103 (MQ=255)
                                                                                ccAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGtt                                                                                                             <  3:178482/103‑1 (MQ=255)
                                                                                   gcagGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGt                                                                           >  8:209861/1‑134 (MQ=255)
                                                                                   gcagGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGt                                                                           <  7:209861/134‑1 (MQ=255)
                                                                                    cagGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGtgc                                                           <  3:281714/149‑1 (MQ=255)
                                                                                    cagGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGtgc                                                           <  8:151161/149‑1 (MQ=255)
                                                                                    cagGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGtgc                                                           >  6:258289/1‑149 (MQ=255)
                                                                                     agGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGtgct                                                          <  3:226103/149‑1 (MQ=255)
                                                                                          cccAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTAGGCGGGGTGATCGACAGGTTCGAATCGCTc                                                                                                   <  8:113267/103‑1 (MQ=255)
                                                                                            cAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGttt                                                                         >  4:223076/1‑127 (MQ=255)
                                                                                            cAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCctgct                                                   >  4:182929/1‑149 (MQ=255)
                                                                                             aaGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCctgctg                                                  >  4:35826/1‑149 (MQ=255)
                                                                                              aGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCctgctg                                                  >  5:216133/1‑148 (MQ=255)
                                                                                                    aCTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTc                                            >  7:90742/1‑148 (MQ=255)
                                                                                                            gCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAAtctc                                   <  4:216054/149‑1 (MQ=255)
                                                                                                                   aGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATg                            <  7:183938/149‑1 (MQ=255)
                                                                                                                     cgcgGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGcc                          <  2:245989/149‑1 (MQ=255)
                                                                                                                              gTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTAATGCCCTTCTgcgc                 >  5:265837/1‑149 (MQ=255)
                                                                                                                                ggCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACGGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCctgctg                                                  <  3:198150/114‑1 (MQ=255)
                                                                                                                                ggCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACGGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCctgctg                                                  >  4:198150/1‑114 (MQ=255)
                                                                                                                                    gAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGc           <  6:216133/149‑1 (MQ=255)
                                                                                                                                         ccGCCCAATGCCGAAATCCGGCAATGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAg        <  6:72524/146‑1 (MQ=255)
                                                                                                                                         ccGCCCAATGCCGAAATCCGGCAATGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAg        >  5:72524/2‑147 (MQ=255)
                                                                                                                                           gCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATg    >  1:270338/1‑149 (MQ=255)
                                                                                                                                           gCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATg    >  4:107531/1‑149 (MQ=255)
                                                                                                                                             ccAATGCCGAAATCCGGCAATGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGc  >  6:172406/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TGAAAGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATA‑GCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGC  >  NC_002947/6040984‑6041272

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