Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,266,039 G→A 5.9% intergenic (+123/+166) rdlC → / ← chaA sRNA antisense regulator affects LdrC translation; proposed addiction module in LDR‑C repeat, with toxic peptide LdrC/calcium/sodium:proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,266,0390GA 5.9% 2266.5 / 76.4 732intergenic (+123/+166)rdlC/chaAsRNA antisense regulator affects LdrC translation; proposed addiction module in LDR‑C repeat, with toxic peptide LdrC/calcium/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base G (235/428);  new base A (29/13);  total (265/467)
Fisher's exact test for biased strand distribution p-value = 2.33e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 748 total aligned reads displayed.

TAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCA  >  CP009273/1265972‑1266106
                                                                   |                                                                   
tAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATg                                                                     >  1:15946433/1‑68 (MQ=40)
 aaCGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGa                                                                    <  1:32108815/68‑1 (MQ=40)
 aaCGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGa                                                                    <  1:17022656/68‑1 (MQ=40)
 aaCGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGa                                                                    <  1:28697063/68‑1 (MQ=40)
 aaCGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGa                                                                    <  1:27736133/68‑1 (MQ=40)
    gCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGGAATGAaga                                                                 >  1:33263596/1‑68 (MQ=255)
    gCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAaga                                                                 >  1:35750592/1‑68 (MQ=40)
     ccGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAagag                                                                >  1:18073250/1‑68 (MQ=40)
      cGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGt                                                               <  1:7721361/68‑1 (MQ=40)
              tgtgGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCg                                                       >  1:17182988/1‑68 (MQ=40)
               gtgGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCg                                                       >  1:24156089/1‑67 (MQ=40)
                   aTACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCttt                                                  >  1:3464016/1‑68 (MQ=40)
                   aTACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCttt                                                  >  1:3391730/1‑68 (MQ=40)
                   aTACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCttt                                                  >  1:11711623/1‑68 (MQ=40)
                    tactacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCttt                                                  >  1:34212437/1‑67 (MQ=40)
                    tactacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCttt                                                  >  1:35279166/1‑67 (MQ=40)
                    tactacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCttt                                                  >  1:7269448/1‑67 (MQ=40)
                    tactacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCttt                                                  >  1:16725014/1‑67 (MQ=40)
                    tactacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCttt                                                  >  1:32382250/1‑67 (MQ=40)
                        acGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAAc                                             >  1:31109577/1‑68 (MQ=40)
                        acGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAAc                                             >  1:31496419/1‑68 (MQ=40)
                        acGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAAc                                             >  1:8562816/1‑68 (MQ=40)
                        acGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAAc                                             >  1:1096508/1‑68 (MQ=40)
                        acGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCAGGCTTTTCAAc                                             >  1:9743090/1‑68 (MQ=255)
                         cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACa                                            <  1:14656942/68‑1 (MQ=40)
                          gTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACa                                            >  1:6347092/1‑67 (MQ=40)
                            ctcAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCt                                         <  1:2938000/66‑1 (MQ=14)
                            ctcAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCt                                         <  1:9527125/66‑1 (MQ=14)
                            cGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCt                                         <  1:36360314/68‑1 (MQ=39)
                            cGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCt                                         <  1:35445663/68‑1 (MQ=39)
                            cGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCt                                         <  1:19662425/68‑1 (MQ=39)
                             gCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACa                                            >  1:21279121/1‑64 (MQ=37)
                               aaTTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGtt                                      <  1:5045247/68‑1 (MQ=35)
                               aaTTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGtt                                      >  1:13365001/1‑68 (MQ=35)
                                aTTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGtt                                      >  1:26586747/1‑67 (MQ=35)
                                aTTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGtt                                      >  1:26119878/1‑67 (MQ=35)
                                 ttCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGc                                    <  1:29182857/68‑1 (MQ=35)
                                  tCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCa                                   <  1:20941694/68‑1 (MQ=35)
                                  tCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCa                                   >  1:15033703/1‑68 (MQ=35)
                                  tCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCa                                   >  1:13922002/1‑68 (MQ=35)
                                    aTCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACCGCTGTTGCAGt                                 <  1:5777251/68‑1 (MQ=255)
                                      cttcacttcaTCCCTGAAATGTTTGCAATTAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTgg                               <  1:17782728/68‑1 (MQ=18)
                                       ttcacttcaTCCCTGAAATGTTTGCAATTAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGt                              <  1:17485894/68‑1 (MQ=25)
                                       ttcacttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGt                              <  1:10325281/68‑1 (MQ=40)
                                          acttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTc                           >  1:30553756/1‑68 (MQ=40)
                                           cttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCg                          >  1:25601214/1‑68 (MQ=255)
                                           cttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCg                          <  1:23117948/68‑1 (MQ=255)
                                           cttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCg                          >  1:16404921/1‑68 (MQ=14)
                                           cttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCg                          <  1:30156302/68‑1 (MQ=14)
                                            ttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGt                         >  1:27528156/1‑68 (MQ=255)
                                             tcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTg                        <  1:13224048/68‑1 (MQ=255)
                                               aTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCa                      >  1:2980002/1‑68 (MQ=255)
                                                tCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAg                     <  1:11956052/68‑1 (MQ=255)
                                                tCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAg                     <  1:15396262/68‑1 (MQ=255)
                                                 cccTGAAATGTTTGCAATGAAGAGTCCATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAgg                    <  1:37887340/68‑1 (MQ=255)
                                                 cccTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAgt                    <  1:34386117/68‑2 (MQ=255)
                                                 cccTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAgg                    <  1:5476527/68‑1 (MQ=255)
                                                 cccTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAgg                    <  1:4805943/68‑1 (MQ=255)
                                                 cccTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAgg                    <  1:2005823/68‑1 (MQ=255)
                                                 cccTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAgg                    <  1:10577601/68‑1 (MQ=255)
                                                 cccTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAgg                    <  1:3692578/68‑1 (MQ=255)
                                                 cccTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAgg                    <  1:38396892/68‑1 (MQ=255)
                                                  ccTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAggg                   >  1:39097417/1‑68 (MQ=255)
                                                  ccTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAggg                   <  1:22857936/68‑1 (MQ=255)
                                                  ccTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAggg                   >  1:22496451/1‑68 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:38149776/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:8032011/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:9038095/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:11285075/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:4236729/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:10025767/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:25163298/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:18691875/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:2365690/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:21594412/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:32769292/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGAACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:37216020/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATCAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:4400853/68‑1 (MQ=255)
                                                      aaaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGtgtg               <  1:18537475/68‑1 (MQ=255)
                                                       aaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGtgtg               <  1:8381136/67‑1 (MQ=255)
                                                       aaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGtgtg               <  1:19256532/67‑1 (MQ=255)
                                                       aaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTgg              >  1:37615936/1‑68 (MQ=255)
                                                       aaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTgg              >  1:25818232/1‑68 (MQ=255)
                                                       aaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTgg              >  1:38752213/1‑68 (MQ=255)
                                                       aaTGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGt                         >  1:2459463/1‑57 (MQ=16)
                                                        aTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGt             <  1:20312075/68‑1 (MQ=255)
                                                        aTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGt             <  1:27305156/68‑1 (MQ=255)
                                                        aTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGt             <  1:12416031/68‑1 (MQ=255)
                                                        aTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGt             <  1:6667899/68‑1 (MQ=255)
                                                        aTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGt             <  1:18228757/68‑1 (MQ=255)
                                                         tGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGt             >  1:14057947/1‑67 (MQ=255)
                                                         tGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGt             >  1:34984444/1‑67 (MQ=255)
                                                          gTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTAc           >  1:12800536/1‑68 (MQ=255)
                                                             tGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAgg        <  1:22087509/68‑1 (MQ=255)
                                                                aaTTAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTc     <  1:26362579/68‑1 (MQ=255)
                                                                aaTGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTc     <  1:28111620/68‑1 (MQ=255)
                                                                aaTGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTc     <  1:24800293/68‑1 (MQ=255)
                                                                 aTGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTc     >  1:28465863/1‑67 (MQ=255)
                                                                 aTGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTc     >  1:26843140/1‑67 (MQ=255)
                                                                   gAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCa  <  1:30795631/68‑1 (MQ=255)
                                                                   |                                                                   
TAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCA  >  CP009273/1265972‑1266106

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