Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,892,882 (C)6→7 coding (354/1911 nt) yoaA ← putative 5' to 3' DNA helicase implicated in DNA repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,892,8761.C97.1% 109.0 / ‑3.4 35coding (360/1911 nt)yoaAputative 5' to 3' DNA helicase implicated in DNA repair
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (16/18);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCG‑‑CCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAG  >  NC_000913/1892753‑1893001
                                                                                                                            |                                                                                                                              
gACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGcc                                                                                                                <  1:841901/140‑1 (MQ=255)
        gCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTc                                                                                                        >  1:705505/1‑140 (MQ=255)
        gCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTc                                                                                                        >  1:629296/1‑140 (MQ=255)
                tgtgAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGtt                                                                                                <  2:617411/140‑1 (MQ=255)
                 gtgAATCTTCCGCCACGCTGACGCAGGTGCTGATATCCCCATCGACTGTTTGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTc                                                                                               <  1:144416/140‑1 (MQ=255)
                     aTCTTCCGCCACGCTGACGCAGGTGCTGATTTCACCATCGACTGTTTGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAgg                                                                                           <  2:629296/140‑1 (MQ=255)
                     aTCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAgg                                                                                           <  2:705505/140‑1 (MQ=255)
                       cTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCa                                                                                         <  2:663103/140‑1 (MQ=255)
                           cGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAgg                                                                                     <  1:6367/140‑1 (MQ=255)
                                     aCGCAGGTGCTGATATCACCATCGACTGTTTGATTAGGCCAGGAGCGCAGCAGGATCAAATCGCTTAAGATTTGTCCCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAgc                                                                           <  1:87653/140‑1 (MQ=255)
                                           gTGCTGATATCACCATCGACTGTTTGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCtt                                                                     >  2:827266/1‑140 (MQ=255)
                                           gTGCTGATATCACCATCGACTGTTTGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCtt                                                                     >  1:142777/1‑140 (MQ=255)
                                                      aCCATCGACTGTTTGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGCGcc                                                          >  2:573860/1‑140 (MQ=255)
                                                       ccATCGACTGTTTGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTTCCGGCAGATCG‑‑CCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGCGCCAc                                                        <  2:519212/140‑1 (MQ=255)
                                                             aCTGTTTGATTAGACAAGGAGCCCATCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGCGCCACGTTGc                                                   <  1:351177/140‑1 (MQ=255)
                                                              cTGTTTGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGCGCCACGTTGcc                                                  >  2:665997/1‑140 (MQ=255)
                                                               tGTTTGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCATCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGCGCCACGTTGccc                                                 >  2:58349/1‑140 (MQ=255)
                                                                  ttGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGACCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGCGCCACGTTGCCCGta                                              <  1:15704/140‑1 (MQ=255)
                                                                     aTTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGCGCCACGTTGCCCGTatat                                           <  1:78920/140‑1 (MQ=255)
                                                                         gACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGCGCCACGTTGCCCGTATATTTCa                                       <  1:665997/140‑1 (MQ=255)
                                                                           ccAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCATCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGCGCCACGTTGCCCGTATATTTCAAt                                     <  1:58349/140‑1 (MQ=255)
                                                                           ccAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGCGCCACGTTGCCCGTATATTTCAAt                                     <  1:573860/140‑1 (MQ=255)
                                                                                    cagcagGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGTAGCGCCACGTTGCCCGTATATTTCAATGCCTTTgag                            >  2:560388/1‑140 (MQ=255)
                                                                                    cagcagGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTgag                            >  2:860995/1‑140 (MQ=255)
                                                                                       cagGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGGCCTTTCAGAAGCGCCACGTTGCCCGGATATTTCAATGCCTTTGAGACt                         >  2:301388/1‑140 (MQ=255)
                                                                                          gATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGGGCCTCGTGGCCCGTATGTTTCTATGCCTTTGAGACTGtt                      >  2:283879/1‑140 (MQ=255)
                                                                                            tCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTgg                    >  1:603564/1‑140 (MQ=255)
                                                                                            tCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTCCAGCAGCGTCACGTTGACCGTATATTTCAAGGCCTTTGAGACTGTTgg                    >  2:262537/1‑140 (MQ=255)
                                                                                                  cGCTTAAGATTTGTACCGGCAGGTCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTgaga                           <  2:166566/127‑1 (MQ=255)
                                                                                                  cGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTgaga                           >  1:166566/1‑127 (MQ=255)
                                                                                                   gCTTAAGATTTGTACCGGCAGAACGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATcg             <  1:283879/140‑1 (MQ=255)
                                                                                                    cTTAAGATTTGTACCGGCAGATCGCCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATcg             >  2:25595/1‑140 (MQ=255)
                                                                                                      tAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCgg          <  1:301388/140‑1 (MQ=255)
                                                                                                         gATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTg       <  1:827266/140‑1 (MQ=255)
                                                                                                              gTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGCGCCACGTTGACCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTagag  >  2:730871/1‑140 (MQ=255)
                                                                                                                            |                                                                                                                              
GACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCG‑‑CCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAG  >  NC_000913/1892753‑1893001

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: