Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 3,725,043 C→A A42A (GCG→GCT)  U712_19205 ← Uncharacterized protein ywcI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,725,0430CA96.8% 100.5 / NA 32A42A (GCG→GCT) U712_19205Uncharacterized protein ywcI
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (5/25);  total (6/26)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAA‑CGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGT‑CCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAA  >  CP006881/3724909‑3725172
                                                                                                                                       |                                                                                                                                  
tcCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAA‑CGAAACGAGTCCTGCGGCCTGCAcagc                                                                                                                                  <  1:129988/137‑1 (MQ=255)
  cTGCTTTGGCGATTCCTTGAAGGATTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACagccgcc                                                                                                                              <  2:299038/139‑2 (MQ=255)
     cTTTGGCGATTCCTTGAAGGCTTGCAAACATCAACTTTAATAAAGCGCGGCGGCTTTCGCGTTAATGTTTTATTACTCAATCACACCCCTCCCCCAGTCATAAAA‑CGAAACGAGTCCTGCGGCCCGCAAcaccgcat                                                                                                                             >  8:403301/1‑135 (MQ=255)
          gCGATTCCTTGAAGGGTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCCTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                                                                                       <  7:416567/138‑1 (MQ=255)
                    gAAGGCTTGAAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGt                                                                                                                <  3:256480/135‑1 (MQ=255)
                     aaGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                                                                                           <  7:21900/139‑1 (MQ=255)
                     aaGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCCTCCTCATTCACAGCCATCCGCCATTCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                                                                                           <  8:43493/139‑1 (MQ=255)
                       cgCTTGAAACCTTCTTCTTTATTAACGGGGGGCGGCTTTGGATTTGTGCTTTGCTTCCTCATTCCCCCCCACCCCCATTTCAT‑AAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                                                        <  1:101369/138‑1 (MQ=255)
                        gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                                                                                           <  5:61767/136‑1 (MQ=255)
                        gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGa                                                                                                         <  2:69459/138‑1 (MQ=255)
                             cAAACAGCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGcc                                                                                                      <  5:359578/136‑1 (MQ=255)
                                             taACGCGCGGGGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCTTCCGCCATTCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACTGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAg                                                                                      <  5:435830/136‑1 (MQ=255)
                                                          cTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGT‑c                                                                     <  2:48684/139‑1 (MQ=255)
                                                          cTTTCGCTTTCTGCTTTTCTTCCTCAATCACAGCCATCCGCCATTCATAAAC‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCt                                                                         <  5:257692/136‑1 (MQ=255)
                                                                      gcttttCTTCCTCATTCACAGCCATCCGCCATTCATAAAA‑CGAAACGAGTCCTGCGGCCTGCCCCGCCGCCGTTGGTTTGGATGTTCGCGGAGGCGGATAGACCCCGGCCGCCCCCGGTTTGTGGT‑CCGGCggccc                                                            >  5:435630/3‑132 (MQ=255)
                                                                             ttcCTCATTCACAGCCATCCGCCATTCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCGTTATTTGGTTTTTTCGCTGAGGCGCATATAACCCGGCAGCCCCAGTTTTTCTGT‑GCAGCGGCACCCAGtt                                                      >  6:248871/1‑135 (MQ=255)
                                                                                   aTTCACAGCCATCCGCCATTCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGT‑CCGGCGTCACCCAGTTCATCaaaaa                                             <  1:114203/138‑1 (MQ=255)
                                                                                   aTTCACAGCCATCCGCCATTCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGT‑CCGGCGTCACCCAGTTCATCaaaa                                              <  1:108639/137‑1 (MQ=255)
                                                                                    ttCACAGCCATCCGCCATTCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGT‑CCGGCGTCACCCAGTTCATCaaaa                                              <  6:425690/136‑1 (MQ=255)
                                                                                         aGCCATCCGCCATTCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGT‑CCGGCGTCACCCAGTTCATCAAAAAGAcc                                         <  5:270749/136‑1 (MQ=255)
                                                                                         aGCCATCCGCCATTCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGT‑CCGGCGTCACCCAGTTCATCAAAAAGAcc                                         <  3:104669/136‑1 (MQ=255)
                                                                                           cactCCGCCATTCATAAAA‑CGAACCGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGT‑CCGGCGTCACCCAGTTCATCAAAAAGACCATCca                                    <  8:373078/136‑1 (MQ=255)
                                                                                            cctCCCCCATTAATAAAACCAAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGT‑CCGGCGTCACCCAGTTCATCAAAAAGACCATCca                                    <  1:352659/137‑1 (MQ=255)
                                                                                              tCCGCCATTCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGT‑CCGGCGTCACCCAGTTCATCAAAAAGa                                           >  1:61959/1‑129 (MQ=255)
                                                                                              tCCGCCATTCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGT‑CCGGCGTCACCCAGTTCATCAAAAAGa                                           <  2:61959/129‑1 (MQ=255)
                                                                                               ccGCCATTCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCGTTGTTTGGTTTGTACGCTGAGCCGCATAGACCGCGGCAGCCGCAGTTTTTCTGTGCCGGGGGC‑CCCAGGTGATCCAAAACACCAACCCCAcc                                >  7:149195/1‑139 (MQ=255)
                                                                                                   cattCATAAAA‑CGAAACGAGTCCTGCGGCCTGGACAGCCGGCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGT‑CCGGCGTCTCCCAGTTCATCAAAAAGACCATCCTCACCCGTa                            <  8:146262/136‑1 (MQ=255)
                                                                                                   cattCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGT‑CCGGCGTCACCCAGTTCATCAAAAAGACCATCCACAccc                               <  4:345412/133‑1 (MQ=255)
                                                                                                   cattCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGT‑CCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTa                            <  1:232107/136‑1 (MQ=255)
                                                                                                    attCATAAAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGT‑CCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCg                              <  4:113986/134‑1 (MQ=255)
                                                                                                         tAAAA‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGT‑CCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAg                         <  4:164361/136‑1 (MQ=255)
                                                                                                          aaaa‑CGAAACGAGTCCTGCGGCCTGCACAGCCGCCGTTCTTTGGTTTGTACGCGGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGT‑CCGGCCGCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAg                         >  3:185941/1‑135 (MQ=255)
                                                                                                                        tCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGCGCCGAATAGACCGCGGCAGCCGCAGGTTTTCTGT‑GCGGCGTCACCCAGTTCATAAAAAAGACCCTCCCAACCCGGAAATATACAAGCAAACCttt         >  7:85545/1‑138 (MQ=255)
                                                                                                                            gcgGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGT‑CCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCATACGTTTCATAg    <  8:154652/139‑1 (MQ=255)
                                                                                                                                 cTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGT‑CCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGaa  <  3:164087/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                  
TCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAA‑CGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGT‑CCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAA  >  CP006881/3724909‑3725172

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑