Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 707,879 A→T L11I (TTA→ATA)  mntB ← ABC superfamily ATP binding cassette transporter, ABC protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730707,8790AT94.1% 103.1 / ‑1.8 34L11I (TTA→ATA) mntBABC superfamily ATP binding cassette transporter, ABC protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (15/17);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CCGGTAGCATTAAATTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAG  >  CP000730/707743‑708008
                                                                                                                                        |                                                                                                                                 
ccGGTAGCATTAAATTCACCAATTACAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAAc                                                                                                                               <  2:19913/141‑1 (MQ=255)
ccGGTAGCATTAAATTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCTGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAAc                                                                                                                               <  2:149401/141‑1 (MQ=255)
          tAAATTCATAAAATAAAGACTTGATAATGGAAGATTTACCAGCACTATTCGGTCCCATGATACAAATTATTTCGCCGCTTACTGATATCGATAAAGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAAt                                                                                                                     <  2:79204/141‑1 (MQ=255)
          tAAATTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAAt                                                                                                                     <  2:55740/141‑1 (MQ=255)
           aaaTTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATc                                                                                                                    <  2:232254/141‑1 (MQ=255)
            aaTTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCt                                                                                                                   >  2:146563/1‑141 (MQ=255)
              tccacaAATTACAGACTTCATAAGGGACGATTTACCCGCGCCGTTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCttt                                                                                                                 <  2:127323/139‑1 (MQ=255)
                 aCCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGt                                                                                                              >  2:155494/1‑141 (MQ=255)
                  ccAATTAAAGACTTGATAAGGGAAGATTTACCAGCACAATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGtt                                                                                                             >  1:16192/1‑141 (MQ=255)
                   cAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATAGATAAGGAAATGTTTTTAAGTACATGACTATTAACTATAAACAGATTTAAATCTTTTGttt                                                                                                            >  1:213578/1‑141 (MQ=255)
                          aGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACaa                                                                                                     >  1:300035/1‑141 (MQ=255)
                          aGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTATAACaa                                                                                                     >  2:191519/1‑141 (MQ=255)
                           gACTTTATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACaaa                                                                                                    <  2:247474/141‑1 (MQ=255)
                           gACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACaaa                                                                                                    <  2:259291/141‑1 (MQ=255)
                                      ggAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATAcct                                                                                         <  2:146273/141‑1 (MQ=255)
                                        aaGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATAcctcc                                                                                       <  2:16192/141‑1 (MQ=255)
                                          gATTTACCAGCACCATTGGGGCCCATGATACCAGTTATTTCGCCGCGTACTGGTATCGACGAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTa                                                                                     >  2:39516/1‑141 (MQ=255)
                                          gATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTa                                                                                     >  2:30964/1‑141 (MQ=255)
                                              tACCAGCACCATTAGGACCAAGTGCACCAAGTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCT‑ATAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTaa                                                                                >  2:306045/1‑141 (MQ=255)
                                                      ccATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTa                                                                         >  1:173758/1‑141 (MQ=255)
                                                                     aTACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAATTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCtaattaa                                                          >  2:1818/1‑141 (MQ=255)
                                                                          aaTTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGta                                                     <  2:80476/141‑1 (MQ=255)
                                                                                    ccgcGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTg                                           <  1:155494/141‑1 (MQ=255)
                                                                                       cgTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTgaga                                        <  1:229783/141‑1 (MQ=255)
                                                                                             ggTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATAttt                                  >  2:94860/1‑141 (MQ=255)
                                                                                              aaaacGATAAGTAAAGTTACTTAACTACATGCATATCACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACATTAATACATCGGAATTAAAAATTTAGGCTAACCTAATTAATTTTATAATAAACTGAGAAGAATTa                                 <  2:227176/137‑1 (MQ=255)
                                                                                                        ggAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAg                       >  2:39793/1‑141 (MQ=255)
                                                                                                        ggAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAg                       <  1:81210/141‑1 (MQ=255)
                                                                                                            aTGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTaaa                   >  2:178196/1‑141 (MQ=255)
                                                                                                             tGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAAt                  <  1:146563/141‑1 (MQ=255)
                                                                                                                        aCATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGatataa       >  1:213211/1‑141 (MQ=255)
                                                                                                                             cTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGtata                                                   <  2:113680/92‑1 (MQ=255)
                                                                                                                             cTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGtata                                                   >  1:113680/1‑92 (MQ=255)
                                                                                                                             cTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAg  <  1:1818/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                 
CCGGTAGCATTAAATTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAG  >  CP000730/707743‑708008

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: ‑