Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 980,236 C→A 11.2% P23P (CCG→CCT)  cspI ← cold shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE980,2360CA11.2% 331.2 / 16.4 108P23P (CCG→CCT) cspIcold shock protein
Reads supporting (aligned to +/- strand):  ref base C (45/50);  new base A (7/5);  total (53/55)
Fisher's exact test for biased strand distribution p-value = 5.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

Only 100 of 111 total aligned reads displayed.

TTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAA  >  minE/980169‑980306
                                                                   |                                                                      
ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                                                     <  1:1635304/71‑1 (MQ=255)
 tAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCg                                                                    >  1:3901107/1‑71 (MQ=255)
  aaTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                                                   >  1:2682005/1‑71 (MQ=255)
   aTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                                                   <  1:4822840/70‑1 (MQ=255)
      ttttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgat                                                                <  1:805262/70‑1 (MQ=255)
      ttttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgat                                                                <  1:2931050/70‑1 (MQ=255)
      ttttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgat                                                                <  1:2937666/70‑1 (MQ=255)
      ttttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgat                                                                <  1:5266302/70‑1 (MQ=255)
      ttttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                                                     <  1:1655403/65‑1 (MQ=255)
        ttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                                              <  1:4530798/70‑1 (MQ=255)
        ttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                                              <  1:3551699/70‑1 (MQ=255)
        ttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                                              <  1:2261872/70‑1 (MQ=255)
            aaTCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCa                                                         <  1:5691882/71‑1 (MQ=255)
             aTCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCa                                                         >  1:4960457/1‑70 (MQ=255)
             aTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                                                        <  1:4837128/71‑1 (MQ=255)
                   gCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCt                                                   <  1:3456392/70‑1 (MQ=255)
                    ctctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttt                                                 >  1:4420352/1‑71 (MQ=255)
                     tctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                                                   <  1:2962676/52‑1 (MQ=255)
                      ctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtga                                                                 >  1:4610993/1‑53 (MQ=255)
                      ctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtga                                                                 >  1:2937930/1‑53 (MQ=255)
                      ctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtga                                                                 >  1:767543/1‑53 (MQ=255)
                      ctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                               >  1:104868/1‑71 (MQ=255)
                       tGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTc                                                                        <  1:3852159/45‑1 (MQ=255)
                       tGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTc                                                                        <  1:3741904/45‑1 (MQ=255)
                       tGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                               >  1:128643/1‑70 (MQ=255)
                       tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                                                   <  1:1798200/50‑1 (MQ=255)
                       tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                                                   <  1:5022575/50‑1 (MQ=255)
                       tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                                                   <  1:1197165/50‑1 (MQ=255)
                        gAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaa                                                       >  1:4012266/1‑61 (MQ=255)
                         aaTTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                               <  1:3675894/68‑1 (MQ=255)
                         aaTTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                               <  1:1508089/68‑1 (MQ=255)
                         aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                                              <  1:2949205/69‑1 (MQ=255)
                         aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                            <  1:4446119/71‑1 (MQ=255)
                           ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatg                                                               >  1:2411473/1‑50 (MQ=255)
                           ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                             <  1:4578964/68‑1 (MQ=255)
                           ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                             <  1:5582488/68‑1 (MQ=255)
                           ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                             <  1:2919994/68‑1 (MQ=255)
                             gCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAc                                                     >  1:5680592/1‑58 (MQ=255)
                             gCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGt                                         >  1:3491141/1‑70 (MQ=255)
                             gCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGt                                         >  1:3022718/1‑70 (MQ=255)
                              cTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                               <  1:5457685/63‑1 (MQ=255)
                              cTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTa                                       <  1:3922659/71‑1 (MQ=255)
                              cTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                             >  1:2389092/1‑65 (MQ=255)
                                gagaAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTaaa                                     >  1:1216871/1‑71 (MQ=255)
                                 agaAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                                                   <  1:5135465/40‑1 (MQ=255)
                                  gaAATGGACAAACTCATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAc                                                     >  1:2788336/1‑53 (MQ=37)
                                  gaAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                                                   <  1:4193400/39‑1 (MQ=255)
                                   aaaTGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCa                                  <  1:2167906/71‑1 (MQ=255)
                                   aaaTAGACAAACACATCTTTGCTGCCATCTTTCGGCg                                                                    >  1:4403533/1‑37 (MQ=255)
                                   aaaTAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                    >  1:220450/1‑69 (MQ=255)
                                     atggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAtt                                >  1:1039078/4‑71 (MQ=255)
                                     atggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAtt                                >  1:1490823/4‑71 (MQ=255)
                                      tggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAttt                               <  1:1379447/69‑1 (MQ=255)
                                        gACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                            >  1:2381801/1‑56 (MQ=255)
                                        caCAAACACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCAtttta                             >  1:1453682/2‑69 (MQ=255)
                                        caCAAACACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACAAtttta                             >  1:1298922/2‑69 (MQ=255)
                                         aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaa                                                       >  1:3443239/1‑44 (MQ=255)
                                         aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaa                                                       >  1:4190745/1‑44 (MQ=255)
                                         aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                                              <  1:4785105/53‑1 (MQ=255)
                                         aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                                              <  1:4048922/53‑1 (MQ=255)
                                         aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                             <  1:395464/54‑1 (MQ=255)
                                           aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg                                           >  1:4306035/1‑54 (MQ=255)
                                           aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                          <  1:5510559/71‑1 (MQ=255)
                                           aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                          <  1:2315796/71‑1 (MQ=255)
                                           aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                          <  1:720516/71‑1 (MQ=255)
                                           aaaCACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCATTTTACt                           >  1:5153911/1‑70 (MQ=255)
                                             acacATCTTTAGAGCCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTAcga                          >  1:4621677/1‑67 (MQ=9)
                                             acacATCTTTAGAGCCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTAcga                          >  1:1052620/1‑67 (MQ=9)
                                              cacaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAttt                               >  1:1567424/1‑63 (MQ=255)
                                              cacaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTc                              >  1:1156252/1‑64 (MQ=255)
                                              cacaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAAc                       <  1:1513145/71‑1 (MQ=255)
                                                caTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                            >  1:631642/1‑48 (MQ=255)
                                                caTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                            >  1:1659681/1‑48 (MQ=255)
                                                caTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                            >  1:5185732/1‑48 (MQ=255)
                                                   cTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCa                     >  1:5511648/1‑68 (MQ=255)
                                                    tttGCTGCCATCTTTCGGCGTGATGAAACCAAAAcc                                                    >  1:1263479/1‑36 (MQ=255)
                                                    tttGCTGCCATCTTTCGGCGTGATGAAACCAAAAcc                                                    >  1:1587214/1‑36 (MQ=255)
                                                    tttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                    <  1:5616151/52‑1 (MQ=255)
                                                    tttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCa                 >  1:3707685/1‑71 (MQ=255)
                                                     ttGCTGCCATCTTTCGGCGTGATGAAACCAAAAcc                                                    >  1:4633623/1‑35 (MQ=255)
                                                     ttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGt                   >  1:3644293/1‑68 (MQ=255)
                                                     ttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAt                >  1:5221626/1‑71 (MQ=255)
                                                      tgctgcCATCATCAGGAGTGATAAAGCCAAAACCTTTATCTGCGTTAAACCATTTTACTAAAc                       >  1:4022803/1‑63 (MQ=9)
                                                       gctgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAtt               <  1:4479250/70‑1 (MQ=255)
                                                       gctgcCATCATCAGGAGTGATAAAGCCAAAACCTTTATCTGCGTTAAACCATTTTACTAAACCAGTCAtt               <  1:2100745/70‑1 (MQ=12)
                                                        gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAttt              <  1:82126/68‑1 (MQ=255)
                                                        gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAttt              <  1:3892167/68‑1 (MQ=255)
                                                        gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAttt              <  1:2003423/68‑1 (MQ=255)
                                                        ctgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAttt              <  1:1591862/70‑1 (MQ=255)
                                                          gcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                    <  1:3448645/46‑1 (MQ=255)
                                                           ccATCTTTGGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt           >  1:2737445/1‑70 (MQ=255)
                                                           ccATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt           <  1:4351269/70‑1 (MQ=255)
                                                             aTCTTTCGGTGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                          <  1:3328034/53‑1 (MQ=255)
                                                             aTCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                          <  1:3092487/53‑1 (MQ=255)
                                                             aTCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                          <  1:160521/53‑1 (MQ=255)
                                                                tttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTa         >  1:5261761/1‑67 (MQ=255)
                                                                 ttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACa     >  1:4028751/1‑70 (MQ=255)
                                                                 ttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACa     >  1:2488971/1‑70 (MQ=255)
                                                                  tCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACAt    >  1:2333624/1‑70 (MQ=255)
                                                                   cGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATaa  >  1:151545/1‑71 (MQ=255)
                                                                   |                                                                      
TTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAA  >  minE/980169‑980306

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: ‑