Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,126,130 A→T 10.2% G227G (GGT→GGA)  priA ← Primosome factor n' (replication factor Y)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,126,1300AT10.2% 155.8 / 7.8 49G227G (GGT→GGA) priAPrimosome factor n' (replication factor Y)
Reads supporting (aligned to +/- strand):  ref base A (23/21);  new base T (2/3);  total (25/24)
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.82e-01

ACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCAC  >  NC_000913/4125993‑4126278
                                                                                                                                         |                                                                                                                                                    
aCGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGTTCCGGTaaaacg                                                                                                                                           <  3:449837/149‑5 (MQ=255)
aCGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAAcgcc                                                                                                                                           >  1:355870/1‑149 (MQ=255)
 cGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAAcgccc                                                                                                                                          <  3:469446/149‑1 (MQ=255)
  gggggCATTAAAACGTTCACGAAAACGGGCGATAGATTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGTTCCGGTaaaacgga                                                                                                                                         <  3:324934/149‑7 (MQ=255)
    gggCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAAcgcccgcc                                                                                                                                       >  4:264656/1‑149 (MQ=255)
            aaaCGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCa                                                                                                                               >  2:160639/1‑149 (MQ=255)
                  tCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGa                                                                                                                         >  3:287012/1‑149 (MQ=255)
                   cACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGaa                                                                                                                        >  4:466072/1‑149 (MQ=255)
                      gAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGaaaaa                                                                                                                     <  1:14090/149‑1 (MQ=255)
                      gAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGaaaaa                                                                                                                     <  3:361641/149‑1 (MQ=255)
                         aaCGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTa                                                                                                                  <  1:124586/149‑1 (MQ=255)
                           cGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATc                                                                                                                <  4:69557/149‑1 (MQ=255)
                              gCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGc                                                                                                             <  4:386203/149‑1 (MQ=255)
                                    gTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTAt                                                                                                                 <  4:136037/139‑1 (MQ=255)
                                    gTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTAt                                                                                                                 >  3:136037/1‑139 (MQ=255)
                                    gTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCt                                                                                                       >  1:306228/1‑149 (MQ=255)
                                     tttGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTa                                                                                                      >  1:183114/1‑149 (MQ=255)
                                      ttGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTAt                                                                                                     >  2:265078/1‑149 (MQ=255)
                                             gtCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTgc                                                                                              <  2:183114/149‑1 (MQ=255)
                                              tCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTgcg                                                                                             <  3:466072/149‑1 (MQ=255)
                                                  gCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTAt                                                                                                                 <  1:48311/125‑1 (MQ=255)
                                                  gCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTAt                                                                                                                 >  2:48311/1‑125 (MQ=255)
                                                       tttCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAAcggcg                                                                                    >  2:181429/1‑149 (MQ=255)
                                                         tCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGt                                                                                  <  2:385308/149‑1 (MQ=255)
                                                                   tCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCg                                                                        <  1:265078/149‑1 (MQ=255)
                                                                    cACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCgcgc                                                                                                        >  1:133918/1‑116 (MQ=255)
                                                                    cACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCgcgc                                                                                                        <  2:133918/116‑1 (MQ=255)
                                                                      ccAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTa                                                                     >  2:280326/1‑149 (MQ=255)
                                                                       cAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTAt                                                                    >  2:338461/1‑149 (MQ=255)
                                                                             cTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATc                                                              <  4:287012/149‑1 (MQ=255)
                                                                                     ccTGAGCGACCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGc                                                      >  1:408933/1‑149 (MQ=255)
                                                                                          gcgagcACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAg                                                                                                                              <  4:58756/72‑1 (MQ=255)
                                                                                          gcgagcACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAg                                                                                                                              >  3:58756/1‑72 (MQ=255)
                                                                                                  gTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACgg                                         >  2:376693/1‑149 (MQ=255)
                                                                                                  gTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACgg                                         >  2:178949/1‑149 (MQ=255)
                                                                                                        cAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAgt                                   <  1:160639/149‑1 (MQ=255)
                                                                                                          gTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAgttc                                 >  3:335393/1‑149 (MQ=255)
                                                                                                            aCGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAgttcgt                               <  1:338461/149‑1 (MQ=255)
                                                                                                             cGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAgttcgtt                              >  2:14914/1‑149 (MQ=255)
                                                                                                                    aTAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGt                       >  2:157507/1‑149 (MQ=255)
                                                                                                                       aaCCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGc                    <  1:178949/149‑1 (MQ=255)
                                                                                                                       aaCCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGc                    <  2:355870/149‑1 (MQ=255)
                                                                                                                          cTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTaa                 >  4:250133/1‑149 (MQ=255)
                                                                                                                          cTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTaa                 >  2:31146/1‑149 (MQ=255)
                                                                                                                                gttACCGGTTCCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTa                                                                     >  3:287763/2‑91 (MQ=255)
                                                                                                                                gttACCGGTTCCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTa                                                                     <  4:287763/90‑1 (MQ=255)
                                                                                                                                gttACCGGTTCCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTg           >  3:328998/2‑149 (MQ=255)
                                                                                                                                  tACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGgt         <  1:376693/149‑1 (MQ=255)
                                                                                                                                      ggAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGttt     <  3:264656/149‑1 (MQ=255)
                                                                                                                                         aCCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCAc  <  4:335393/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
ACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCAC  >  NC_000913/4125993‑4126278

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