Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,302,995 T→G 12.4% intergenic (‑93/‑50) clsA ← / → yciY cardiolipin synthase 1/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,302,9950TG12.4% 2001.8 / 52.3 800intergenic (‑93/‑50)clsA/yciYcardiolipin synthase 1/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (379/303);  new base G (2/113);  total (383/417)
Fisher's exact test for biased strand distribution p-value = 1.01e-32
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 809 total aligned reads displayed.

GTTATGAGTTTACGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGAGTAGAACG  >  CP009273/1302931‑1303062
                                                                |                                                                   
gTTATGAGTTTACGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCgg                                                                  >  1:25537442/1‑68 (MQ=255)
gTTATGAGTTTACGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCgg                                                                  <  1:35050678/68‑1 (MQ=255)
  tATGAGTTTACGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCgggg                                                                <  1:28725529/68‑1 (MQ=255)
      aGTTTACGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCgggggtgg                                                            >  1:16910115/1‑68 (MQ=255)
      aGTTTACGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCgggggtgg                                                            >  1:19210602/1‑68 (MQ=255)
       gTTTACGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCgggggtggg                                                           >  1:49790569/1‑68 (MQ=255)
       gTTTACGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCgggggtggg                                                           >  1:44736735/1‑68 (MQ=255)
       gTTTACGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCgggggtggg                                                           >  1:48779337/1‑68 (MQ=255)
       gTTTACGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCgggggtggg                                                           >  1:47733689/1‑68 (MQ=255)
       gTTTACGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCgggggtggg                                                           >  1:32141503/1‑68 (MQ=255)
       gTTTACGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCgggggtggg                                                           >  1:3739963/1‑68 (MQ=255)
       gTTTACGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCgggggtggg                                                           >  1:12543560/1‑68 (MQ=255)
       gTTTACGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCgggggtggg                                                           >  1:12221939/1‑68 (MQ=255)
       gTTTACGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCgggggtggg                                                           >  1:19837070/1‑68 (MQ=255)
       gTTTACGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCgggggtggg                                                           >  1:50689898/1‑68 (MQ=255)
       gTTTACGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCgggggtggg                                                           >  1:37704940/1‑68 (MQ=255)
       gTTTACGCATATGGATTCATCTGGCAAATAAAAACGCGAGAAAAGCACTGGTCAGCGTCgggggtggg                                                           >  1:40951529/1‑68 (MQ=255)
       gTTTACGCATATGGAATCATCTGGCAAATAAAAACGCGGGAAAAGAACTGGTCAGCGTCgggggtggg                                                           >  1:10977821/1‑68 (MQ=255)
         ttAGGCATATGGATTCATCTGGCAAATAAAAAGGCGGGAAAAGCACTGGTCAGCGTCgggggtgggg                                                          <  1:29801966/67‑1 (MQ=255)
         ttAGGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCgggggtgggg                                                          <  1:18856786/67‑1 (MQ=255)
         ttACGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCgggggtggg                                                           <  1:12974165/66‑1 (MQ=255)
         ttACGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCgggggtggg                                                           <  1:44142356/66‑1 (MQ=255)
         ttACGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGGCgggggtggg                                                           <  1:11896028/66‑1 (MQ=255)
           aCGCATATTGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGta                                                       <  1:24459834/68‑1 (MQ=255)
           aCGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGta                                                       <  1:34124817/68‑1 (MQ=255)
            cGCATTTGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCTTCGGGGGGGGGGGtat                                                      <  1:50502299/68‑1 (MQ=255)
            cGCATATGGATTCATCTGTCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGcat                                                      <  1:23660551/68‑1 (MQ=255)
            cGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGtat                                                      <  1:36002677/68‑1 (MQ=255)
            cGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGtat                                                      >  1:51226536/1‑68 (MQ=255)
            cGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGtat                                                      <  1:5466894/68‑1 (MQ=255)
            cGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGtat                                                      >  1:43657393/1‑68 (MQ=255)
            cGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGtat                                                      >  1:13459721/1‑68 (MQ=255)
            cGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGtat                                                      <  1:8703881/68‑1 (MQ=255)
             gCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGtata                                                     >  1:2334153/1‑68 (MQ=255)
             gCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGtata                                                     >  1:6053947/1‑68 (MQ=255)
             gCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGtata                                                     >  1:50987102/1‑68 (MQ=255)
              cATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGtata                                                     >  1:7960802/1‑67 (MQ=255)
              cATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATaa                                                    >  1:26097075/1‑68 (MQ=255)
              cATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATaa                                                    >  1:31620567/1‑68 (MQ=255)
               atatTGATTCTTCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATaa                                                    <  1:28060626/67‑1 (MQ=255)
                  tGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGa                                                 >  1:33517979/1‑67 (MQ=255)
                  tGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGa                                                 >  1:39863918/1‑67 (MQ=255)
                   ggATTCATCTGGCAAATAAAAACGCGGGAAAAGCTCTGGTCAGCGTCGGGGGTGGGGGTATAATGAcc                                               <  1:24172758/68‑1 (MQ=255)
                   ggATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGAcc                                               >  1:15015860/1‑68 (MQ=255)
                   ggATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGAcc                                               <  1:22999234/68‑1 (MQ=255)
                     aTTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCAt                                             >  1:15365687/1‑68 (MQ=255)
                     aTTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCAt                                             >  1:14704889/1‑68 (MQ=255)
                     aTTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCAt                                             >  1:41883172/1‑68 (MQ=255)
                     aTTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCAt                                             >  1:7286042/1‑68 (MQ=255)
                     aTTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCAt                                             >  1:34613362/1‑68 (MQ=255)
                     aTTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCAt                                             >  1:22067260/1‑68 (MQ=255)
                      ttCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCAtt                                            <  1:9146917/68‑1 (MQ=255)
                       tCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCAt                                             >  1:1641493/1‑66 (MQ=255)
                       tCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCAt                                             >  1:39291995/1‑66 (MQ=255)
                          tCTGGCAAATAAAAACGCGGGAAACGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCATTCTGt                                        <  1:40143008/68‑1 (MQ=255)
                          tCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCATTCTGt                                        >  1:11587333/1‑68 (MQ=255)
                          tCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCATTATGt                                        >  1:32626581/1‑68 (MQ=255)
                           cTGGCAACAAAAAACGCCGGACCCGCACTGGTCCGCGGCGGGGGGGGGCGTATAATGACCATTCTGtt                                       <  1:36473061/68‑1 (MQ=255)
                           cTGGCAAATAAAAACGCGGGAAAAGCACTGGTCCGCGTCGGGGGGGGGGGTATAATGACCATTCTGtt                                       <  1:31194473/68‑1 (MQ=255)
                           cTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCATTCTGtt                                       >  1:29192462/1‑68 (MQ=255)
                           cTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGGGGGCGTATAATGACCATTCTGtt                                       <  1:21375721/68‑1 (MQ=255)
                           cTGGCAAATAAAAACGCGGGAAAAGCACTGGGCCGCGGCGGGGGTGGGGGTATAATGACCATTCTGtt                                       <  1:6405947/68‑1 (MQ=255)
                           cTGGCAAATAAAAACGCGGGAAAAGCACGGGGCAGCGTCGGGGGTGGGGGTATAATGACCATTCTGtt                                       <  1:25845073/68‑1 (MQ=255)
                           cTGCCAAATAAAACCGCGGGCAAAGCCCAGGGCAGCGTCGGGGGTGGGCGTATAATGACCATTCTGtt                                       <  1:4735494/68‑1 (MQ=255)
                                aaaTAAAACCGCGGGAAAAGCACTGGTCAGCGGCGGGGGTGGGGGTATAATGACCATTCTGTTATTGc                                  <  1:17237860/68‑1 (MQ=255)
                                aaaTAAAAACGCGTGAAAAGCCCTGGGCAGCGGCGGGGGTGGGGGTATAATGACCATTCTGTTATTGc                                  <  1:37225865/68‑1 (MQ=255)
                                aaaTAAAAACGCGGTAACAGCACTGGTCAGCGTCGGGGGGGGGGGTATAATGACCATTCTGTTATTGc                                  <  1:37901175/68‑1 (MQ=255)
                                aaaTAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGCTATAATGACCATTCTGTTATTGc                                  <  1:41401725/68‑1 (MQ=255)
                                aaaTAAAAACGCGGGAAAAGCACTGGGCTGCGTCTGGGGTGGGGGTATAATGACCATTCTGTTATTGc                                  <  1:13834999/68‑1 (MQ=255)
                                     taaaCGCGGGAAAAGCCCTGGTCTGCGTCGGGGGGGGGGGTATAATGACCATTCTGTTATTGCATaga                             <  1:10341450/67‑1 (MQ=38)
                                     cacaCGCGGGAAAAGCACTGGCCCGCGGCGGGGGGGGGGGTATAATGACCATTCTGTTATTGCATaga                             <  1:6974676/65‑1 (MQ=255)
                                     aaacCGCGGGAAAAGGCCTGGTCCGCGTCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATaga                             <  1:19602792/68‑1 (MQ=255)
                                     aaacCGCGGCAAAAGCACTGGGCCGCGGCGGGGGGGGGGGTATAATGACCATTCTGTTATTGCATaga                             <  1:45742921/68‑1 (MQ=255)
                                     aaaaCGCGGTAAAACCCCTGGTCAGCGGGGGGGGGGGGGGAATAATGACCATTCTGTTATTGCATaga                             <  1:28326879/68‑1 (MQ=255)
                                     aaaaCGCGGGAAACGCACTGGCCCGCGGGGGGGGTGGGGCTATAATGACCATTCTGTTATTGCATaga                             <  1:27416498/68‑1 (MQ=255)
                                     aaaaCGCGGGAAAAGCTCTGGTCAGCGTCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATaga                             <  1:17694335/68‑1 (MQ=255)
                                     aaaaCGCGGGAAAAGCACTGGTCCGCGTCGGCGGGGGGGGTATAATGACCATTCTGTTATTGCATaga                             <  1:40414794/68‑1 (MQ=255)
                                     aaaaCGCGGGAAAAGCACTGGTCAGCGTCGGGGGGGGGGGTATAATGACCATTCTGTTATTGCATaga                             <  1:7672821/68‑1 (MQ=255)
                                     aaaaCGCGGGAAAAGCACTGGTCAGCGGCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATaga                             <  1:48220506/68‑1 (MQ=255)
                                     aaaaCGCGGGAAAACCACTGGTCAGCGGCGGGGGGGGGGGTATAATGACCATTCTGTTATTGCATaga                             <  1:43393649/68‑1 (MQ=255)
                                     aaaaCGCGCTAAAAGCACTGGGCAGCGTCGGGGGGGGGGGGATAATGACCATTCTGTTATTGCATaga                             <  1:2492036/68‑1 (MQ=255)
                                         cgctGGAATAGCACTGGTCCGCGTCGGGGGTGGTGGTATAATGACCATTCTGTTATTGCATAGagtag                         <  1:6744244/68‑1 (MQ=255)
                                         cgcgGGAAAAGCACTGGTCAGCGGCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATAGagtag                         <  1:15994806/68‑1 (MQ=255)
                                          gcgGTAAAAGGCCTGGTCAGCGTCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATAGagtagt                        <  1:20641701/68‑1 (MQ=255)
                                          gcgGGAAAAGCACTGGTCAGAGACCGGCGTGCGGGTATAAAGACAATTCTCATAATGCATAGagtagt                        >  1:10008748/1‑68 (MQ=255)
                                                   gcccTGGGCAGCGTCGGGGGGGGGGGTATAAAGACCATTCTGTTATTGCATAGAGTAGTTAACATGaa               <  1:42701874/65‑1 (MQ=255)
                                                     aCTGGTCAGCGTCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGc             >  1:18394213/1‑68 (MQ=255)
                                                     aCTGGTCAGCGTCGGGGGGGGGGGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGc             <  1:46714605/68‑1 (MQ=255)
                                                        gggcAGCGTGGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGa          <  1:48489535/65‑1 (MQ=255)
                                                        gggcAGCGGCGGGGGGGGGGGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGa          <  1:22342888/65‑1 (MQ=255)
                                                        ggTCAGCGTCGGGGGGGGGGGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGa          <  1:5807275/68‑1 (MQ=255)
                                                        ggTCAGCGTCGGGGGGGGGGGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGa          <  1:45143097/68‑1 (MQ=255)
                                                         ggcAGCGGCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGa          <  1:30709501/65‑1 (MQ=255)
                                                         gTCAGCGGCGGGGGCGGGGCTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGa          <  1:4623244/67‑1 (MQ=255)
                                                          tCAGCGTCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGAGt        >  1:26605236/1‑68 (MQ=255)
                                                          tCAGCGTCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGAGt        >  1:3046711/1‑68 (MQ=255)
                                                           ccgCGTCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGAGt        <  1:4938942/65‑1 (MQ=255)
                                                           ccgCGTCGGGGGGGGGGGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGAGt        <  1:9420113/65‑1 (MQ=255)
                                                           catcGTCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGAGt        <  1:25011735/64‑1 (MQ=255)
                                                             gcggcgGGGGGGGGGGTAGAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGAGTAGa     <  1:38655486/68‑1 (MQ=255)
                                                                |                                                                   
GTTATGAGTTTACGCATATGGATTCATCTGGCAAATAAAAACGCGGGAAAAGCACTGGTCAGCGTCGGGGGTGGGGGTATAATGACCATTCTGTTATTGCATAGAGTAGTTAACATGAAGCGGAGTAGAACG  >  CP009273/1302931‑1303062

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