Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 707,910 A→G 100% intergenic (‑1/‑121) mntB ← / → sirR ABC superfamily ATP binding cassette transporter, ABC protein/possible metal‑dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730707,9100AG94.1% 105.9 / ‑1.4 34intergenic (‑1/‑121)mntB/sirRABC superfamily ATP binding cassette transporter, ABC protein/possible metal‑dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (11/21);  total (12/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAACT  >  CP000730/707782‑708039
                                                                                                                                |                                                                                                                                 
gAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATAcctc                                                                                                                       >  2:216034/1‑141 (MQ=255)
 aaGATTTCCCAGCACCGTTCGGGCGCATTATAGCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACAGGCTTAGTACCTAAAAACAGATTTAAATCTGGTGTTTATAGCAAGCGTTTACAcctcc                                                                                                                      <  2:66990/141‑1 (MQ=255)
 aaGATTTACCAGCTCCATTCGGGCCCATGATACTAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTATCTTATAACAGATTTAAATCTGTTGTTTCTAACAAGTGTTTATAcctcc                                                                                                                      <  2:268988/141‑1 (MQ=255)
 aaGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACATACGTTTATAcctcc                                                                                                                      <  2:173793/141‑1 (MQ=255)
 aaGAGTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATAcctcc                                                                                                                      <  2:305724/141‑1 (MQ=255)
     tttAGCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAAt                                                                                                                  <  1:68333/141‑1 (MQ=255)
      ttACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAAtt                                                                                                                 <  2:25632/141‑1 (MQ=255)
        aCCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTaa                                                                                                               <  1:306485/141‑1 (MQ=255)
          cAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTaaaa                                                                                                             >  1:212243/1‑141 (MQ=255)
                 aTTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAgg                                                                                                      <  1:216034/141‑1 (MQ=255)
                       gCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAAcc                                                                                                <  2:167816/141‑1 (MQ=255)
                       gCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAAcc                                                                                                <  2:53143/141‑1 (MQ=255)
                       gCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAAcc                                                                                                <  2:221048/141‑1 (MQ=255)
                          cATGATACCAATTATTTCTTGGTTTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCtaa                                                                                             <  1:253785/141‑1 (MQ=255)
                                aCCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTaattaatt                                                                                       >  1:116034/1‑141 (MQ=255)
                                             ccgcGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTg                                                                          <  1:68424/141‑1 (MQ=255)
                                                  taatgatataGATAAGGAAATGTTTTTAAGTACATGCTAAATAAAAAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTCATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAat                                                                     <  2:249493/135‑1 (MQ=255)
                                                                 ggAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAg                                                      <  2:212243/141‑1 (MQ=255)
                                                                   aaaTGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTa                                                    <  2:230281/141‑1 (MQ=255)
                                                                   aaaTGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTAGCATGTCAAGTa                                                    <  2:67567/141‑1 (MQ=255)
                                                                          tttAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAGATTCGt                                             <  1:82329/141‑1 (MQ=255)
                                                                          tttAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGt                                             >  1:67257/1‑141 (MQ=255)
                                                                           ttAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTg                                            >  2:141845/1‑141 (MQ=255)
                                                                           ttAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTg                                            <  1:119879/141‑1 (MQ=255)
                                                                               gTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGatat                                        <  1:173008/141‑1 (MQ=255)
                                                                                        tattaCCTAAAAACAAATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTAGCATGTCAAGTAGATTCGTTTAATAATATAGAc                               <  1:245654/141‑1 (MQ=255)
                                                                                                 aaaaCAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGTCAatgtatgt                      >  2:127500/1‑141 (MQ=255)
                                                                                                      aGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTgag                                                                        >  2:180812/1‑86 (MQ=255)
                                                                                                      aGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTgag                                                                        <  1:180812/86‑1 (MQ=255)
                                                                                                      aGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGt                 >  1:137937/1‑141 (MQ=255)
                                                                                                       gATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGAGATAATATAGACAATGTATGTGAGGTg                >  1:178199/1‑141 (MQ=255)
                                                                                                           tAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAg            >  1:31473/1‑141 (MQ=255)
                                                                                                                    tGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAAc   >  2:301841/1‑141 (MQ=255)
                                                                                                                     gTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAACt  >  2:319786/1‑141 (MQ=255)
                                                                                                                                |                                                                                                                                 
GAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAACT  >  CP000730/707782‑708039

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