Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,091,167 T→C 8.4% E178E (GAA→GAG sufB ← component of SufBCD complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,091,1670TC 8.4% 199.5 / 8.9 71E178E (GAA→GAGsufBcomponent of SufBCD complex
Reads supporting (aligned to +/- strand):  ref base T (24/41);  new base C (3/3);  total (27/44)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

CCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGAT  >  Exported/1091100‑1091237
                                                                   |                                                                      
ccgcATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTcc                                                                     <  1:597576/71‑1 (MQ=255)
ccgcATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTcc                                                                     <  1:3368319/71‑1 (MQ=255)
   cATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTcc                                                                     <  1:265448/68‑1 (MQ=255)
    aTTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTcc                                                                     <  1:222361/67‑1 (MQ=255)
     ttAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTg                                                                <  1:3340312/71‑1 (MQ=255)
        aGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGAt                                                              <  1:799874/70‑1 (MQ=255)
        aGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGAt                                                              <  1:2267286/70‑1 (MQ=255)
         gcgcGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCg                                                            <  1:3425612/71‑1 (MQ=255)
          cgcgGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGt                                                           <  1:4218521/71‑1 (MQ=255)
                aaaGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCg                                                     >  1:2084387/1‑71 (MQ=255)
                 aagaagTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTcc                                                                     <  1:1764664/54‑1 (MQ=255)
                  agaagTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGa                                                               >  1:2293328/1‑59 (MQ=255)
                  agaagTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGAt                                                              >  1:1122207/1‑60 (MQ=255)
                     agTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCa                                                <  1:3105040/71‑1 (MQ=255)
                     agTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCa                                                <  1:540256/71‑1 (MQ=255)
                       ttGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATc                                                      <  1:1712841/63‑1 (MQ=255)
                       ttGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATc                                                      <  1:3327180/63‑1 (MQ=255)
                         gTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCg                                                            >  1:242711/1‑55 (MQ=255)
                         gTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCg                                                            >  1:2432777/1‑55 (MQ=255)
                            aTTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAgg                                         <  1:680632/71‑1 (MQ=255)
                            aTTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAg                                          <  1:2949760/70‑1 (MQ=255)
                            aTTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAg                                          >  1:2838370/1‑70 (MQ=255)
                              tCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTgg                                                         <  1:1832820/53‑1 (MQ=255)
                              tCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTgg                                                         <  1:2858472/53‑1 (MQ=255)
                                   ggCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAg                                    >  1:4270373/1‑69 (MQ=255)
                                   ggCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGa                                   >  1:3832116/1‑70 (MQ=255)
                                   ggCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGa                                   >  1:4746705/1‑70 (MQ=255)
                                   ggCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGa                                   >  1:309768/1‑70 (MQ=255)
                                   ggCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGa                                   >  1:1133739/1‑70 (MQ=255)
                                     caccacGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGaaaa                                <  1:3899393/71‑1 (MQ=255)
                                      accacGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGaaaa                                >  1:1800203/1‑70 (MQ=255)
                                      accacGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaat                               <  1:507586/71‑1 (MQ=255)
                                      accacGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaat                               <  1:3562289/71‑1 (MQ=255)
                                      accacGGTGACGATATATTTACGCACCAGTTCCGGGTTAt                                                              <  1:2058878/40‑1 (MQ=25)
                                        cacGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaataa                             >  1:2812267/1‑71 (MQ=255)
                                        cacGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaataa                             <  1:4216535/71‑1 (MQ=255)
                                         acGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGaaa                                 >  1:2258823/1‑66 (MQ=255)
                                         acGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaataat                            <  1:683617/71‑1 (MQ=255)
                                          cGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACa                                     >  1:2100885/1‑61 (MQ=255)
                                          cGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaataat                            >  1:133924/1‑70 (MQ=255)
                                                cGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGc                     <  1:2256171/71‑1 (MQ=255)
                                                 gagaTATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCt                    <  1:3535621/71‑1 (MQ=255)
                                                 gagaTATTTACGCACCAGTTCCGGGTGATCGGGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCt                    <  1:2131610/71‑1 (MQ=255)
                                                  agaTATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaat                               <  1:4214459/59‑1 (MQ=255)
                                                  agaTATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaat                               <  1:1649766/59‑1 (MQ=255)
                                                  agaTATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaat                               <  1:3755994/59‑1 (MQ=255)
                                                  agaTATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaat                               <  1:473039/59‑1 (MQ=255)
                                                    atatTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTcc                  >  1:474965/1‑70 (MQ=255)
                                                      atTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGc                >  1:4752184/1‑70 (MQ=255)
                                                      atTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGc                >  1:677471/1‑70 (MQ=255)
                                                         tACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTg                      >  1:3205025/1‑61 (MQ=255)
                                                         tACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTg                      >  1:2342086/1‑61 (MQ=255)
                                                         tACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAg             <  1:3644013/70‑1 (MQ=255)
                                                         tACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAg             <  1:3675767/70‑1 (MQ=255)
                                                         tACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGt            >  1:4145124/1‑71 (MQ=255)
                                                         tACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGt            >  1:272152/1‑71 (MQ=255)
                                                         tACGCACCAGTT‑‑GGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTg                      >  1:2374930/1‑59 (MQ=255)
                                                            gCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAGCAGAAAATAATTCCATGc                     <  1:2416109/59‑1 (MQ=255)
                                                            gCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGc                     <  1:1971759/59‑1 (MQ=255)
                                                             cACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGaaaa                                <  1:4636224/47‑1 (MQ=255)
                                                             cACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGaaaa                                <  1:4020177/47‑1 (MQ=255)
                                                             cACCAGTTCCGGGTGATCGTGAATCGCCTCACCAAAGGAACAGaaaa                                <  1:3953104/47‑1 (MQ=255)
                                                              aCCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAAtt                           <  1:2715858/51‑1 (MQ=255)
                                                              aCCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAAtt                           <  1:390298/51‑1 (MQ=255)
                                                              aCCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAAtt                           <  1:2274847/51‑1 (MQ=255)
                                                                cAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTcgc     <  1:1456337/71‑1 (MQ=255)
                                                                cAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTcgc     >  1:2807026/1‑71 (MQ=255)
                                                                cAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTcgc     <  1:881238/71‑1 (MQ=255)
                                                                 aGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTcgcg    <  1:2901376/71‑1 (MQ=255)
                                                                  gTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCGAGTTTTTCGCGa   >  1:3634092/1‑71 (MQ=255)
                                                                  gTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGa   >  1:2313721/1‑71 (MQ=255)
                                                                   ttCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGaaaa                                >  1:4030721/1‑41 (MQ=255)
                                                                   ttCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGaaaa                                >  1:3070067/1‑41 (MQ=255)
                                                                   ttCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGaaaa                                >  1:55424/1‑41 (MQ=255)
                                                                   ttCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGaaaa                                >  1:1318966/1‑41 (MQ=255)
                                                                   ttCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaataa                             <  1:3482526/44‑1 (MQ=255)
                                                                   ttCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGAt  <  1:1354566/71‑1 (MQ=255)
                                                                   ttCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGAt  <  1:1276540/71‑1 (MQ=255)
                                                                   |                                                                      
CCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGAT  >  Exported/1091100‑1091237

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: