Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,126,130 A→T 6.3% 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,1300AT 6.3% 275.1 / 6.7 79G227G (GGT→GGA) priAPrimosome factor n' (replication factor Y)
Reads supporting (aligned to +/- strand):  ref base A (41/33);  new base T (2/3);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

CCGAGATTTTTAAACGGCGTAAACAGCGCGGAGCGGGTGCCGATCACAATCGCCGCCTCACCATTTTTCGCTTTCAGCCACGCCGAAAGACGCTCGCTGTCGTTCAGGCCGGAATGCAGAACTTCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGATTTTGCCTTGCCGTAACGCCGCCAGCGCCTGTTG  >  NC_000913/4125867‑4126396
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                          
ccGAGATTTTTAAACGGCGTAAACAGCGCGGAGCGGATGCCGCTCACAATCGCGGCCTCACCATTTTTCGCTTTCAGCCACGCCGAAAGACGCTCGCTGTCGTTCAGGCCGGAATGCAGAACTTCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAAcgcc                                                                                                                                                                                                                                                                 <  2:144838/275‑1 (MQ=255)
                         gcgcGGAGCGGGTGCCGATCACAATCGCCGCCTCACCATTTTTCGCTTTCAGCCACGCCGAAAGACGCTCGCTGTCGTTCAGGCCGGAATGCAGAACTTCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCCGGCAGAAAAAGTa                                                                                                                                                                                                                                        >  2:237004/1‑275 (MQ=255)
                          cgcgGAGCGGGTGCCGATCACAATCGCCGCCTCACCATTTTTCGCTTTCAGCCAAGCCGAAAGACGCTCGCTGTCGTTCAGGCCGGAATGCAGAACTTCCACGGGGGCATTAAAACGTTAACGAAAACGGGCGATAGTTTGCGGTTGCAGGCCGATTTCCGGCACAATCCCCAGCGACTGTTTGCCCTGAGCGAGCATGGTTTCCAGTACGACGAGATGAACCCCAGTTTTACCGGTACCGGTCAACCCCGTCGGCCGCCCGGCAAGAAAATTAt                                                                                                                                                                                                                                       >  4:384960/1‑275 (MQ=255)
                                gCGGGTGCCGATCACAATCGCCGCCTCACCATTTTTCGCTTTCAGCCACGCCGAAAGACGCTCGCTGTCGTTCAGGCCGGAATGCTGAACATCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGTTAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGcc                                                                                                                                                                                                                                  <  2:360106/274‑1 (MQ=255)
                                gCGGGTGCCGATCACAATCGCCGCCTCACCAGTTTTCGCTTTCAGCCACGCCGAAAGACGCTCGCTGTCGTTCAGGCCGGAATGCAGAACTTCCACGTGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGcc                                                                                                                                                                                                                                  <  2:231371/274‑1 (MQ=255)
                                                   gccgccTCACCATTTTTCGCTTTCAGCCACGCCGAAAGACGCTCGCTGTCGTTCAGGCCGGAATGCAGAACTTCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCaa                                                                                                                                                                                                               <  2:242381/274‑1 (MQ=255)
                                                            ccATTTTGCGCTTTCAGCCAAGCCGAAAGACGCTCGCTGTCGTCCAGGCCGGAATGCAGAACTTCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGc                                                                                                                                                                                                     <  1:54216/275‑1 (MQ=255)
                                                                       tttCAGCCACGCCGAAAGACGCTCGCTGTCGTTCAGGCCGGAATGCAGAACTTCCACGGGGGCATTAAACCGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACCGCGGTGGCCTGTTCGGTa                                                                                                                                                                                           >  2:106758/1‑274 (MQ=255)
                                                                                                           gCCGGAATGCAGAACTTCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTACCGCCCGCCAGCAGCCAGGCAAAAAAAGTATCCGCCGCCCTATGCATTGACGCAACCGCGGTGGACTGTTCGGTATTTCACCTCAACCCCTCCCCACAACCggaccagttc                                                                                                                                                       >  4:358819/1‑266 (MQ=255)
                                                                                                               gAATGCAGAACTTCCACGGGGGCATTAAACCGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATa                                                                                                                                                           <  2:358824/266‑1 (MQ=255)
                                                                                                                          ttCCACGGGTGCAGTAAAACGTTCACGAAAACGGGCGATAGTTGGGGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGGTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCGCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACAGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCAGAGTTCGTTCGCCAGTCGCTa                                                                                                                                        <  2:94377/274‑1 (MQ=255)
                                                                                                                          ttCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTTTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCt                                                                                                                                         <  2:384965/273‑1 (MQ=255)
                                                                                                                              aCGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGTTCCGGTaaaacg                                                                                                                                                                                                                                                                 <  7:449837/149‑5 (MQ=255)
                                                                                                                              aCGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAAcgcc                                                                                                                                                                                                                                                                 >  5:355870/1‑149 (MQ=255)
                                                                                                                               cGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAAcgccc                                                                                                                                                                                                                                                                <  7:469446/149‑1 (MQ=255)
                                                                                                                                gggggCATTAAAACGTTCACGAAAACGGGCGATAGATTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGTTCCGGTaaaacgga                                                                                                                                                                                                                                                               <  7:324934/149‑7 (MQ=255)
                                                                                                                                  gggCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAAcgcccgcc                                                                                                                                                                                                                                                             >  8:264656/1‑149 (MQ=255)
                                                                                                                                      aTTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGtt                                                                                                                            >  2:10774/1‑274 (MQ=255)
                                                                                                                                          aaaCGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCa                                                                                                                                                                                                                                                     >  6:160639/1‑149 (MQ=255)
                                                                                                                                                tCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGa                                                                                                                                                                                                                                               >  7:287012/1‑149 (MQ=255)
                                                                                                                                                 cACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGaa                                                                                                                                                                                                                                              >  8:466072/1‑149 (MQ=255)
                                                                                                                                                 cACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTaaa                                                                                                                >  2:285651/1‑275 (MQ=255)
                                                                                                                                                 cACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACCNCCGCCAGCAGCCAGGCCGAAAAAGTATCTGCAGCGCTATGAATTgcg                                                                                                                                                                                                                   >  3:52485/1‑176 (MQ=255)
                                                                                                                                                  aCGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCAAAGTTCTTTCCCCATTCCCTAAACTCTGGTTTTTCACTTGTCAAAt                                                                                                               >  4:402123/1‑275 (MQ=255)
                                                                                                                                                    gAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGaaaaa                                                                                                                                                                                                                                           <  5:14090/149‑1 (MQ=255)
                                                                                                                                                    gAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGaaaaa                                                                                                                                                                                                                                           <  7:361641/149‑1 (MQ=255)
                                                                                                                                                     aaaaCGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCaa                                                                                                                                                                                                               <  1:52485/176‑1 (MQ=255)
                                                                                                                                                       aaCGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTa                                                                                                                                                                                                                                        <  5:124586/149‑1 (MQ=255)
                                                                                                                                                         cGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATc                                                                                                                                                                                                                                      <  8:69557/149‑1 (MQ=255)
                                                                                                                                                            gCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGc                                                                                                                                                                                                                                   <  8:386203/149‑1 (MQ=255)
                                                                                                                                                                  gTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTAt                                                                                                                                                                                                                                       <  8:136037/139‑1 (MQ=255)
                                                                                                                                                                  gTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTAt                                                                                                                                                                                                                                       >  7:136037/1‑139 (MQ=255)
                                                                                                                                                                  gTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCt                                                                                                                                                                                                                             >  5:306228/1‑149 (MQ=255)
                                                                                                                                                                   tttGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTa                                                                                                                                                                                                                            >  5:183114/1‑149 (MQ=255)
                                                                                                                                                                    ttGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTAt                                                                                                                                                                                                                           >  6:265078/1‑149 (MQ=255)
                                                                                                                                                                     tGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGc                                                                                             >  2:315303/1‑274 (MQ=255)
                                                                                                                                                                           gtCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTgc                                                                                                                                                                                                                    <  6:183114/149‑1 (MQ=255)
                                                                                                                                                                            tCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTgcg                                                                                                                                                                                                                   <  7:466072/149‑1 (MQ=255)
                                                                                                                                                                                gCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTAt                                                                                                                                                                                                                                       >  6:48311/1‑125 (MQ=255)
                                                                                                                                                                                gCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTAt                                                                                                                                                                                                                                       <  5:48311/125‑1 (MQ=255)
                                                                                                                                                                                    aTTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCCCGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTCCGCCAACGGCGGCGGCCTGTTCGGTATTCCATCGCAACCGCTCACCCGAACCGGCATAGTTCGTTCGCCCGTCCCTAAACTCTGGTGTTTCCCCTGCTCAGTCACCCAGCGCTTTTTTGCTAACCGCATTCAAcg                                                                              >  4:321596/1‑274 (MQ=255)
                                                                                                                                                                                    aTTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAAcggcgg                                                                                                                                                                                                         >  2:521229/1‑151 (MQ=255)
                                                                                                                                                                                    aTTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGTCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGATAAATCACACAGACCTTTTTTTCCTAGCCTCTGCAAccc                                                                             >  2:202112/1‑273 (MQ=255)
                                                                                                                                                                                     tttCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAAcggcg                                                                                                                                                                                                          >  6:181429/1‑149 (MQ=255)
                                                                                                                                                                                       tCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGt                                                                                                                                                                                                        <  6:385308/149‑1 (MQ=255)
                                                                                                                                                                                        ccGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGc                                                                                                                                                                                                     <  4:521224/151‑1 (MQ=255)
                                                                                                                                                                                                 tCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCg                                                                                                                                                                                              <  5:265078/149‑1 (MQ=255)
                                                                                                                                                                                                  cACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCgcgc                                                                                                                                                                                                                              <  6:133918/116‑1 (MQ=255)
                                                                                                                                                                                                  cACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCgcgc                                                                                                                                                                                                                              >  5:133918/1‑116 (MQ=255)
                                                                                                                                                                                                    ccAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTa                                                                                                                                                                                           >  6:280326/1‑149 (MQ=255)
                                                                                                                                                                                                     cAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTAt                                                                                                                                                                                          >  6:338461/1‑149 (MQ=255)
                                                                                                                                                                                                           cTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATc                                                                                                                                                                                    <  8:287012/149‑1 (MQ=255)
                                                                                                                                                                                                             gTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAACAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGCGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGCTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCCCTTGCTAACTCACCCAGACCTTTTTGGCGTATCGCTCGCAACGCGGCATCATTAAATTCGAGCGTGGc                                                     >  4:433383/1‑274 (MQ=255)
                                                                                                                                                                                                                   ccTGAGCGACCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGc                                                                                                                                                                            >  5:408933/1‑149 (MQ=255)
                                                                                                                                                                                                                        gcgagcACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAg                                                                                                                                                                                                                                                    <  8:58756/72‑1 (MQ=255)
                                                                                                                                                                                                                        gcgagcACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAg                                                                                                                                                                                                                                                    >  7:58756/1‑72 (MQ=255)
                                                                                                                                                                                                                                gTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACgg                                                                                                                                                               >  6:178949/1‑149 (MQ=255)
                                                                                                                                                                                                                                gTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACgg                                                                                                                                                               >  6:376693/1‑149 (MQ=255)
                                                                                                                                                                                                                                      cAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCATGCAGAAAAAGTATCTGCCGAGATATGAATTGCGCAAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGCCGCTAAACTCCGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGATTttgcct                            <  1:140478/274‑1 (MQ=255)
                                                                                                                                                                                                                                      cAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAgt                                                                                                                                                         <  5:160639/149‑1 (MQ=255)
                                                                                                                                                                                                                                        gTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAgttc                                                                                                                                                       >  7:335393/1‑149 (MQ=255)
                                                                                                                                                                                                                                          aCGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAgttcgt                                                                                                                                                     <  5:338461/149‑1 (MQ=255)
                                                                                                                                                                                                                                           cGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAgttcgtt                                                                                                                                                    >  6:14914/1‑149 (MQ=255)
                                                                                                                                                                                                                                              tGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCATCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGGTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCGGGTGTTTCACTTGCCAAATAACACAGACCCTTTTTGCGGAGGGCCAGCCAAGCGGCATCATTAAATCCGAGCGTGGGGAACCGCGTGCGCCCGATTTTGCCTTGCCggacc                   >  4:372606/1‑271 (MQ=255)
                                                                                                                                                                                                                                              tGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAATCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGAGTTGGCCTTGCCg                       >  1:117589/1‑271 (MQ=255)
                                                                                                                                                                                                                                              tGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCTCGCCCGATTTTGCCTTGCCGTAAc                   >  2:193284/1‑275 (MQ=255)
                                                                                                                                                                                                                                                  aTAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGt                                                                                                                                             >  6:157507/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     aaCCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGc                                                                                                                                          <  5:178949/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     aaCCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGc                                                                                                                                          <  6:355870/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        cTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTaa                                                                                                                                       >  8:250133/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        cTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTaa                                                                                                                                       >  6:31146/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              gttACCGGTTCCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTa                                                                                                                                                                                           <  8:287763/90‑1 (MQ=255)
                                                                                                                                                                                                                                                              gttACCGGTTCCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTa                                                                                                                                                                                           >  7:287763/2‑91 (MQ=255)
                                                                                                                                                                                                                                                              gttACCGGTTCCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTg                                                                                                                                 >  7:328998/2‑149 (MQ=255)
                                                                                                                                                                                                                                                                tACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGgt                                                                                                                               <  5:376693/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                tACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGCCGCCAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTCCGTAGCCCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCCAACTGGTCGCGCCAGCATTTGCCTTGCCGTAACGCCCCCAACCCCcgtt   >  4:105222/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                 aCCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGATTTTGCCTTGCCGTAACGCCGCCAGCGCCtgttg  <  2:402128/273‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ggAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGttt                                                                                                                           <  7:264656/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCCACCGCt                                                                                                                                                                           >  4:481150/1‑100 (MQ=255)
                                                                                                                                                                                                                                                                       aCCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCAc                                                                                                                        <  8:335393/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                          
CCGAGATTTTTAAACGGCGTAAACAGCGCGGAGCGGGTGCCGATCACAATCGCCGCCTCACCATTTTTCGCTTTCAGCCACGCCGAAAGACGCTCGCTGTCGTTCAGGCCGGAATGCAGAACTTCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGATTTTGCCTTGCCGTAACGCCGCCAGCGCCTGTTG  >  NC_000913/4125867‑4126396

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