Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 185,090 C→T 8.0% D175N (GAC→AAC)  PP_0161 ← putative transmembrane sensor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947185,0900CT 8.0% 142.9 / 2.1 50D175N (GAC→AAC) PP_0161putative transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base C (23/23);  new base T (2/2);  total (25/25)
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.59e-01

CCAGTTCTGGCCACCTTCGATCCACTGAAGTTGGCCAGCACCAGGCCGGTGTATCGCCACTTTGCCGTGGCTGACGCTGACACGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAGGTCCACGGCACTGCAGGTATTGAGTTGCATCAGCGAGCCATCGGG  >  NC_002947/184947‑185231
                                                                                                                                               |                                                                                                                                             
ccAGTTCTGGCCACCTTCGATCCACTGAAGTTGGCCAGCACCAGGCCGGTGTATCGCCACTTTGCCGTGGCTGACGCTGACACGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGgggg                                                                                                                                          >  6:214133/1‑149 (MQ=255)
            aCCTTCGATCCACTGAAGTTGGCCAGCACCAGGCCGGTGTATCGCCACTTTGCCGTGGCTGACGCTGACACGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAg                                                                                                                                >  7:128167/1‑147 (MQ=255)
              cTTCGATCCACTGAAGTTGGCCAGCACCAGGCCGGTGTATCGCCACTTTGCCGTGGCTGACGCTGACACGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAgg                                                                                                                            <  5:193695/149‑1 (MQ=255)
                       aCTGAAGTTGGCCAGCACCAGGCCGGTGTATCGCCACTTTGCCGTGGCTGACGCTGACACGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCgg                                                                                                                   >  7:146836/1‑149 (MQ=255)
                        cTGAAGTTGGCCAGCACCAGGCCGGTGTATCGCCACTTTGCCGTGGCTGACGCTGACACGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCggg                                                                                                                  <  8:128167/149‑1 (MQ=255)
                            aGTTGGCCAGCACCAGGCCGGTGTATCGCCACTTTCCCGTGGCTGACGCTGACACGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGtt                                                                                                              <  8:241247/149‑1 (MQ=255)
                                ggCCAGCACCAGGCCGGTGTATCGCCACTTTGCCGTGGCTGACGCTGACACGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCa                                                                                                           <  2:117807/148‑1 (MQ=255)
                                  ccAGCACCAGGCCGGTGTATCGCCACTTTGCCGTGGCTGACGCTGACACGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGtt                                                                                                        >  6:263417/1‑149 (MQ=255)
                                      cACCAGGCCGGTGTATCGCCACTTTGCCGTGGCTGACGCTGACACGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGatc                                                                                                    <  2:208639/149‑1 (MQ=255)
                                      cACCAGGCCGGTGTATCGCCACTTTGCCGTGGCTGACGCTGACACGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGat                                                                                                     >  7:194165/1‑148 (MQ=255)
                                                     tCGCCACTTTGCCGTGGCTGACGCTGACACGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGc                                                                                     <  7:292265/149‑1 (MQ=255)
                                                       gCCACTTTGCCGTGGCTGACGCTGACACGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCtt                                                                                   >  6:212548/1‑149 (MQ=255)
                                                        ccACTTTGCCGTGGCTGACGCTGACACGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTTGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCa                                                                                                           <  5:58613/124‑1 (MQ=255)
                                                        ccACTTTGCCGTGGCTGACGCTGACACGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTTGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCa                                                                                                           >  6:58613/1‑124 (MQ=255)
                                                               gCCGTGCCTTCCGATTACACGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCATTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATTCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCgg                                                                            <  8:279685/148‑1 (MQ=255)
                                                                ccGTGGCTGACGCTGACACGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGAc                                                                          >  2:97884/1‑149 (MQ=255)
                                                                    ggCTGACGCTGACACGTGTGCAGTCGCTATACTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGTTGCCGGGGGTTGCAGTTGATCTTCAATTCCCCTAGGTTCAGGCGGACCAGa                                                                      >  1:306126/1‑149 (MQ=255)
                                                                              gacACGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGa                                                             >  7:18925/1‑148 (MQ=255)
                                                                               acacGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGa                                                             <  5:214133/147‑1 (MQ=255)
                                                                               acacGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATc                                                           <  8:18925/149‑1 (MQ=255)
                                                                                  cGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATg                                                        <  5:263417/149‑1 (MQ=255)
                                                                                         aGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCggggg                                                                                                                >  4:100361/1‑86 (MQ=255)
                                                                                         aGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCggggg                                                                                                                <  3:100361/86‑1 (MQ=255)
                                                                                             gCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTc                                                                                  >  4:293785/1‑112 (MQ=255)
                                                                                             gCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTc                                                                                  <  3:293785/112‑1 (MQ=255)
                                                                                             gCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAGGt                                             <  8:194165/149‑1 (MQ=255)
                                                                                              cTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGtt                                                               >  6:92490/1‑130 (MQ=255)
                                                                                              cTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGtt                                                               <  5:92490/130‑1 (MQ=255)
                                                                                                   cTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTTGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAgg                                              <  1:235523/142‑1 (MQ=255)
                                                                                                   cTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTTGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAgg                                              >  2:235523/1‑142 (MQ=255)
                                                                                                         ggCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAGGTCCACGGCACTGc                                 >  7:242193/1‑149 (MQ=255)
                                                                                                          gCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAGGTCCACGGCACTGCa                                <  2:173328/149‑1 (MQ=255)
                                                                                                           ccACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAGGTCCACGGCACTGCAg                               >  7:272410/1‑149 (MQ=255)
                                                                                                            cacaAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAg                                                                                      <  8:232713/93‑1 (MQ=255)
                                                                                                            cacaAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAg                                                                                      >  7:232713/1‑93 (MQ=255)
                                                                                                                aaGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAGGTCCACGGCACTGCAGGTAtt                          <  5:212548/149‑1 (MQ=255)
                                                                                                                   cGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAGGTCCACGGCACTGCAGGTATTGAg                       <  8:146836/149‑1 (MQ=255)
                                                                                                                     tCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAGGTCCACGGCACTGCAGGTATTGAGtt                     >  3:79548/1‑149 (MQ=255)
                                                                                                                       cTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAGGTCCACGGCACTGCAGGTATTGAGTTg                    <  8:272410/148‑1 (MQ=255)
                                                                                                                        ttCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAGGTCCACGGCACTGCAGGTATTGAGTTGCa                  >  1:251447/1‑149 (MQ=255)
                                                                                                                          cGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAGGTCCACTGCACTGCAGGTATTGAGTTGCAt                 >  5:65310/1‑148 (MQ=255)
                                                                                                                          cGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAGGTCCACGGCACTGCAGGTATTGAGTTGCAt                 >  4:260380/1‑148 (MQ=255)
                                                                                                                           gAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAGGTCCACGGCACTGCAGGTATTGAGTTGCATc                <  3:260380/148‑1 (MQ=255)
                                                                                                                             aGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAGGTCCACGGCACTGCAGGTATTGAGTTGCATCAGc             <  2:251447/149‑1 (MQ=255)
                                                                                                                               ccTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAGGTCCACGGCACTGCAGGTATTGAGTTGCATCAGc             >  8:161321/1‑147 (MQ=255)
                                                                                                                               ccTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAGGTCCACGGCACTGCAGGTATTGAGTTGCATCAGCGa           <  1:97884/149‑1 (MQ=255)
                                                                                                                                cTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAGGTCCACGGCACTGc                                 <  3:1518/126‑1 (MQ=255)
                                                                                                                                cTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAGGTCCACGGCACTGc                                 >  4:1518/1‑126 (MQ=255)
                                                                                                                                cTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAGGTCCACGGCACTGCAGGTATTGAGTTGCATCAGCGAg          >  6:176022/1‑149 (MQ=255)
                                                                                                                                          cgcgTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAGGTCCACGGCACTGCAGGTATTGAGTTGCATCAGCGAGCCATCggg  >  6:86497/1‑147 (MQ=255)
                                                                                                                                               cgcgGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAg                                               <  2:212275/97‑1 (MQ=255)
                                                                                                                                               cgcgGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAg                                               >  1:212275/1‑97 (MQ=255)
                                                                                                                                               |                                                                                                                                             
CCAGTTCTGGCCACCTTCGATCCACTGAAGTTGGCCAGCACCAGGCCGGTGTATCGCCACTTTGCCGTGGCTGACGCTGACACGTGTGCAGTCGCTATCCTGATAGGCCACAAAGCGTCCTTCGAAGCCTTCCAGCAGCGCGTCGCGGGTTTCTACCAGTAGGGGATGCCGGGGGTTGCAGTTGATCATCAATTCCCCAAGCTTCAGGCGGACCAGACGCTGTTGATCATGGAAGGCCAGGTCCACGGCACTGCAGGTATTGAGTTGCATCAGCGAGCCATCGGG  >  NC_002947/184947‑185231

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