Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,328,359 T→C 14.7% N71S (AAC→AGC)  greA ← transcription elongation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,328,3590TC14.7% 75.0 / 8.1 34N71S (AAC→AGC) greAtranscription elongation factor
Reads supporting (aligned to +/- strand):  ref base T (12/17);  new base C (3/2);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

GTTACGGTAGCACCAAAAATAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGAAACCCTGCTGTTCAC  >  W3110S.gb/3328293‑3328416
                                                                  |                                                         
gTTACGGTAGCACCAAAAATAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGCTCGa                                                       >  1:2805308/1‑71 (MQ=255)
      gTAGCACCAAAAATAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTc                                                         <  1:898623/63‑1 (MQ=255)
      gTAGCACCAAAAATAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTc                                                         <  1:870698/63‑1 (MQ=255)
         gCACCAAAAATAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGCTCGACAGCTTGGc                                              >  1:1395862/1‑71 (MQ=255)
             cAAAAATAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCg                                          <  1:1003122/71‑1 (MQ=255)
                  aTAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGCTCGACAGCTTGGCTTCGAtctc                                     <  1:1202940/71‑1 (MQ=255)
                      agcgcCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCttt                                  <  1:2346941/69‑1 (MQ=255)
                       gcgcCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTaa                                <  1:2430975/71‑1 (MQ=255)
                       gcgcCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTaa                                <  1:1803855/71‑1 (MQ=255)
                         gcCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAAt                               <  1:673748/70‑1 (MQ=255)
                             aTTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTc                                                         <  1:1999038/40‑1 (MQ=255)
                             aTTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTc                                                         <  1:238919/40‑1 (MQ=255)
                              ttgttgGGCATTTTGGTGACATCTATCACCTGCGCGCTCGACAGCTTGGCTTCGAt                                        <  1:2639888/56‑1 (MQ=39)
                              ttgttgGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTc                                           <  1:895341/53‑1 (MQ=255)
                              ttgttgGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGAt                                        <  1:443812/56‑1 (MQ=255)
                              ttgttgGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGcc                         <  1:1550375/71‑1 (MQ=255)
                                gttgGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCtt                       >  1:1497985/1‑71 (MQ=255)
                                gttgGGCATTTTGGTGACATCAATCACCTGCGCGCTCGACAGCTTGGCTTCGATGt                                      >  1:1029709/1‑56 (MQ=255)
                                   gggCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCg                     <  1:2489809/70‑1 (MQ=255)
                                     gCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTc                                           <  1:2629798/46‑1 (MQ=255)
                                        ttttGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTaa                                <  1:971808/54‑1 (MQ=255)
                                          ttGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGt                                      >  1:25210/1‑46 (MQ=255)
                                          ttGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGAAAcc             >  1:667747/1‑71 (MQ=255)
                                                  atcaatcaCCTGTGCGTTCGACAGCTTGGCTTCGATGTCt                                    >  1:100163/1‑40 (MQ=255)
                                                  atcaatcaCCTGCGCGTTCGACAGCTTGGCTTCGATGTCt                                    >  1:2571228/1‑40 (MQ=255)
                                                  atcaatcaCCTGCGCGTTCGACAGCTTGGCTTCGATGTCt                                    >  1:2018144/1‑40 (MQ=255)
                                                  atcaatcaCCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTaa                                >  1:1282690/1‑44 (MQ=255)
                                                   tcaatcaCCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACg                            <  1:1298847/47‑1 (MQ=255)
                                                     aatcaCCTGCGCGTTCGACAGCTTGGCTTCGATGTCttt                                  <  1:2419374/39‑1 (MQ=255)
                                                     aatcaCCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGAAACCCTGCTGTTCAc  >  1:398022/1‑71 (MQ=255)
                                                     aatcaCCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGAAACCCTGCTGTTCAc  >  1:2264745/1‑71 (MQ=255)
                                                         aCCTGCGCTTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGaaa               >  1:1527082/1‑54 (MQ=255)
                                                         aCCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGc                          >  1:597709/1‑43 (MQ=255)
                                                         aCCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCg                     >  1:441980/1‑48 (MQ=255)
                                                                  ttCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCa                   <  1:2641065/41‑1 (MQ=255)
                                                                  |                                                         
GTTACGGTAGCACCAAAAATAACGCGCCCATTGTTGGGCATTTTGGTGACATCAATCACCTGCGCGTTCGACAGCTTGGCTTCGATGTCTTTAATACGGCCTTCGCAGAAACCCTGCTGTTCAC  >  W3110S.gb/3328293‑3328416

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: