Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,575 C→G 17.6% A41P (GCC→CCC) ‡ argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,5750CG17.6% 148.7 / 25.4 68A41P (GCC→CCC) ‡argIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base C (21/35);  new base G (3/9);  total (24/44)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTC  >  minE/2919506‑2919641
                                                                     |                                                                  
cGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGct                                                                   >  1:3577608/1‑71 (MQ=255)
   aaGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttct                                                                >  1:1460845/1‑71 (MQ=255)
     gagCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttc                                                              <  1:2143960/71‑1 (MQ=255)
     gagCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttc                                                              <  1:1332095/71‑1 (MQ=255)
      agCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttct                                                             >  1:2330690/1‑71 (MQ=255)
      agCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttct                                                             <  1:3049360/71‑1 (MQ=255)
       gCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttctt                                                            >  1:2547126/1‑71 (MQ=255)
        cATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttctt                                                            <  1:2522935/70‑1 (MQ=255)
               tACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCtt                                                    <  1:3063819/71‑1 (MQ=255)
               tACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCtt                                                    <  1:2555359/71‑1 (MQ=255)
                 cgagtcgagtCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCttt                                                   >  1:1144197/1‑70 (MQ=255)
                 cgagtcgagtCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCttt                                                   <  1:1964499/70‑1 (MQ=255)
                  gagtcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttctt                                                            <  1:2133161/60‑1 (MQ=255)
                    gtcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCtt                                               <  1:3491272/71‑1 (MQ=255)
                    gtcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCtt                                               <  1:156901/71‑1 (MQ=255)
                    gtcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCtt                                               <  1:458861/71‑1 (MQ=255)
                    gtcgagtcTTTTTCGAAGATGAGCCCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCt                                                >  1:1430430/1‑70 (MQ=255)
                     tcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTa                                              <  1:2460278/71‑1 (MQ=255)
                      cgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCtt                                                    <  1:932425/64‑1 (MQ=255)
                       gagtcTTTTTGGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATc                                            >  1:2406720/1‑71 (MQ=255)
                       gagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGATTCTGTACTTCCTTGCCATTTTTTTTATc                                            >  1:1382560/1‑71 (MQ=255)
                       gagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTAt                                             <  1:1185952/70‑1 (MQ=255)
                       gagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATc                                            <  1:917204/71‑1 (MQ=255)
                         gtcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCg                                           >  1:737502/1‑70 (MQ=255)
                         gtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCgg                                          <  1:1258608/71‑1 (MQ=255)
                          tcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGc                                         <  1:3193547/71‑1 (MQ=255)
                          tcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttc                                                              >  1:3147016/1‑50 (MQ=255)
                            tttttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttctt                                                            <  1:1690350/50‑1 (MQ=255)
                            tttttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttct                                                                <  1:2084985/46‑1 (MQ=255)
                            tttttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTc                                                 >  1:1812235/1‑61 (MQ=255)
                              tttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCtt                                       <  1:3306119/69‑1 (MQ=255)
                              tttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCtt                                       <  1:2106005/69‑1 (MQ=255)
                              tttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCt                                        <  1:260631/68‑1 (MQ=255)
                              tttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCt                                        <  1:2574720/68‑1 (MQ=255)
                              tttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCt                                        <  1:3732892/68‑1 (MQ=255)
                              tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGct                                                                   <  1:1047821/41‑1 (MQ=255)
                              tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCttt                                      <  1:1167246/70‑1 (MQ=255)
                              tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCttt                                      <  1:1762441/70‑1 (MQ=255)
                               ttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCtttg                                     <  1:875383/70‑2 (MQ=255)
                               ttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTc                                     <  1:918403/70‑1 (MQ=255)
                                tCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGct                                                                   <  1:963369/39‑1 (MQ=255)
                                  gAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCttt                                                   <  1:2989852/53‑1 (MQ=255)
                                  gAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATc                                            >  1:2414015/1‑60 (MQ=255)
                                    aGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTa                                                          >  1:996348/1‑44 (MQ=255)
                                    aGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTa                                                          >  1:1467276/1‑44 (MQ=255)
                                     gATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttc                                                              >  1:3382425/1‑39 (MQ=255)
                                     gATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttct                                                                >  1:3680665/1‑37 (MQ=255)
                                     gATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTc                                                 >  1:3413156/1‑52 (MQ=255)
                                     gATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTcg                              >  1:1783812/1‑71 (MQ=255)
                                     gATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTcg                              <  1:1160/71‑1 (MQ=255)
                                      aTGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                                       <  1:3525389/61‑1 (MQ=255)
                                      aTGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTcgc                             >  1:104555/1‑71 (MQ=255)
                                      aTGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTcgc                             <  1:1868592/71‑1 (MQ=255)
                                       tGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGc                                         <  1:1810735/58‑1 (MQ=255)
                                       tGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCtttt                                                  >  1:3398013/1‑49 (MQ=255)
                                        gAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCttt                                                   >  1:2509847/1‑47 (MQ=255)
                                        gAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCttt                                                   >  1:79207/1‑47 (MQ=255)
                                         aGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGc                          >  1:89350/1‑71 (MQ=255)
                                          gcgcGATGTTTTTACCAGTGTGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCt                         <  1:2880690/71‑1 (MQ=255)
                                          gcgcGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCt                         <  1:2324083/71‑1 (MQ=255)
                                          gcgcGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCt                         <  1:2816680/71‑1 (MQ=255)
                                           cgcgATGTTTTTAGCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTa                        <  1:1296647/71‑1 (MQ=255)
                                           cgcgATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGcaa                        <  1:3578651/71‑3 (MQ=255)
                                            gcgATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTa                                                          <  1:307647/36‑1 (MQ=255)
                                             cgATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTc                                     >  1:726269/1‑56 (MQ=255)
                                               aTGTTTTTACCAGTGAGGTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTa                        <  1:728554/67‑1 (MQ=255)
                                                tGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCg                                           >  1:194072/1‑47 (MQ=255)
                                                tGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCg                                           >  1:3055440/1‑47 (MQ=255)
                                                 gTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCttt                                      >  1:1636867/1‑51 (MQ=255)
                                                  tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCt                                                <  1:2845223/40‑1 (MQ=255)
                                                  tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCt                                                <  1:792723/40‑1 (MQ=255)
                                                  tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCt                                                <  1:451192/40‑1 (MQ=255)
                                                  tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCa                                    <  1:1537903/52‑1 (MQ=255)
                                                  tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcag                 <  1:2830672/71‑1 (MQ=255)
                                                  tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgca                  <  1:2776959/70‑1 (MQ=255)
                                                    tttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcag                 >  1:455139/1‑69 (MQ=255)
                                                         cAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTg          <  1:1053806/71‑1 (MQ=255)
                                                           gTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGtt        <  1:857016/71‑1 (MQ=255)
                                                           gTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGtt        <  1:2841737/71‑1 (MQ=255)
                                                            tGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGtt        >  1:1997174/1‑70 (MQ=255)
                                                            tGAGTTTGGCTTCTTCTTTACCGCTTTTATTATCGGCTTTCAGCTTCGCGGCt                         >  1:2001694/1‑53 (MQ=255)
                                                             gAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGtt        >  1:3628590/1‑69 (MQ=255)
                                                              aGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTg       <  1:683609/69‑1 (MQ=255)
                                                               gTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTct                              <  1:111335/45‑2 (MQ=255)
                                                                     cttcttcttTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTc  <  1:1127622/67‑1 (MQ=255)
                                                                     |                                                                  
CGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTC  >  minE/2919506‑2919641

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: