Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,174,461 T→C A155A (GCA→GCG U712_10875 ← CCA‑adding enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,174,4610TC98.2% 173.8 / NA 55A155A (GCA→GCGU712_10875CCA‑adding enzyme
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (23/31);  minor base A (1/0);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACTCTATTGTTTTTCGTTCAACTGAAACATGGGATAGGAGTGATTTTTCTTTCGTGATGGCTTCTTTTGTTTCTAGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACAGCTCTCAACATACGAAGTGCATCCTCTTGAAATCTGTCCTCAGGCTTTCCGACGGTTCGAATCAGTTTCTGATCAATATCTTTCTTGCCGCCGAAATAATCAAGCACCTTCCCATCCGCTGTCATGGCCATCGCATTAATCGTCAAATCT  >  CP006881/2174329‑2174593
                                                                                                                                    |                                                                                                                                    
aaaCTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCa                                                                                                                                   <  1:432906/136‑1 (MQ=255)
aaaCTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATc                                                                                                                                 <  8:350820/138‑1 (MQ=255)
aaaCTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATc                                                                                                                                 <  5:275786/138‑1 (MQ=255)
aaaCTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATc                                                                                                                                 <  7:206901/138‑1 (MQ=255)
aaaCTCAATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCa                                                                                                                                   <  6:38443/136‑1 (MQ=255)
      tataGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCtt                                                                                                                            <  2:227551/137‑1 (MQ=255)
        tagTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCtt                                                                                                                            <  6:198120/133‑1 (MQ=255)
        tagTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCtt                                                                                                                            >  5:198120/3‑135 (MQ=255)
          gTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATc                                                                                                                      <  4:293236/139‑1 (MQ=255)
            ttttCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCta                                                                                                                    >  4:307610/1‑138 (MQ=255)
             tttCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCa                                                                                                                                   <  6:325259/123‑1 (MQ=255)
             tttCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCa                                                                                                                                   >  5:325259/1‑123 (MQ=255)
             tttCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATc                                                                                                                                 <  3:57882/125‑1 (MQ=255)
             tttCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATc                                                                                                                                 >  4:57882/1‑125 (MQ=255)
              ttcgttcgACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATc                                                                                                                      <  4:28193/135‑1 (MQ=255)
              ttcgttcgACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATc                                                                                                                      >  3:28193/1‑135 (MQ=255)
              ttcgttcgACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCACTCAGCATACGAAGCGCATCCTATTGAAATCtaac                                                                                                                  >  5:41208/1‑136 (MQ=255)
                cGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCTCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCt                                                                                                                >  1:67286/1‑139 (MQ=255)
                 gTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCt                                                                                                                <  3:307610/138‑1 (MQ=255)
                           aaCATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTcca                                                                                                     <  8:151092/139‑2 (MQ=255)
                            aCATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAACGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTccaa                                                                                                    <  1:48576/139‑3 (MQ=255)
                                gggAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGt                                                                                                >  5:29681/1‑139 (MQ=255)
                                 ggAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGtt                                                                                               >  7:108134/1‑139 (MQ=255)
                                     aGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGaa                                                                                           >  4:113982/1‑139 (MQ=255)
                                     aGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGaa                                                                                           >  7:143146/1‑139 (MQ=255)
                                        aGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCa                                                                                        >  2:300282/1‑139 (MQ=255)
                                           gATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACtt                                                                                     >  7:308100/1‑139 (MQ=255)
                                            aTTTTTCTTTAGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAAt                                                                                                                       <  3:43883/104‑1 (MQ=255)
                                            aTTTTTATTTAGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAAt                                                                                                                       >  4:43883/1‑104 (MQ=255)
                                                  ctttcGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGAt                                                                               <  4:36015/138‑1 (MQ=255)
                                                       gcaaTGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAatat                                                                         <  8:217873/136‑1 (MQ=255)
                                                              cttcttCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATAtctttctt                                                                  <  1:300282/139‑1 (MQ=255)
                                                                    ctGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTgccg                                                              <  3:113982/136‑1 (MQ=255)
                                                                        ttccggTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCaaaa                                                        >  3:388875/6‑139 (MQ=255)
                                                                         tccggTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCaaaa                                                        <  6:401310/134‑1 (MQ=255)
                                                                                            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)
                                                                                                             cGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAg                                                <  4:444800/110‑1 (MQ=255)
                                                                                                             cGTACGGCTCTCAGCATACGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAg                                                >  3:444800/1‑110 (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                                  >  5:70671/1‑107 (MQ=255)
                                                                                                                              aCGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCg                                  <  6:70671/107‑1 (MQ=255)
                                                                                                                              aCGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCAtt               <  7:53361/126‑1 (MQ=255)
                                                                                                                              aCGAAGCGCATCCTCTTGAAATCTATCCTCAGGCTTTCCAACGGTTCGAATCACTTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGTCCGCTGTCATGGCCATCGCAtt               >  8:53361/1‑126 (MQ=255)
                                                                                                                              aCGAAGAGCATCCTCTTGAAATCTATCCTCAGGCTTTACAACGGTTCGATTCAATTTCTGATCAATATCTTTCTTGCCGCCAAAATAATCAAGCACCTTCCCGACCGCTGTCATTGCCATCGCGTTGATCGTTAAAtct  >  1:74929/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                    
AAACTCTATTGTTTTTCGTTCAACTGAAACATGGGATAGGAGTGATTTTTCTTTCGTGATGGCTTCTTTTGTTTCTAGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACAGCTCTCAACATACGAAGTGCATCCTCTTGAAATCTGTCCTCAGGCTTTCCGACGGTTCGAATCAGTTTCTGATCAATATCTTTCTTGCCGCCGAAATAATCAAGCACCTTCCCATCCGCTGTCATGGCCATCGCATTAATCGTCAAATCT  >  CP006881/2174329‑2174593

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