Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,127,745 A→T 17.2% intergenic (‑201/‑80) ctaA ← / → ctaB cytochrome oxidase assembly protein/protoheme IX farnesyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,127,7450AT17.2% 66.7 / 6.4 29intergenic (‑201/‑80)ctaA/ctaBcytochrome oxidase assembly protein/protoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base A (13/11);  new base T (3/2);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGCTTAAATATTAATCTCTTTTAGATGCATGATTCTACTTATTAAAATTTTCTAACATATACATTATAGCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTAAGGAGGGGGATTATGAGCAAAGAGCATACTTTGTCACAAAATATTAGCAGAGTTAACTTCAAA  >  CP000730/1127611‑1127875
                                                                                                                                      |                                                                                                                                  
gAAGCTTAAATATTAATCTCTTTTAGATGCATGATTCTACTTATTAAAATTTTCTAACATATACATTATAGCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACaatttaat                                                                                                                              <  1:70877/141‑1 (MQ=255)
 aaGCTTAAATATTAATCTCTTTTAGATGCATGATTCTACTTATTAAAATTTTCTAACATATACATTATAGCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACaatttaatt                                                                                                                             >  2:168124/1‑141 (MQ=255)
 aaGCTTAAATATTAATCTCTTTTAGATGCATGATTCTACTTATTAAAATTTTCTAACATATACATTATAGCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACaatttaatt                                                                                                                             >  1:115340/1‑141 (MQ=255)
 aaGCTTAAATATTAATCTCTTTTAGATGCATGATTCTACTTATTAAAATTTTCTAACATATACATTATAGCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACaatttaatt                                                                                                                             <  2:226265/141‑1 (MQ=255)
  aGCTTAAATATTAATCTCTTTTAGATGCATGATTCTACTTATTAAAATTTTCTAACATATACATTATAGCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACaatttaattt                                                                                                                            >  1:439814/1‑141 (MQ=255)
   gCTTAAATATTAATCTCTTTTAGATGCATGATTCTACTTATTAAAATTTTCTAACATATACATTATAGCAATAAATTAAATGTTCTTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAAtttt                                                                                                                           >  1:169279/1‑141 (MQ=255)
              aaTCTCTTTTAGATGCATGATTCTACTTATTAAAATTTTCTAACATATACATTATAGCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAg                                                                                                                >  1:54368/1‑141 (MQ=255)
                           tGCATGATTCTACTTATTAAAATTTTCTAACATATACATTATAGCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACATTTTAATTTTATTAAACCTAGTCTTTTTATGCaa                                                                                                   >  1:268464/1‑141 (MQ=255)
                             cATGATTCTACTTATTAAAATTTTCTAACATATACATTATAGCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATc                                                                                                 >  1:252956/1‑141 (MQ=255)
                                     tACTTATTAAAATTTTCTAACATATACATTATAGCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGatat                                                                                         <  2:115340/141‑1 (MQ=255)
                                     tACTTATTAAAATTTTCTAACATATACATTATAGCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGatat                                                                                         <  1:136994/141‑1 (MQ=255)
                                             aaaaTTTTCTAACATATACATTATAGCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCt                                                                                 <  2:45599/141‑1 (MQ=255)
                                                          atatACATTATAGCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAAt                                                                    >  2:357719/1‑141 (MQ=255)
                                                              aCATTATAGCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTaa                                                                >  2:400170/1‑141 (MQ=255)
                                                               cATTATAGCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTAag                                                               >  1:110380/1‑141 (MQ=255)
                                                                 ttATAGCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTAagga                                                             >  2:29512/1‑141 (MQ=255)
                                                                  tataGCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTAaggag                                                            <  2:172130/141‑1 (MQ=255)
                                                                    taGCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTAAGGAggg                                                          <  1:283785/141‑1 (MQ=255)
                                                                     aGCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTAAGGAgggg                                                         >  1:210077/1‑141 (MQ=255)
                                                                     aGCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTAAGGAgggg                                                         >  2:406314/1‑141 (MQ=255)
                                                                      gCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTAAGGAggggg                                                        >  1:271580/1‑141 (MQ=255)
                                                                                aaatGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTAAGGAGGGGGATTATGAGCa                                              >  1:382363/1‑141 (MQ=255)
                                                                                     ttATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTAAGGAGGGGGATTATGAGCAAagag                                         <  2:165576/141‑1 (MQ=255)
                                                                                      tATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTAAGGAGGGGGATTATGAGCAAAGAGc                                        <  1:236196/141‑1 (MQ=255)
                                                                                            tgtCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTAAGGAGGGGGATTATGAGCAAAGAGCATACtt                                  >  1:228005/1‑141 (MQ=255)
                                                                                                 caAAATTGACATTTTCTCTTGTCGTTAATTTGTCACATTTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTAAGGAGGGGGATTATGAGCAAAGAGCATACTTAGTca                             >  1:334822/1‑141 (MQ=255)
                                                                                                   aaaTTGACATTTTCTCTTGTCGTTAATTTGTCACATTTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTAAGGAGGGGGATTATGAGCAAAGAGCATACTTTGTcaca                           <  1:48785/141‑1 (MQ=255)
                                                                                                     aTTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTAAGGAGGGGGATTATGAGCAAAGAGCATACTTTGTCACaaa                         <  1:135184/141‑1 (MQ=255)
                                                                                                       tGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTAAGGAGGGGGATTATGAGCAAAGAGCATACTTTGTCACAAAat                       <  1:410300/141‑1 (MQ=255)
                                                                                                                     gTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTAAGGAGGGGGATTATGAGCAAAGAGCATACTTTGTCACAAAATATTAGCAGAGTTaa         <  2:439814/141‑1 (MQ=255)
                                                                                                                       cGTTAATTTGTCACATTTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTACATGAACAAATTTAAGGAGGGGGATTATGAGCAAAGAGCATACTTTGTCACAAACTATTAGCAGAGTTAACt       >  1:222010/1‑141 (MQ=255)
                                                                                                                         ttAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTAAGGAGGGGGATTATGAGCAAAGAGCATACTTTGTCACAAAATATTAGCAGAGTTAACTTc     <  2:225100/141‑1 (MQ=255)
                                                                                                                           aaTTTGTCACATTTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTAAGGAGGGGGATTATGAGCAAAGAGCATACTTTGTCACAAAATATTAGCAGAGTTAACGTCac   <  2:334822/141‑2 (MQ=255)
                                                                                                                            aTTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTAAGGAGGGGGATTATGAGCAAAGAGCATACTTTGTCACAAAATATTAGCAGAGTTAACTTCaaa  >  1:431226/1‑141 (MQ=255)
                                                                                                                                      |                                                                                                                                  
GAAGCTTAAATATTAATCTCTTTTAGATGCATGATTCTACTTATTAAAATTTTCTAACATATACATTATAGCAATAAATTAAATGTTATTTGTGTCACAAAATTGACATTTTCTCTTGTCGTTAATTTGTCACAATTTAATTTTATTAAACCTAGTCTTTTTATGCAATCTAAGATATCATAATCTTATATGAACAAATTTAAGGAGGGGGATTATGAGCAAAGAGCATACTTTGTCACAAAATATTAGCAGAGTTAACTTCAAA  >  CP000730/1127611‑1127875

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