Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,126,129 A→T 7.2% S228T (TCC→ACC)  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,1290AT 7.2% 272.1 / 7.6 79S228T (TCC→ACC) priAPrimosome factor n' (replication factor Y)
Reads supporting (aligned to +/- strand):  ref base A (40/33);  new base T (3/3);  total (43/36)
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.20e-01

ATGACAATCACGCCGAGATTTTTAAACGGCGTAAACAGCGCGGAGCGGGTGCCGATCACAATCGCCGCCTCACCATTTTTCGCTTTCAGCCACGCCGAAAGACGCTCGCTGTCGTTCAGGCCGGAATGCAGAACTTCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGATTTTGCCTTGCCGTAACGCCGCCAGCGCCTGTTG  >  NC_000913/4125855‑4126396
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                           
aTGCCAATCACGCCGAGATTTTTAAACGGCGTAAACAGCGCGGAGCGGGTGCCGATCACAATCGCCGCCTCACCATTTTTCGCTTTCAGCCACGCCGAAAGACGCTCGCTGTCGTTCAGGCCGGAATGCAGAACTTCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGa                                                                                                                                                                                                                                                                             <  2:139518/275‑1 (MQ=255)
aTGACAATCACGCCGAGATTTTTAAACGGCGTAAACAGCGCGGAGCGGGTGCCGATCACAATCGCCGCCTCACCATTTTTCGCTTTCAGCCACGCCGAAAGACGCTCGCTGTCGTTCAGGCCGGAATGCAGAACTTCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGa                                                                                                                                                                                                                                                                             <  2:292692/275‑1 (MQ=255)
            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)
                                                                                                                              tGCAGAACTTCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGa                                                                                                                                                                                                                                                                             >  7:153891/1‑149 (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                                                                                                                                                                                                                                           <  7:361641/149‑1 (MQ=255)
                                                                                                                                                                gAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGaaaaa                                                                                                                                                                                                                                           <  5:14090/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                                                                                                                                                                                                                                       >  7:136037/1‑139 (MQ=255)
                                                                                                                                                                              gTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTAt                                                                                                                                                                                                                                       <  8:136037/139‑1 (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                                                                                                                                                                                                                                                    >  7:58756/1‑72 (MQ=255)
                                                                                                                                                                                                                                    gcgagcACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAg                                                                                                                                                                                                                                                    <  8:58756/72‑1 (MQ=255)
                                                                                                                                                                                                                                            gTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACgg                                                                                                                                                               >  6:376693/1‑149 (MQ=255)
                                                                                                                                                                                                                                            gTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACgg                                                                                                                                                               >  6:178949/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                                                                                                                                                                                           >  7:287763/2‑91 (MQ=255)
                                                                                                                                                                                                                                                                          gttACCGGTTCCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTa                                                                                                                                                                                           <  8:287763/90‑1 (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)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                           
ATGACAATCACGCCGAGATTTTTAAACGGCGTAAACAGCGCGGAGCGGGTGCCGATCACAATCGCCGCCTCACCATTTTTCGCTTTCAGCCACGCCGAAAGACGCTCGCTGTCGTTCAGGCCGGAATGCAGAACTTCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGATTTTGCCTTGCCGTAACGCCGCCAGCGCCTGTTG  >  NC_000913/4125855‑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: