Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 7,925 +A coding (35/1431 nt) yaaJ ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009137,9221.A97.1% 108.4 / ‑3.5 35coding (38/1431 nt)yaaJputative transporter
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (16/18);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCT  >  NC_000913/7793‑8055
                                                                                                                                   |                                                                                                                                     
ccTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                               >  1:656030/1‑139 (MQ=255)
                        ttAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaagga                                                                                                                                >  1:17595/1‑109 (MQ=255)
                        ttAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaagga                                                                                                                                <  2:17595/113‑5 (MQ=255)
                                            gCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCGGATAAATCATTACCGATCCTCAAAAGGACGCTGTTAATGATGGAGAAAAAATCTGGCATGAATATACCTCTTAt                                                                                   >  2:434629/1‑139 (MQ=255)
                                            gCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTAt                                                                                   >  1:735251/1‑139 (MQ=255)
                                             cGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAAt                                                                                                         <  1:42701/116‑1 (MQ=255)
                                             cGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAAt                                                                                                         >  2:42701/1‑116 (MQ=255)
                                                      cGAAACTGCACAAATCCGGTGCGAAAAGTGAACTAACAACCTGAGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                               <  2:573106/85‑1 (MQ=255)
                                                      cGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                               >  1:573106/1‑85 (MQ=255)
                                                       gAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCTCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTcgcg                                                                        <  1:434629/139‑1 (MQ=255)
                                                             gCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACt                                                                  >  2:438180/1‑139 (MQ=255)
                                                                caAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                               <  1:689698/139‑1 (MQ=255)
                                                                   aTCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                               >  1:709227/1‑72 (MQ=255)
                                                                   aTCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                               <  2:709227/72‑1 (MQ=255)
                                                                             aaaaGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgca                                                                                                   <  2:175805/90‑1 (MQ=255)
                                                                             aaaaGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgca                                                                                                   <  2:649764/90‑1 (MQ=255)
                                                                             aaaaGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgca                                                                                                   >  1:175805/1‑90 (MQ=255)
                                                                             aaaaGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgca                                                                                                   >  1:649764/1‑90 (MQ=255)
                                                                                 gTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAATGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCa                                              <  2:207685/139‑1 (MQ=255)
                                                                                          caacCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaa                                     <  2:463490/139‑1 (MQ=255)
                                                                                            acCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                               <  1:344902/47‑1 (MQ=39)
                                                                                            acCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                               >  2:344902/1‑47 (MQ=39)
                                                                                                 cgcCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGtat                              <  1:438180/139‑1 (MQ=255)
                                                                                                   ccGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGtatat                            <  2:541093/139‑1 (MQ=255)
                                                                                                   ccGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGtatat                            >  1:541093/1‑139 (MQ=255)
                                                                                                     gAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATAc                          >  2:473464/1‑139 (MQ=255)
                                                                                                       agagCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                               <  1:569408/100‑1 (MQ=255)
                                                                                                       agagCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                               >  2:569408/1‑100 (MQ=255)
                                                                                                        gagCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCt                       <  2:656030/139‑1 (MQ=255)
                                                                                                           cAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCg                                                                              >  2:294826/1‑81 (MQ=255)
                                                                                                           cAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCg                                                                              <  1:294826/81‑1 (MQ=255)
                                                                                                                aaaTCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATc                                                                                           <  1:124145/63‑1 (MQ=255)
                                                                                                                aaaTCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATc                                                                                           >  2:124145/1‑63 (MQ=255)
                                                                                                                     aTTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGa          <  1:473464/139‑1 (MQ=255)
                                                                                                                      ttACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGAt         <  1:252660/139‑1 (MQ=255)
                                                                                                                           atccccaAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCg                                                                              <  2:382418/61‑1 (MQ=255)
                                                                                                                            tccccaAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCt  >  2:5318/3‑139 (MQ=255)
                                                                                                                            aTCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCg                                                                              >  1:382418/1‑64 (MQ=255)
                                                                                                                            aTCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTc   >  2:767361/1‑139 (MQ=255)
                                                                                                                             ccccaAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGa                                                                       <  2:198291/68‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                     
CCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCT  >  NC_000913/7793‑8055

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: