Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 973,326 C→G 33.3% G26R (GGC→CGC)  elyC ← envelope biogenesis factor; DUF218 superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913973,3260CG33.3% 74.0 / 55.0 54G26R (GGC→CGC) elyCenvelope biogenesis factor; DUF218 superfamily protein
Reads supporting (aligned to +/- strand):  ref base C (19/17);  new base G (5/13);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATATTACACCAGACGGGTTC  >  NC_000913/973179‑973466
                                                                                                                                                   |                                                                                                                                            
aggcatttacgggagaaaatatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCgccaggcc                                                                                                                                              <  8:302529‑M2/56‑1 (MQ=255)
       tGTCGGATAGGTGCTTTCGATGGGACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCa                                                                                                                                       <  3:98343/148‑1 (MQ=255)
              gaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATca                                                                                                                                <  2:288304‑M2/70‑1 (MQ=255)
               aGGTGCTTTCGATGGGACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATca                                                                                                                                <  4:91916/147‑1 (MQ=255)
               aGGTGCTTTCGATGGGACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATcaa                                                                                                                               >  6:286024/1‑148 (MQ=255)
                ggTGCTTTCGATGGGACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAAAAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATcaa                                                                                                                               <  4:35205/147‑1 (MQ=255)
                    atcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAacaaca                                                                                                                          >  7:190160‑M2/73‑148 (MQ=255)
                            gggACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCgg                                                                                                                   >  5:7953/1‑147 (MQ=255)
                            aaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGa                                                                                                                  <  8:190160‑M2/84‑1 (MQ=255)
                             ggACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGa                                                                                                                  <  1:180842/147‑1 (MQ=255)
                             ggACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGa                                                                                                                  <  4:62061/147‑1 (MQ=255)
                             ggACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGa                                                                                                                  <  5:286024/147‑1 (MQ=255)
                             ggACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGa                                                                                                                  <  6:7953/147‑1 (MQ=255)
                             aacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCgg                                                                                                                   >  7:11862‑M2/64‑146 (MQ=255)
                                  gaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTgg                                                                                                             >  4:145359‑M2/59‑147 (MQ=255)
                                  gaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTgg                                                                                                             >  4:140572‑M2/59‑147 (MQ=255)
                                      aaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGgcagc                                                                                                         >  1:64944‑M2/55‑147 (MQ=255)
                                        atgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAAc                                                                                                      <  2:64944‑M2/96‑1 (MQ=255)
                                           gTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGa                                                                                                                  <  3:253876/133‑1 (MQ=255)
                                           gTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGa                                                                                                                  >  4:253876/1‑133 (MQ=255)
                                           gTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACata                                                                                                   >  5:163633/1‑148 (MQ=255)
                                              ggCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTg                                                                                                <  2:244383/148‑1 (MQ=255)
                                                          ccgagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCacttat                                                                                    >  2:121181‑M2/35‑146 (MQ=255)
                                                            gagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACtttttt                                                                                  >  1:212028‑M2/33‑148 (MQ=255)
                                                            aaCTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACtttttt                                                                                  >  2:9042/1‑148 (MQ=255)
                                                             aCTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACtttttt                                                                                  >  2:40453/1‑147 (MQ=255)
                                                                 aaTAACAAGAGCTCCAGACACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAAGAACTTCTGCGGTAGTGGCAGCAACATATTGCCTATCACTTTTTTCAgt                                                                              <  8:83806/147‑1 (MQ=255)
                                                                   tAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAgtgt                                                                            >  7:306366/1‑147 (MQ=255)
                                                                     aCAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTaa                                                                          <  1:17181/147‑1 (MQ=255)
                                                                       ccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCTCAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTaaaa                                                                        <  3:140572‑M2/126‑1 (MQ=255)
                                                                       ccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTaaaa                                                                        <  3:145359‑M2/126‑1 (MQ=255)
                                                                          cgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTACACCACAGCAGATCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGACAATCACTTTTTTCAGTGTAAAAAGCa                                                                    <  5:39181‑M2/130‑1 (MQ=255)
                                                                           gcgcCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCAt                                                                   <  1:9042/148‑1 (MQ=255)
                                                                                agccaCCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCttt                                                              <  5:104762/148‑1 (MQ=255)
                                                                                 gccaCCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTc                                                             <  6:163633/148‑1 (MQ=255)
                                                                                  ccaCCCAATGCTAATAAAAATCCTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCt                                                            >  7:44345/1‑148 (MQ=255)
                                                                                      ccAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAGCATCCtt                                                               >  8:235919/1‑140 (MQ=255)
                                                                                      ccAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAGCATCCtt                                                               <  7:235919/140‑1 (MQ=255)
                                                                                          tGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCtt                                                    >  2:320679/1‑148 (MQ=255)
                                                                                            cTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACttttt                                                                                   <  5:263825/115‑1 (MQ=255)
                                                                                            cTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACttttt                                                                                   >  6:263825/1‑115 (MQ=255)
                                                                                            cTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCtttt                                                  >  8:141549/1‑148 (MQ=255)
                                                                                              aataaAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCttttt                                                 >  1:18158/1‑147 (MQ=255)
                                                                                               ataaAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGAGTAAAAAGCATCCTTTCTGGTTCCTTTTTTa                                               <  1:40453/148‑1 (MQ=255)
                                                                                               ataaAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCtttttt                                                >  2:12115/1‑147 (MQ=255)
                                                                                                 aaaaaTCTTCCCGGGTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTACTGTTTGAAc                                             >  6:223336/1‑148 (MQ=255)
                                                                                                   aaaTCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTAACCat                                          <  7:141549/148‑1 (MQ=255)
                                                                                                       cTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCata                                         >  2:153364/1‑146 (MQ=255)
                                                                                                             ggTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAgg                                 <  8:306366/148‑1 (MQ=255)
                                                                                                                 ttttGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGAt                              <  7:286277/147‑1 (MQ=255)
                                                                                                                        aaCGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGc                       >  4:50976/1‑147 (MQ=255)
                                                                                                                             cTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACCTCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTCAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATAtt                 <  5:223336/148‑1 (MQ=255)
                                                                                                                             cTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATAtt                 <  8:11862/148‑1 (MQ=255)
                                                                                                                                    acaGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATATTACACCAg          <  1:320679/148‑1 (MQ=255)
                                                                                                                                    acaGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATATTACACCAg          <  1:153364/148‑1 (MQ=255)
                                                                                                                                            gccaggccagCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATATTACACCAGACGGGTTc  <  3:50976/148‑1 (MQ=255)
                                                                                                                                                ggccagCGCCAATAATCAACAACATCAGCGGGAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGata                   <  1:135457/127‑1 (MQ=255)
                                                                                                                                                ggccagCGCCAATAATCAACAACATCAGCGGGAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGata                   >  2:135457/1‑127 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
TATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATATTACACCAGACGGGTTC  >  NC_000913/973179‑973466

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: