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

TAATTGAATACTCATTGACACTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑C‑GGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCT  >  NC_000913/903720‑904009
                                                                                                                                                     |                                                                                                                                              
tAATTGAATACTCATTGACACTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑C‑ggg                                                                                                                                             >  1:674709/1‑151 (MQ=255)
 aaTTGAATACTCATTGACACTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑C‑ggag                                                                                                                                            <  2:600843/151‑3 (MQ=255)
  aTTGAATACTCATTGACACTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑C‑ggggg                                                                                                                                           <  2:590667/151‑1 (MQ=255)
   ttGAATACTCATTGACACTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑C‑gggggg                                                                                                                                          >  1:151514/1‑150 (MQ=255)
    tGAATACTCATTGACACTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑C‑gggggg                                                                                                                                          <  2:6768/150‑2 (MQ=255)
     gAATACTCATTGACACTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑C‑ggggggcg                                                                                                                                        >  2:622781/1‑148 (MQ=255)
     gAATACTCATTGACACTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑C‑ggggggcg                                                                                                                                        >  1:422593/1‑148 (MQ=255)
                gACACTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGa                                                                                                                              <  1:177529/151‑1 (MQ=255)
                   acTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGc                                                                                                                            >  2:361028/1‑150 (MQ=255)
                            cTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCg                                                                                                                  >  2:351921/1‑151 (MQ=255)
                                  gTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTc                                                                                                            <  2:15066/151‑1 (MQ=255)
                                      gATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTc                                                                                                         <  2:181604/150‑1 (MQ=255)
                                      gATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCa                                                                                                        <  2:138056/151‑1 (MQ=255)
                                           cTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGaa                                                                                                    >  1:204485/1‑150 (MQ=255)
                                             aTTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAAt                                                                                                  <  1:584609/150‑1 (MQ=255)
                                                gAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCg                                                                                              >  2:421434/1‑151 (MQ=255)
                                                     cGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACa                                                                                          <  2:395090/150‑1 (MQ=255)
                                                     cGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACa                                                                                          >  2:9927/1‑150 (MQ=255)
                                                        aGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAAttt                                                                                      >  2:512852/1‑151 (MQ=255)
                                                                   ttGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTAt                                                                           >  1:99731/1‑151 (MQ=255)
                                                                   ttGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTAt                                                                           >  1:269752/1‑151 (MQ=255)
                                                                     gACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTAtt                                                                          <  1:347816/150‑1 (MQ=255)
                                                                        aGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGtt                                                                      <  1:74580/151‑1 (MQ=255)
                                                                        aGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGtt                                                                      >  1:27293/1‑151 (MQ=255)
                                                                                tATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAgg                                                              >  2:28882/1‑151 (MQ=255)
                                                                                 aTGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAgg                                                              <  2:586174/150‑1 (MQ=255)
                                                                                 aTGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGc                                                             <  2:114429/151‑1 (MQ=255)
                                                                                  tGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGAGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTTGTTGATCa                                                                 <  1:130952/146‑1 (MQ=255)
                                                                                    cAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGtt                                                          <  2:178292/151‑1 (MQ=255)
                                                                                     aaTATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTa                                                         <  2:638507/151‑1 (MQ=255)
                                                                                               cTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCg                                               >  1:351919/1‑151 (MQ=255)
                                                                                               cTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCg                                               >  2:600729/1‑151 (MQ=255)
                                                                                                  taaaagttaaaagCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATg                                            >  2:487630/1‑151 (MQ=255)
                                                                                                         taaaagCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGAc                                     >  1:477957/1‑151 (MQ=255)
                                                                                                           aaagcaaagcGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACgg                                   >  2:283048/1‑151 (MQ=255)
                                                                                                           aaagcaaagcGCATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACgg                                   <  2:509861/151‑1 (MQ=255)
                                                                                                               cacagcGCATTATTCAATAAAGATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCTAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGCTCAAGGCGATAGCGATTCTCGATCGACTGACGGAGc                                >  1:617024/4‑150 (MQ=255)
                                                                                                                    cgcATTATTCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTAccc                          <  1:569791/151‑1 (MQ=255)
                                                                                                                       attattCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACccgcc                       >  1:418934/1‑151 (MQ=255)
                                                                                                                          attCAATAAACATAGCACAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGt                     >  2:327251/1‑150 (MQ=255)
                                                                                                                              aaaaaaCATAGCAAAAAAAAAGG‑GGGGGCGGTGGTGGGCGAGAATAAATGTCGGCATTCCTCACAAAATGCCGGAAATTTTGCGGGGTTTATTGTTTGATAAGGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCAt                <  1:143838/148‑1 (MQ=255)
                                                                                                                                   aCATAGCAAAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGTTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACg           <  2:208451/151‑1 (MQ=255)
                                                                                                                                    cATAGCAAAAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCTTCCCTCACGAAATGCCGAACATTTTACGGGTTTTTTTGTTTGATCAAGGCTTTAGCGATTCTCGTTGGACTGCCGGAGCGTACCCGCCGTGGCATGAACGc          <  2:572665/151‑1 (MQ=255)
                                                                                                                                         cacaAAATAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCg     >  1:552695/1‑151 (MQ=255)
                                                                                                                                            aaaaTAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCt  >  1:27214/1‑151 (MQ=255)
                                                                                                                                             aaaTAA‑CGGGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCt  >  2:43411/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
TAATTGAATACTCATTGACACTCGTATACTGGCAGTCTGATAGCTATTGAGGTCGAAGATAGTACCTTTGACAGAATAATTATGCAATATTTCTGCTTTAAAAGTTAAAAGCAAAGCGCATTATTCAATAAACATAGCACAAAATAA‑C‑GGGGGCGGTGGTCGGCGAGCATAAATGTCGGCATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCT  >  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: ‑