Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,126,130 A→T 6.6% 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.6% 247.8 / 7.4 97G227G (GGT→GGA) priAPrimosome factor n' (replication factor Y)
Reads supporting (aligned to +/- strand):  ref base A (46/43);  new base T (2/4);  total (50/47)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.29e-01

GCAGAACTTCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACT  >  NC_000913/4125982‑4126279
                                                                                                                                                    |                                                                                                                                                     
gCAGAACTTCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAAc                                                                                                                                                      >  2:1439600/1‑150 (MQ=255)
  aGAACTTCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCg                                                                                                                                                    <  2:718474/150‑1 (MQ=255)
      cTTCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTTTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTaa                                                                                                                                                <  1:139360/150‑1 (MQ=255)
      cTTCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTaa                                                                                                                                                >  2:1497918/1‑150 (MQ=255)
      cTTCCACGGGGGCATTAAAACGGTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTaa                                                                                                                                                <  2:187417/150‑1 (MQ=255)
         ccACGTGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGTTCCGGTaaaac                                                                                                                                             <  2:1122503/150‑4 (MQ=255)
         ccACGGGGGCATTTAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAACCCTCCGTTTTTCCGGAACCGGTAAcgc                                                                                                                                             >  2:1654397/1‑150 (MQ=255)
         ccACGGGGGCATTAACACGTGCACGAAAACGGGCGATAGTTTGCGGTGTCAGTGCGATTTACGGCACCATCACCAGCTTATTTTTTCCCTTAGCTAGCACTTTTTCCAGTACGATGAGATAAACCTCCGTTTTACCGGTTCCGGTaaaac                                                                                                                                             <  2:284593/150‑4 (MQ=255)
         ccACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGTTCCGGTaaaac                                                                                                                                             <  2:785540/150‑4 (MQ=255)
          cACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGTTCCGGTaaaacg                                                                                                                                            <  2:758297/150‑5 (MQ=255)
            cgtggtcATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTACAGTACGCTTAGATAAACCTCCGTTTTACCGGAACCGGTAAcgcccg                                                                                                                                          <  2:1172354/144‑1 (MQ=255)
               gggCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGGCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGGTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCca                                                                                                                                       <  2:1646332/150‑1 (MQ=255)
                ggCATTAAAACGTTCTCGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCcag                                                                                                                                      <  2:818495/150‑1 (MQ=255)
                 gCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCcagc                                                                                                                                     <  2:537613/150‑1 (MQ=255)
                  cATTAAAACGTTCACGAATACGGGCGATAGTTTGCGGTGTCTGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCcagca                                                                                                                                    <  1:1539256/150‑1 (MQ=255)
                  cATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCcagca                                                                                                                                    <  1:1053308/150‑1 (MQ=255)
                     tAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGcc                                                                                                                                 <  2:132153/150‑1 (MQ=255)
                      aaaaCGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCa                                                                                                                                <  2:81336/150‑1 (MQ=255)
                      aaaaCGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCa                                                                                                                                >  2:518470/1‑150 (MQ=255)
                         aCGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGGGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGc                                                                                                                             <  1:1213742/150‑1 (MQ=255)
                         aCGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGc                                                                                                                             <  2:1653778/150‑1 (MQ=255)
                              cACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGaaa                                                                                                                        >  2:935715/1‑150 (MQ=255)
                               aCGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGCGCACGTTTTCCAGGACGCTGAGCAAAACCTCCGTATTTCCGGAACCGGTAACGCCCGCCAGCCGACAGGCAGaaaa                                                                                                                       >  2:1212721/1‑150 (MQ=255)
                               aCGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTACCGCCCGCCAGCAGCCAGGCAGaaaa                                                                                                                       >  1:630136/1‑150 (MQ=255)
                               aCGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGaaaa                                                                                                                       >  1:1077084/1‑150 (MQ=255)
                                  aaaaCGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCCGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGt                                                                                                                    >  2:292249/1‑150 (MQ=255)
                                  aaaaCGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGACCCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGt                                                                                                                    >  2:481790/1‑150 (MQ=255)
                                   aaaCGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTa                                                                                                                   <  1:762097/150‑1 (MQ=255)
                                          cGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGCGCGAGCACGTTTTCCAGTACGCTGAGATACACCTCCGTTTTACCGGAACCGGTAACGACCGCCAGAAGCCAGGCAGCAACAATATCttgac                                                                                                            >  2:352447/1‑146 (MQ=255)
                                           gATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCcgc                                                                                                           >  2:151634/1‑150 (MQ=255)
                                                tttGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTAt                                                                                                      >  2:374094/1‑150 (MQ=255)
                                                tttGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTAt                                                                                                      >  1:714827/1‑150 (MQ=255)
                                                    cGGTGTCAGGCCGATTTCCGGCTCCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAAt                                                                                                  <  2:1302562/150‑1 (MQ=255)
                                                            ggCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTACCGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCaa                                                                                          >  1:1520849/1‑150 (MQ=255)
                                                            ggCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCaa                                                                                          >  2:1395245/1‑150 (MQ=255)
                                                             gCCGATTTCCGGCACCTTCACCAGCGCCTGTTAGCCCTGAGCGAGCTCGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACAGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAATAGTATCTGCCGCGCTATGAATTGCGCCAAc                                                                                         <  1:267266/150‑1 (MQ=255)
                                                              ccGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAAcg                                                                                        >  1:1292482/1‑150 (MQ=255)
                                                               cGATTTCCGGCACCATCACCAGCGCCTGTTTGGGCTGAGCGAGCACGTTTTCCCGGAGGCTGAGATAAACCTCCGTTTTTGCGGAGCCGGTAACGCCCGCCAGCTGCCAGGCAGCAGAAGTATCTGCCGCGCGATGAATTGAGCAATcgg                                                                                       >  1:994803/1‑150 (MQ=255)
                                                               cGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAAcgg                                                                                       >  2:1647509/1‑150 (MQ=255)
                                                                      cGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGc                                                                                <  2:217613/150‑1 (MQ=255)
                                                                      cGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGGACGCTGAGATTAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCCGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGGGGc                                                                                >  1:981886/1‑150 (MQ=255)
                                                                       ggCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGcc                                                                               <  1:535673/150‑1 (MQ=255)
                                                                       ggCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGcc                                                                               >  1:1462127/1‑150 (MQ=255)
                                                                       ggCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGcc                                                                               <  1:151634/150‑1 (MQ=255)
                                                                        gCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTAACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCt                                                                              <  2:1714424/150‑1 (MQ=255)
                                                                         cACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTg                                                                             <  2:5840/150‑1 (MQ=255)
                                                                            cATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTc                                                                          <  2:1315042/150‑1 (MQ=255)
                                                                              tCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCgg                                                                        <  2:1761812/150‑1 (MQ=255)
                                                                              tCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCgg                                                                        >  1:740465/1‑150 (MQ=255)
                                                                                aCCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTa                                                                      <  2:1555935/150‑1 (MQ=255)
                                                                                 ccAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTAt                                                                     >  1:1518852/1‑150 (MQ=255)
                                                                                 ccAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTAt                                                                     <  2:493874/150‑1 (MQ=255)
                                                                                  cAGCGACTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTAACTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGtttt                                                                    >  2:260305/1‑150 (MQ=255)
                                                                                   cgcgcCTGTTTGCCCGGAGCGAGCACGTTCTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTc                                                                   <  1:1585389/149‑1 (MQ=255)
                                                                                        cTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCg                                                              <  1:1569000/150‑1 (MQ=255)
                                                                                        cTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCg                                                              >  2:1020145/1‑150 (MQ=255)
                                                                                           tttGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATACACTCGCaa                                                           >  2:352074/1‑150 (MQ=255)
                                                                                             tGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAAcc                                                         >  1:1697566/1‑150 (MQ=255)
                                                                                              gCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCg                                                        <  2:689299/150‑1 (MQ=255)
                                                                                               cccTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGc                                                       <  2:582980/150‑1 (MQ=255)
                                                                                                 cTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTc                                                     <  1:1396631/150‑1 (MQ=255)
                                                                                                 cTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTc                                                     >  1:900117/1‑150 (MQ=255)
                                                                                                        agcACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGaa                                              >  1:1463592/1‑150 (MQ=255)
                                                                                                         gcACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGaaa                                             >  1:1514496/1‑150 (MQ=255)
                                                                                                            cGTTTTCCAGTACGCTGAGATAAACCTCCGATTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGAATTCAATCGCAACCGCTCACCAGAAAccc                                          >  2:491073/1‑148 (MQ=255)
                                                                                                               tttCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCAt                                       >  1:200628/1‑150 (MQ=255)
                                                                                                                ttCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATa                                      <  1:1770128/150‑1 (MQ=255)
                                                                                                                ttCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATa                                      <  2:1583927/150‑1 (MQ=255)
                                                                                                                   cAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAgtt                                   >  1:1489751/1‑150 (MQ=255)
                                                                                                                   cAGTACGCTGAGATAAACCTCCGATTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAgtt                                   <  2:552961/150‑1 (MQ=255)
                                                                                                                     gTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAgttcg                                 <  1:1178125/150‑1 (MQ=255)
                                                                                                                      tACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAgttcgt                                <  1:1682890/150‑1 (MQ=255)
                                                                                                                        cGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAgttcgttc                              >  2:1040535/1‑150 (MQ=255)
                                                                                                                        cGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAgttcgttc                              <  1:1605722/150‑1 (MQ=255)
                                                                                                                        cGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTAACTGCCGATCTATGAAATGCGCCAACGGCGGTGGCCTGTTCGGTTTTCAATCGCAACCGATCACCAGAACCGGCCTAgttcgttc                              >  2:111970/1‑150 (MQ=255)
                                                                                                                         gCTGAGCTCCACCTCAGGTTTACCGGAACCTGTAACGCCCGCCAGCAGCCTGTCAGAAAACGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCACCCGCTCACCAGAAACTGCATAGttcgttcg                             <  2:318252/150‑1 (MQ=255)
                                                                                                                         gCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGttcgttcg                             <  1:518470/150‑1 (MQ=255)
                                                                                                                         gCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGttcgttcg                             <  1:125297/150‑1 (MQ=255)
                                                                                                                          cTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAGCGGCATAGTTCGTTCGc                            >  1:1046401/1‑150 (MQ=255)
                                                                                                                          cTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTAGTCAATCGCAACCGCTCACCAGAAACGGCCTAGTTCGTTCGc                            >  1:791821/1‑150 (MQ=255)
                                                                                                                            gagaTAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCa                          >  1:988996/1‑150 (MQ=255)
                                                                                                                                tAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCg                      >  1:1353816/1‑150 (MQ=255)
                                                                                                                                     cTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTaaa                 >  2:1449552/1‑150 (MQ=255)
                                                                                                                                     cTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTaaa                 <  2:1294698/150‑1 (MQ=255)
                                                                                                                                      tCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAAc                <  2:1077084/150‑1 (MQ=255)
                                                                                                                                       ccGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAAct               <  2:1462127/150‑1 (MQ=255)
                                                                                                                                       ccGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAAct               >  1:823923/1‑150 (MQ=255)
                                                                                                                                         gcgttACCGGTTCCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAACctct             >  1:1775437/4‑150 (MQ=255)
                                                                                                                                         gcgttACCGGTTCCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAActct             >  2:459505/4‑150 (MQ=255)
                                                                                                                                         gTTTTACCGGAACCGGTAACGCCCGCCACCACCCAGGCAGAAAAAGTATCTGACGCGCTATGAATTGCGCCACCGGTGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAGCTGCATAGTTCGTTCGCCAGACGCTAAAatct             >  2:1767139/1‑150 (MQ=255)
                                                                                                                                               ccGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGtt       >  1:1771246/1‑150 (MQ=255)
                                                                                                                                               ccGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGtt       >  1:1809238/1‑150 (MQ=255)
                                                                                                                                                cGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGttt      >  1:1807294/1‑150 (MQ=255)
                                                                                                                                                cGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGttt      >  1:549500/1‑150 (MQ=255)
                                                                                                                                                 ggAACCGTGAACGCCCGCCAGCAGCCAGTCAGAAAATGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTc     <  1:337674/150‑1 (MQ=255)
                                                                                                                                                 ggAACCGGTAACGCCCGCCAGCAGCCCGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTc     >  2:1412736/1‑150 (MQ=255)
                                                                                                                                                 ggAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCCCCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTGTc     >  2:851034/1‑150 (MQ=255)
                                                                                                                                                    aCCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGGGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACCGCATAGTTCGTACTCCAGTCGCTAAACTCTGGTGATTCACt  >  1:96170/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCAGAACTTCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACT  >  NC_000913/4125982‑4126279

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