Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,448,548 A→C T7T (ACT→ACG)  rpsQ ← 30S ribosomal subunit protein S17

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,448,5480AC98.0% 177.3 / ‑4.1 51T7T (ACT→ACG) rpsQ30S ribosomal subunit protein S17
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (22/28);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

GATACCGCATTCGTTGTTCTCGTCATGTACGTGCAGTTTGGTCGTACGCTTGATGAATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGAGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTGAGACTGTTGCAGCTGGCCACTTGCAGCCTGCATACGCAGGTTGAACTGCTCACGCA  >  NC_000913/3448404‑3448695
                                                                                                                                                |                                                                                                                                                   
gATACCGCATTCGTTGTTCTCGTCATGTACGTGCAGTTTGGTCGTACGCTTGATGAATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACgg                                                                                                                                               <  2:477149/151‑1 (MQ=255)
 aTACCGCATTCGTTGTTCTCGTCATGTACGTGCAGTTTGGTCGTACGCTTGATGAATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGa                                                                                                                                              <  2:461060/151‑1 (MQ=255)
   aCCGCATTCGTTGTTCTCGTCATGTACGTGCAGTTTGGTCGTACGCTTGATGAATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGCTAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGAtt                                                                                                                                            >  1:480799/1‑151 (MQ=255)
   aCCGCATTCGTTGTTCTCGTCATGTACGTGCAGTTTGGTCGTACGCTTGATGAATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGAtt                                                                                                                                            <  1:819048/151‑1 (MQ=255)
     cGCATTCGTTGTTCTCGTCATGTACGTGCAGTTTGGTCGTACGCTTGATGAATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGGACGGAtttt                                                                                                                                          >  1:595897/1‑151 (MQ=255)
          tCGTTGTTCTCGTCATGTACGTGCAGTTTGGTCGTACGCTTGATGAATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCgg                                                                                                                                     >  2:786734/1‑151 (MQ=255)
            gttgttCTCGTCATGTACGTGCAGTTTGGTCGTACGCTTGATGAATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTc                                                                                                                                   <  2:586012/151‑1 (MQ=255)
                 tctcGTCATGTACGTGCAGTTTGGTCGTACGCTTGATGAATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTa                                                                                                                               <  1:814770/150‑1 (MQ=255)
                 tctcGTCATGTACGTGCAGTTTGGTCGTACGCTTGATGAATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTAc                                                                                                                              >  2:793512/1‑151 (MQ=255)
                 tctcGTCATGTACGTGCAGTTTGGTCGTACGCTTGATGAATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTAc                                                                                                                              >  2:625945/1‑151 (MQ=255)
                  ctcGTCATGTACGTGCAGTTTGGTCGTACGCTTGATGAATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACg                                                                                                                             >  2:320124/1‑151 (MQ=255)
                    cGTCATGTACGTGCAGTTTGGTCGTACGCTTGATGAATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACg                                                                                                                             <  1:910498/149‑1 (MQ=255)
                       cATGTACGTGCAGTTTGGTCGTACGCTTGATGAATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCAccc                                                                                                                        >  2:212566/1‑151 (MQ=255)
                            aCGTGCAGTTTGGTCGTACGCTTGATGAATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCtt                                                                                                                   >  2:661056/1‑151 (MQ=255)
                              gTGCAGTTTGGTCGTACGCTTGATGAATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTtct                                                                                                                 >  2:872162/1‑151 (MQ=255)
                               tGCAGTTTGGTCGTACGCTTGATGAATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTtctc                                                                                                                >  2:329032/1‑151 (MQ=255)
                                    tttGGTCGTACGCTTGATGAATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTc                                                                                                            >  2:130236/1‑150 (MQ=255)
                                     ttGGTCGTACGCTTGATGAATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCa                                                                                                           <  2:257786/150‑1 (MQ=255)
                                      tGGTCGTACGCTTGATGAATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAg                                                                                                          <  2:699927/150‑1 (MQ=255)
                                                  tgatgaATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAAc                                                                                              <  1:53868/150‑1 (MQ=255)
                                                      gaATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTg                                                                                         >  2:894128/1‑151 (MQ=255)
                                                         tttACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCg                                                                                       >  2:837166/1‑150 (MQ=255)
                                                              cGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATcg                                                                                 <  2:256982/151‑1 (MQ=255)
                                                                 aGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATcgcg                                                                               <  2:379246/150‑1 (MQ=255)
                                                                        gtgtTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACtt                                                                       <  2:199916/151‑1 (MQ=255)
                                                                        gtgtTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACtt                                                                       <  1:133814/151‑1 (MQ=255)
                                                                             tCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTc                                                                  >  1:387731/1‑151 (MQ=255)
                                                                              cacaAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCa                                                                 <  2:46437/151‑1 (MQ=255)
                                                                               acaAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCa                                                                 >  1:840897/1‑150 (MQ=255)
                                                                                    aCGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGt                                                           <  1:640904/151‑1 (MQ=255)
                                                                                     cGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTg                                                          <  1:796184/151‑1 (MQ=255)
                                                                                      gTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTga                                                         <  1:32968/151‑1 (MQ=255)
                                                                                        tCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTg                                                          >  1:667538/1‑148 (MQ=255)
                                                                                         cGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTGAGAc                                                      <  2:155296/151‑1 (MQ=255)
                                                                                         cGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTGAGAc                                                      >  2:88765/1‑151 (MQ=255)
                                                                                          gATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTGAGACt                                                     <  2:769210/151‑1 (MQ=255)
                                                                                                      aaTGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTGAGACTGTTGCAGCTGGc                                         <  1:50087/151‑1 (MQ=255)
                                                                                                       aTGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTGAGACTGTTGCAGCTGGcc                                        >  2:387327/1‑151 (MQ=255)
                                                                                                         ggATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTGAGACTGTTGCAGCTGGCCa                                       <  2:692434/150‑1 (MQ=255)
                                                                                                         ggATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTGAGACTGTTGCAGCTGGCCa                                       <  2:799728/150‑1 (MQ=255)
                                                                                                         ggATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGAATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTGAGACTGTTGCAGCTGGCCAc                                      >  2:680021/1‑151 (MQ=255)
                                                                                                           aTTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTGAGACTGTTGCAGCTGGCCACtt                                    <  1:244199/151‑1 (MQ=255)
                                                                                                             ttCTCCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTGAGACTGTTGCAGCTGGCCACTTGc                                  >  2:117550/1‑151 (MQ=255)
                                                                                                                tcCATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTGAGACTGTTGCAGCTGGCCACTTGCAGc                               >  1:300549/1‑151 (MQ=255)
                                                                                                                 ccATTTTGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTGAGACTGTTGCAGCTGGCCACTTGCAGcc                              <  1:625925/151‑1 (MQ=255)
                                                                                                                  cATTTTGTCGCTAACAACGCGACCTTGCAGAGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTGAGACTGTTGCAGCTGGCCACTTGCAGCCt                             <  2:262659/151‑1 (MQ=255)
                                                                                                                     tttGTCGCTAACAACGCGACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTGAGACTGTTGCAGCTGGCCACTTGCAGCCTGc                           <  2:906312/150‑1 (MQ=255)
                                                                                                                                    gcgACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTGAGACTGTTGCAGCTGGCCACTTGCAGCCTGCATACGCAGGTTGAACt           <  2:416725/151‑1 (MQ=255)
                                                                                                                                     cgACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTGAGACTGTTGCAGCTGGCCACTTGCAGCCTGCATACGCAGGTTGAACt           <  1:501757/150‑1 (MQ=255)
                                                                                                                                      gACCTTGCAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTGAGACTGTTGCAGCTGGCCACTTGCAGCCTGCATACGCAGGTTGAACTGc         >  1:171061/1‑151 (MQ=255)
                                                                                                                                             cAGCGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTGAGACTGTTGCAGCTGGCCACTTGCAGCCTGCATACGCAGGTTGAACTGCTCACgca  <  1:661036/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GATACCGCATTCGTTGTTCTCGTCATGTACGTGCAGTTTGGTCGTACGCTTGATGAATTTACCGTAGATCGGGTGTTTCACAAAACGTTCGATAGCAACAACAATGGATTTCTCCATTTTGTCGCTAACAACGCGACCTTGCAGAGTACGGATTTTATCGGTCATTACGCACCCGCCTTCTCGTTCAGTAAAGTCTTAACGCGTGCGACATCGCGACGCACTTGCTTCAACAGGTGAGACTGTTGCAGCTGGCCACTTGCAGCCTGCATACGCAGGTTGAACTGCTCACGCA  >  NC_000913/3448404‑3448695

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑