Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 903,872 (G)5→6 intergenic (‑138/+80) artP ← / ← ybjP arginine ABC transporter ATPase/lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913903,8671.G97.4% 138.4 / ‑2.5 39intergenic (‑133/+85)artP/ybjParginine ABC transporter ATPase/lipoprotein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (21/17);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

TAATTGAATACTCATTGACACTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCT  >  NC_000913/903720‑904009
                                                                                                                                                     |                                                                                                                                              
tAATTGAATACTCATTGACACTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACggg                                                                                                                                             >  1:674709/1‑151 (MQ=255)
 aaTTGAATACTCATTGACACTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACggag                                                                                                                                            <  2:600843/151‑3 (MQ=255)
  aTTGAATACTCATTGACACTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACggggg                                                                                                                                           <  2:590667/151‑1 (MQ=255)
   ttGAATACTCATTGACACTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACgggggg                                                                                                                                          >  1:151514/1‑150 (MQ=255)
    tGAATACTCATTGACACTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACgggggg                                                                                                                                          <  2:6768/150‑2 (MQ=255)
     gAATACTCATTGACACTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACggggggcg                                                                                                                                        >  2:622781/1‑148 (MQ=255)
     gAATACTCATTGACACTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACggggggcg                                                                                                                                        >  1:422593/1‑148 (MQ=255)
                gACACTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGa                                                                                                                              <  1:177529/151‑1 (MQ=255)
                   acTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGc                                                                                                                            >  2:361028/1‑150 (MQ=255)
                            cTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCg                                                                                                                  >  2:351921/1‑151 (MQ=255)
                                  gTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTc                                                                                                            <  2:15066/151‑1 (MQ=255)
                                      gATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTc                                                                                                         <  2:181604/150‑1 (MQ=255)
                                      gATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCa                                                                                                        <  2:138056/151‑1 (MQ=255)
                                           cTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGaa                                                                                                    >  1:204485/1‑150 (MQ=255)
                                             aTTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAAt                                                                                                  <  1:584609/150‑1 (MQ=255)
                                                gAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCg                                                                                              >  2:421434/1‑151 (MQ=255)
                                                     cGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACa                                                                                          <  2:395090/150‑1 (MQ=255)
                                                     cGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACa                                                                                          >  2:9927/1‑150 (MQ=255)
                                                        aGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAAttt                                                                                      >  2:512852/1‑151 (MQ=255)
                                                                   ttGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTAt                                                                           >  1:99731/1‑151 (MQ=255)
                                                                   ttGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTAt                                                                           >  1:269752/1‑151 (MQ=255)
                                                                     gACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTAtt                                                                          <  1:347816/150‑1 (MQ=255)
                                                                        aGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGtt                                                                      <  1:74580/151‑1 (MQ=255)
                                                                        aGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGtt                                                                      >  1:27293/1‑151 (MQ=255)
                                                                                tATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAgg                                                              >  2:28882/1‑151 (MQ=255)
                                                                                 aTGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAgg                                                              <  2:586174/150‑1 (MQ=255)
                                                                                 aTGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGc                                                             <  2:114429/151‑1 (MQ=255)
                                                                                  tGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGAGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTTGTTGATCa                                                                 <  1:130952/146‑1 (MQ=255)
                                                                                    cAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGtt                                                          <  2:178292/151‑1 (MQ=255)
                                                                                     aaTATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTa                                                         <  2:638507/151‑1 (MQ=255)
                                                                                               cTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCg                                               >  1:351919/1‑151 (MQ=255)
                                                                                               cTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCg                                               >  2:600729/1‑151 (MQ=255)
                                                                                                  taaaagttaaaagCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATg                                            >  2:487630/1‑151 (MQ=255)
                                                                                                         taaaagCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGAc                                     >  1:477957/1‑151 (MQ=255)
                                                                                                           aaagcaaagcGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACgg                                   >  2:283048/1‑151 (MQ=255)
                                                                                                           aaagcaaagcGCATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACgg                                   <  2:509861/151‑1 (MQ=255)
                                                                                                               cacagcGCATTATTCAATAAAGATAGCACAAAATAACGGGGGGCGGTGGTCGGCTAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGCTCAAGGCGATAGCGATTCTCGATCGACTGACGGAGc                                >  1:617024/4‑150 (MQ=255)
                                                                                                                    cgcATTATTCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTAccc                          <  1:569791/151‑1 (MQ=255)
                                                                                                                       attattCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACccgcc                       >  1:418934/1‑151 (MQ=255)
                                                                                                                          attCAATAAACATAGCACAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGt                     >  2:327251/1‑150 (MQ=255)
                                                                                                                              aaaaaaCATAGCAAAAAAAAAGGGGGGGCGGTGGTGGGCGAGAATAAATGTCGGCATTCCTCACAAAATGCCGGAAATTTTGCGGGGTTTATTGTTTGATAAGGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCAt                <  1:143838/148‑1 (MQ=255)
                                                                                                                                   aCATAGCAAAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGTTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACg           <  2:208451/151‑1 (MQ=255)
                                                                                                                                    cATAGCAAAAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCTTCCCTCACGAAATGCCGAACATTTTACGGGTTTTTTTGTTTGATCAAGGCTTTAGCGATTCTCGTTGGACTGCCGGAGCGTACCCGCCGTGGCATGAACGc          <  2:572665/151‑1 (MQ=255)
                                                                                                                                         cacaAAATAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCg     >  1:552695/1‑151 (MQ=255)
                                                                                                                                            aaaaTAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCt  >  1:27214/1‑151 (MQ=255)
                                                                                                                                             aaaTAACGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCt  >  2:43411/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
TAATTGAATACTCATTGACACTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAACGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCT  >  NC_000913/903720‑904009

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