Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,266,039 G→A 6.4% 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 6.4% 1927.1 / 72.4 616intergenic (+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 (195/362);  new base A (28/10);  total (225/391)
Fisher's exact test for biased strand distribution p-value = 3.31e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 627 total aligned reads displayed.

TAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCA  >  CP009273/1265972‑1266106
                                                                   |                                                                   
tAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATg                                                                     >  1:23694253/1‑68 (MQ=40)
 aaCGCCGTCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGa                                                                    <  1:18997048/68‑1 (MQ=255)
 aaCGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGa                                                                    <  1:11396179/68‑1 (MQ=40)
 aaCGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGa                                                                    <  1:5389687/68‑1 (MQ=40)
  aCGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGaa                                                                   >  1:19393411/1‑68 (MQ=40)
  aCGCCGGCGTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCCTCCCTGAAATGTTTGCAATGaa                                                                   >  1:22609034/1‑68 (MQ=255)
   cGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAag                                                                  <  1:931208/68‑1 (MQ=40)
    gCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAaga                                                                 >  1:21471569/1‑68 (MQ=40)
    gCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAaga                                                                 <  1:20832984/68‑1 (MQ=40)
    gCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAaga                                                                 >  1:15263624/1‑68 (MQ=40)
      cGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAagag                                                                <  1:7736184/67‑1 (MQ=40)
      cGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAagag                                                                <  1:27961196/67‑1 (MQ=40)
      cGGCTTGGTGTGGATACTACGTCGAAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGt                                                               <  1:22095483/68‑1 (MQ=255)
               gtgGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCg                                                       >  1:11660704/1‑67 (MQ=40)
                 ggATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCt                                                    <  1:18733898/68‑1 (MQ=40)
                  gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCtt                                                   >  1:2839341/1‑68 (MQ=40)
                   aTACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCttt                                                  >  1:16643859/1‑68 (MQ=40)
                    tactacGTCGCAATTCATCTTCCCTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCtttt                                                 <  1:20089763/68‑1 (MQ=255)
                    tactacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCtttt                                                 <  1:27664073/68‑1 (MQ=40)
                    tactacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCttt                                                  >  1:29832710/1‑67 (MQ=40)
                    tactacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCttt                                                  >  1:25889594/1‑67 (MQ=40)
                    tactacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCttt                                                  >  1:21762006/1‑67 (MQ=40)
                    tactacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCttt                                                  >  1:11953688/1‑67 (MQ=40)
                    tactacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCttt                                                  >  1:16063269/1‑67 (MQ=40)
                    tactacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCttt                                                  >  1:18306243/1‑67 (MQ=40)
                    tactacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCttt                                                  >  1:21232410/1‑67 (MQ=40)
                    tactacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCttt                                                  >  1:17534609/1‑67 (MQ=40)
                     actacGTCGCAATTCATCTTCACTTCCTCACTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTc                                                >  1:5218448/1‑68 (MQ=255)
                      ctacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCa                                               >  1:2651553/1‑68 (MQ=40)
                      ctacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCa                                               >  1:5678402/1‑68 (MQ=40)
                      ctacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCa                                               <  1:25220050/68‑1 (MQ=40)
                        acGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAAc                                             >  1:14317801/1‑68 (MQ=40)
                        acGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAAc                                             >  1:17800297/1‑68 (MQ=40)
                        acGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAAc                                             >  1:25637212/1‑68 (MQ=40)
                        acGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAAc                                             >  1:2888172/1‑68 (MQ=40)
                        acGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAAc                                             >  1:5876441/1‑68 (MQ=40)
                             gCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACa                                            >  1:13845506/1‑64 (MQ=37)
                             gCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACa                                            >  1:14636288/1‑64 (MQ=37)
                             gCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACa                                            >  1:6074132/1‑64 (MQ=37)
                               aaTTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGtt                                      <  1:3883905/68‑1 (MQ=35)
                               aaTTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGtt                                      <  1:26909746/68‑1 (MQ=35)
                                  tCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCa                                   <  1:864671/68‑1 (MQ=35)
                                  tCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCa                                   <  1:11118226/68‑1 (MQ=35)
                                  tCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCa                                   <  1:1481079/68‑1 (MQ=35)
                                  tCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCa                                   <  1:14059402/68‑1 (MQ=35)
                                     tCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTg                                >  1:22418493/1‑68 (MQ=37)
                                       ttcacttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGt                              <  1:24057529/68‑1 (MQ=40)
                                       ttcacttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGt                              <  1:10049291/68‑1 (MQ=40)
                                       ttcacttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGt                              <  1:19687486/68‑1 (MQ=40)
                                       ttcacttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGt                              <  1:16285599/68‑1 (MQ=40)
                                       ttcacttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGt                              <  1:9014841/68‑1 (MQ=40)
                                         cacttcaTCCCTGAAATGTTTGCATTTAAGAGTACATTCCGTCTTTTCAACAGCTGTTGCAGTGGttt                            <  1:8312968/68‑1 (MQ=11)
                                         cacttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGttt                            <  1:1816045/68‑1 (MQ=40)
                                         cacttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGttt                            <  1:15734267/68‑1 (MQ=40)
                                          acttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTc                           >  1:27061366/1‑68 (MQ=40)
                                           cttcaTCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCg                          <  1:2972920/68‑1 (MQ=255)
                                            ttcaTCCCTGAAATGTTTGCAATTAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGt                         >  1:16402438/1‑68 (MQ=255)
                                             tcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGt                         >  1:5539982/1‑67 (MQ=18)
                                               aTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCa                      >  1:18432149/1‑68 (MQ=255)
                                                 cccTGAACTCTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAgg                    <  1:22862579/68‑1 (MQ=14)
                                                 cccTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAgg                    <  1:26748292/68‑1 (MQ=255)
                                                 cccTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAgg                    <  1:12140315/68‑1 (MQ=255)
                                                 cccTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAgg                    <  1:6688568/68‑1 (MQ=255)
                                                 cccTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAgg                    <  1:24303019/68‑1 (MQ=255)
                                                 cccTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAgg                    <  1:15894489/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCTATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:9758234/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATTAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:25033525/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:17857614/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:12754663/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:7260513/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:12989810/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:13339958/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:15019579/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:15510672/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:17085571/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:4900773/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:17970160/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:3197346/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:1928514/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:29498881/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:27422328/68‑1 (MQ=255)
                                                   cTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgt                  <  1:2821315/68‑1 (MQ=255)
                                                     gAAATGTTTGCAATGAAGATTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgtgt                <  1:4245306/68‑1 (MQ=255)
                                                     gAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgtgt                <  1:6961472/68‑1 (MQ=255)
                                                     gAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGgtgt                <  1:10588970/68‑1 (MQ=255)
                                                       aacgTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTgg              >  1:24845429/4‑68 (MQ=255)
                                                       aaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACATCTGTTGCAGTGGTTTCGTGCAGGGTGTgg              <  1:10874965/68‑1 (MQ=255)
                                                       aaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGtgtg               <  1:3678941/67‑1 (MQ=255)
                                                       aaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGtgtg               <  1:11791781/67‑1 (MQ=255)
                                                       aaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGtgtg               <  1:14131144/67‑1 (MQ=255)
                                                       aaTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGtgtg               <  1:23233507/67‑1 (MQ=255)
                                                       aaTGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTgg              >  1:6409429/1‑68 (MQ=255)
                                                       aaTGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTgg              >  1:20646544/1‑68 (MQ=255)
                                                       aaTGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTgg              >  1:20344780/1‑68 (MQ=255)
                                                        aTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGt             <  1:8128928/68‑1 (MQ=255)
                                                        aTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGt             <  1:26225739/68‑1 (MQ=255)
                                                        aTGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGt             <  1:10288046/68‑1 (MQ=255)
                                                                 aTGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTc     >  1:23441882/1‑67 (MQ=255)
                                                                 aTGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTc     >  1:4612219/1‑67 (MQ=255)
                                                                   gAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCa  >  1:1227528/1‑68 (MQ=255)
                                                                   |                                                                   
TAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTGGTTTCGTGCAGGGTGTGGTACAGGCTCGCA  >  CP009273/1265972‑1266106

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