Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,174,557 A→G D123D (GAT→GAC)  U712_10875 ← CCA‑adding enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,174,5570AG98.3% 182.2 / NA 58D123D (GAT→GAC) U712_10875CCA‑adding enzyme
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (29/28);  minor base C (1/0);  total (30/28)
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

CAAGCTGGCTCATAAAGCGTACAGCTCTCAACATACGAAGTGCATCCTCTTGAAATCTGTCCTCAGGCTTTCCGACGGTTCGAATCAGTTTCTGATCAATATCTTTCTTGCCGCCGAAATAATCAAGCACCTTCCCATCCGCTGTCATGGCCATCGCATTAATCGTCAAATCTCTGCGTTTTAAATCCTCTTCTAATGATGTGATAAATTGCACGTCTGACGGTCTTCTATAATCTTCATACTCTGATTCAGCCCGGAATGTAGTGACTT  >  CP006881/2174421‑2174690
                                                                                                                                        |                                                                                                                                     
cAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTc                                                                                                                                     <  6:41208/139‑1 (MQ=255)
  aGCTGGCTCATAAATCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCg                                                                                                                                   <  8:383888/139‑1 (MQ=255)
  aGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCg                                                                                                                                   <  8:308100/139‑1 (MQ=255)
        cTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTc                                                                                                                              <  1:253228/138‑1 (MQ=255)
        cTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTc                                                                                                                              >  2:253228/1‑138 (MQ=255)
           aTAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATAACTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCAt                                                                                                                            <  1:120101/137‑1 (MQ=255)
                 cGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCg                                                                                                                    >  2:277787/1‑139 (MQ=255)
                   tACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTACTGATCATAATCTTTCATGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCa                                                                                                                  >  4:317872/1‑139 (MQ=255)
                      ggCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCAttg                                                                                                               <  4:388875/138‑2 (MQ=255)
                       gCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCAttga                                                                                                              <  2:67286/139‑3 (MQ=255)
                        ctctCAGCCTACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGAt                                                                                                             <  6:29681/139‑1 (MQ=255)
                           tcAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGt                                                                                                          >  7:421926/1‑139 (MQ=255)
                                  aCGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCCCCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAAtct                                                                                                   >  3:425963/1‑139 (MQ=255)
                                  aCGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCg                                                                                                                                   <  6:70671/107‑1 (MQ=255)
                                  aCGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCg                                                                                                                                   >  5:70671/1‑107 (MQ=255)
                                  aCGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCAtt                                                                                                                <  7:53361/126‑1 (MQ=255)
                                  aCGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCAtt                                                                                                                >  8:53361/1‑126 (MQ=255)
                                  aCGAAGAGCATCCTCTTGAAATCTATCCTCAGGCTTTACAACGGTTCGATTCAATTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGACCGCTGTCATTGCCATCGCGTTGATCGTTAAAtct                                                                                                   >  1:74929/1‑139 (MQ=255)
                                       gcgCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCAt                                                                                                                      >  7:380549/3‑115 (MQ=255)
                                       gcgCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCAt                                                                                                                      <  8:380549/113‑1 (MQ=255)
                                         gCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGtt                                                                                            >  5:238198/1‑139 (MQ=255)
                                         gCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCAGCGtt                                                                                            >  3:375124/1‑139 (MQ=255)
                                          cATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGttt                                                                                           >  4:262654/1‑139 (MQ=255)
                                              ctctTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGTTTTaga                                                                                       <  2:60865/139‑3 (MQ=255)
                                              ctctTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGTTTTaga                                                                                       <  8:143146/139‑3 (MQ=255)
                                                         tatCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCAttg                                                                                                               <  3:58828/102‑2 (MQ=255)
                                                         tatCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCAttg                                                                                                               >  4:58828/3‑103 (MQ=255)
                                                         tatCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGAt                                                                                                             >  5:267376/3‑106 (MQ=255)
                                                         tatCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGAt                                                                                                             <  6:267376/104‑1 (MQ=255)
                                                         tatCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGTTTTAGATCCTCTTCTaa                                                                            >  5:338009/3‑139 (MQ=255)
                                                         tatCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGTTTTAGATCCTCTTCTaa                                                                            >  5:404332/3‑139 (MQ=255)
                                                         tatCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTATCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGATGTCATGGCCATCGCATTGATCGTTAAATATCTGCGTGTTAGATCCTCTTCTa                                                                             >  1:199537/3‑138 (MQ=255)
                                                           tCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGTTTTAGATCCTCTTCTAat                                                                           >  4:75016/1‑138 (MQ=255)
                                                                     ttCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGTTTTAGATCCTCTTCTAATGATGAGATaaa                                                                >  3:161212/1‑139 (MQ=255)
                                                                       ccaaCGGTTCGAATCACTTTCTGATCAATATCATTCTTGCAGCCAAAATAATCAAGCACCTTCCCCTCCGCTCTCAGTGCCATCCAAGTGATCGTTAAATCTATGCGTTTTAGATCCTCTTCTAATGATTAGATATAtt                                                              >  6:370323/4‑139 (MQ=255)
                                                                          aCGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGTTTTAGATCCTCTTCTAATGATGAGATAAATTGCa                                                           <  1:277787/139‑1 (MQ=255)
                                                                          aCGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGTTTTAGATCCTCTTCTAATGATGAGATAAATTGCa                                                           <  8:108134/139‑1 (MQ=255)
                                                                              ttCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGATTTATATACTCTTCTAATTATGAGATAAATTGAACTTc                                                       >  1:173784/1‑139 (MQ=255)
                                                                                 gAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGTTTTAGATCCTCTTCTAATGATGAGATAAATTGCACTTCTGa                                                    <  4:161212/139‑1 (MQ=255)
                                                                                        tttCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGtta                                                                                                        >  1:195924/1‑78 (MQ=255)
                                                                                        tttCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGtta                                                                                                        <  2:195924/80‑3 (MQ=255)
                                                                                        tttCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTTCGTTTTAGATCCTCTTCTAATGATGAGATAAATTGCACTTCTGACGGtctt                                             >  3:257117/1‑139 (MQ=255)
                                                                                             gATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGtt                                                                                                         >  1:178692/1‑73 (MQ=255)
                                                                                             gATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGtt                                                                                                         <  2:178692/74‑2 (MQ=255)
                                                                                              aTCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGTTTTAGATCCTCTTCTAATGATGAGATAAATTGCACTTCTGACGGtcttctgaaa                                       <  3:75016/139‑5 (MQ=255)
                                                                                                  atatCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGTTTTAGATc                                                                                     >  6:421419/1‑89 (MQ=255)
                                                                                                  atatCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGTTTTAGATc                                                                                     <  5:421419/89‑1 (MQ=255)
                                                                                                     tctttcttGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGTTTTAGATCCTCTTCTAATGATGAGATAAATTGCACTTCTGACGGTCTTCTGAAATcaaca                                 <  6:238198/138‑5 (MQ=255)
                                                                                                     tctttcttGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCAGCGTTTTAGATCCTCTTCTAATGATGAGATAAATTGCACTTCTGACGGTCTTCTGAAATcaaca                                 <  4:375124/138‑5 (MQ=255)
                                                                                                         tcttGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGGCATGGCCATCGCATTGATCGTTAAATCTCTGCGTTTTAGATCCTCTTCTAATGATGAGATAAATTGCACTTCTGACGGTCTTCTGAAATCAACATaat                             <  8:421926/138‑3 (MQ=255)
                                                                                                             gccgccAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGTTTTAGATCCTCTTCTAATGATGAGATAAATTGCACTTCTGACGGTCTTCTGAAATCAACATAATCAGAt                        >  7:335567/1‑139 (MQ=255)
                                                                                                             gccgccAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGTTTTAGATCCTCTTCTAATGATGAGATAAATTGCACTTCTGACGGTCTTCTGAAATCAACATAATCAGAt                        >  5:220833/1‑139 (MQ=255)
                                                                                                              ccgccaACATAATCAAGGACCTTCACGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGAGTTTTAGATCCTCTTCTAATGATGAGATAAATTGCACTTCTGACGGTCTTCTGAAATCAACATAATCAGAtt                       <  5:370323/139‑1 (MQ=255)
                                                                                                                              gCACCTTCCAGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGTTTTAGATCCTCTTCTAATGATGAGATAAATTGCACTTCTGACGGTCTTCTGAAATCAACATAATCAGATTCAGTCCGGAATGTCGt       <  6:404332/139‑1 (MQ=255)
                                                                                                                                    tCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGTTTTAGATCCTCTTCTAATGATGAGATAAATTGCACTTCTGACGGTCTTCTGAAATcaaca                                 >  5:424959/1‑103 (MQ=255)
                                                                                                                                    tCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGTTTTAGATCCTCTTCTAATGATGAGATAAATTGCACTTCTGACGGTCTTCTGAAATcaaca                                 <  6:424959/107‑5 (MQ=255)
                                                                                                                                    tCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGTTTTAGATCCTCTTCTAATGATGAGATAAATTGCACTTCTGACGGTCTTCTGAAATCAACATAATCAGATTCAGTCCGGAATGTCGTGACtt  >  7:231264/1‑138 (MQ=255)
                                                                                                                                    tCCCGTCCGCTGTCATGGCCATCGCATTGATCGTTAAATCTCTGCGTTTTAGATCCTCTTCTAATGATGAGATAAATTGCACTTCTGACGGTCTTCTGAAATCAACATAATCAGATTCAGTCCGGAATGTCGTGACtt  >  6:257434/1‑138 (MQ=255)
                                                                                                                                        |                                                                                                                                     
CAAGCTGGCTCATAAAGCGTACAGCTCTCAACATACGAAGTGCATCCTCTTGAAATCTGTCCTCAGGCTTTCCGACGGTTCGAATCAGTTTCTGATCAATATCTTTCTTGCCGCCGAAATAATCAAGCACCTTCCCATCCGCTGTCATGGCCATCGCATTAATCGTCAAATCTCTGCGTTTTAAATCCTCTTCTAATGATGTGATAAATTGCACGTCTGACGGTCTTCTATAATCTTCATACTCTGATTCAGCCCGGAATGTAGTGACTT  >  CP006881/2174421‑2174690

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