Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 68,344 A→C 15.5% intergenic (+32/+143) PP_0057 → / ← PP_0058 major facilitator family transporter/1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294768,3440AC15.5% 135.7 / 18.6 58intergenic (+32/+143)PP_0057/PP_0058major facilitator family transporter/1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (25/24);  new base C (4/5);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

TACCGGGGATATCTTCTACGGGCTGTGGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACAAGTTATCCACAGAAGGTCAAGCC  >  NC_002947/68198‑68493
                                                                                                                                                  |                                                                                                                                                      
taCCGGGGATATCTTCTACGGGCTGTGGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGc                                                                                                                                                     <  1:1077309/150‑1 (MQ=255)
  ccGGGGATATCTTCTACGGGCTGTGGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTGGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCg                                                                                                                                                   <  1:452123/150‑1 (MQ=255)
   cGGGGATATCTTCTACGGGCTGTGGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCgg                                                                                                                                                  <  1:354695/150‑1 (MQ=255)
           tcttctACGGGCTGTGGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTa                                                                                                                                          <  1:849383/150‑1 (MQ=255)
            cttctACGGGCTGTGGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTTTTACAGCCGGTTTTTTTAt                                                                                                                                         <  2:1032361/150‑1 (MQ=255)
             ttctACGGGCTGTGGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATg                                                                                                                                        >  1:1068358/1‑150 (MQ=255)
                     gCTGTGGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCaaaa                                                                                                                               >  2:515027/1‑150 (MQ=255)
                            tACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTa                                                                                                                        >  1:714595/1‑150 (MQ=255)
                                 ggTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCAc                                                                                                                   <  1:515027/150‑1 (MQ=255)
                                             cGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTTTTACAGCCGGttttt                                                                                                                                             >  2:967100/1‑113 (MQ=39)
                                             cGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTTTTACAGCCGGttttt                                                                                                                                             <  1:967100/113‑1 (MQ=39)
                                             cGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGaaa                                                                                                       >  2:1125748/1‑150 (MQ=255)
                                             cGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGaaa                                                                                                       >  2:270476/1‑150 (MQ=255)
                                                    aGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCAATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCa                                                                                                <  2:94802/150‑1 (MQ=255)
                                                            ggtTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGa                                                                                        <  1:828972/150‑1 (MQ=255)
                                                              ttGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGaaa                                                                                      <  2:1068358/150‑1 (MQ=255)
                                                                  gATCTTCTGTCTGAGAGAAACCAAGGATTTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAAc                                                                                  >  1:1032544/1‑150 (MQ=255)
                                                                  gATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAAc                                                                                  <  2:714595/150‑1 (MQ=255)
                                                                           tctgAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAAtt                                                                         >  1:963587/1‑150 (MQ=255)
                                                                             tgAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCa                                                                       >  2:482523/1‑150 (MQ=255)
                                                                              gagagaAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCAt                                                                      >  2:114681/1‑150 (MQ=255)
                                                                               agagaAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATaa                                                                                   <  1:951177/136‑1 (MQ=255)
                                                                               agagaAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATaa                                                                                   >  2:951177/1‑136 (MQ=255)
                                                                                gagaAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTc                                                                        <  1:925723/146‑1 (MQ=255)
                                                                                gagaAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTc                                                                        >  2:925723/1‑146 (MQ=255)
                                                                                 agaAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACCTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATaca                                                                   <  1:733796/150‑1 (MQ=255)
                                                                                 agaAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACCTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATaca                                                                   >  2:733796/1‑150 (MQ=255)
                                                                                                    atCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGc                                                >  2:1179956/1‑150 (MQ=255)
                                                                                                    atCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGc                                                >  2:1179955/1‑150 (MQ=255)
                                                                                                    atCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCTCAAATGAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGc                                                <  1:927956/150‑1 (MQ=255)
                                                                                                     tCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCaa                                                  >  1:1080197/1‑147 (MQ=255)
                                                                                                     tCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCaa                                                  <  2:1080197/147‑1 (MQ=255)
                                                                                                        aCAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCAttt                                            >  2:340220/1‑150 (MQ=255)
                                                                                                          aaaGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTaa                                          >  2:176018/1‑150 (MQ=255)
                                                                                                          aaaGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTaa                                          <  1:114681/150‑1 (MQ=255)
                                                                                                           aaGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCAAAAAAAAAAAAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAAt                                         <  2:997223/150‑1 (MQ=255)
                                                                                                              tCTGATGTCCACATGCCATAAAAAAACCGGCTTTTACAGCCGGttt                                                                                                                                               <  1:810250/46‑1 (MQ=14)
                                                                                                              tCTGATGTCCACATGCCATAAAAAAACCGGCTTTTACAGCCGGttt                                                                                                                                               >  2:810250/1‑46 (MQ=14)
                                                                                                                tGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCaa                                     >  1:212183/1‑149 (MQ=255)
                                                                                                                tGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCaa                                     <  2:212183/149‑1 (MQ=255)
                                                                                                                  aTGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAg                                                 <  1:414533/135‑1 (MQ=255)
                                                                                                                  aTGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAg                                                 >  2:414533/1‑135 (MQ=255)
                                                                                                                        acaTGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGcac                            >  1:892648/1‑150 (MQ=255)
                                                                                                                            gCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACaa                        <  1:1125748/150‑1 (MQ=255)
                                                                                                                             gcATAAAAAAACCGGCTTTTACAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGatat                                                          >  1:596979/2‑115 (MQ=255)
                                                                                                                             gcATAAAAAAACCGGCTTTTACAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGatat                                                          <  2:596979/114‑1 (MQ=255)
                                                                                                                               aTAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTc                                                    >  1:900368/1‑119 (MQ=255)
                                                                                                                               aTAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTc                                                    <  2:900368/119‑1 (MQ=255)
                                                                                                                                  aaaaaaCCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACAAGTTATc                  <  2:1032544/150‑1 (MQ=255)
                                                                                                                                   aaaaaCCGGCTTTTACAGCCGGTTTTTTTATGGCTCAAAAAAACTTATGCATGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAAtt                                                                         <  1:549764/94‑1 (MQ=37)
                                                                                                                                   aaaaaCCGGCTTTTACAGCCGGTTTTTTTATGGCTCAAAAAAACTTATGCATGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAAtt                                                                         >  2:549764/1‑94 (MQ=25)
                                                                                                                                      aaCCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTc                                                                        >  2:147436/1‑92 (MQ=255)
                                                                                                                                      aaCCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTc                                                                        <  1:147436/92‑1 (MQ=255)
                                                                                                                                       aCCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACAAGTTATCCACAg             <  1:340220/150‑1 (MQ=255)
                                                                                                                                             tGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCCAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACAAGTTATCCACAGAAGGTc       >  2:342940/1‑150 (MQ=255)
                                                                                                                                             tGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACAAGTTATCCACAGAAGGTc       >  2:1141377/1‑150 (MQ=255)
                                                                                                                                               tAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACAAGTTATCcaca              <  1:451838/141‑1 (MQ=255)
                                                                                                                                               tAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACAAGTTATCcaca              >  2:451838/1‑141 (MQ=255)
                                                                                                                                                  aaGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACAAGTTATCCACAGAAGGTCAAGcc  >  1:197491/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
TACCGGGGATATCTTCTACGGGCTGTGGTACCCGGTGCTGGTGACCGGGATCAGCCTGGTGGTTGGGATCTTCTGTCTGAGAGAAACCAAGGATGTAGATATCGACAAAGTCTGATGTCCACATGCCATAAAAAAACCGGCTGTAAAAGCCGGTTTTTTTATGCCTCAAAAAAACTTATGCACGATTTTCAGAAAAAATTCATACAAAGAAATAACTAGCAAATTCATACACTTAGATATCAATTCAAGCATTTAATCCAAAAACTGCACACAAGTTATCCACAGAAGGTCAAGCC  >  NC_002947/68198‑68493

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