Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,432,990 G→C 7.0% intergenic (‑48/+22) folC ← / ← accD bifunctional folylpolyglutamate synthase/ dihydrofolate synthase/acetyl‑CoA carboxylase, beta (carboxyltransferase) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,432,9900GC 7.0% 304.8 / 11.1 86intergenic (‑48/+22)folC/accDbifunctional folylpolyglutamate synthase/ dihydrofolate synthase/acetyl‑CoA carboxylase, beta (carboxyltransferase) subunit
Reads supporting (aligned to +/- strand):  ref base G (40/40);  new base C (3/3);  total (43/43)
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.84e-01

AAGGCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGGCGCATT  >  NC_000913/2432847‑2433132
                                                                                                                                               |                                                                                                                                              
aaGGCCGAGTTCGTTATTTTTACTGTGCAGTTTTTCAAGTTAAGAAAGCCACGAAGCCAGAGGGAACGCGGCTTGAGGAGGGCGTTTGATAACCATGGTACCGGCGGATTCGTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgccc                                                                                                                                          <  3:1369817/149‑1 (MQ=255)
aaGGCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgcc                                                                                                                                           <  3:1124528/149‑1 (MQ=255)
 aGGCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgccc                                                                                                                                          >  2:433822/1‑149 (MQ=255)
 aGGCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgccc                                                                                                                                          >  2:1050968/1‑149 (MQ=255)
 aGGCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgccc                                                                                                                                          >  4:902130/1‑149 (MQ=255)
   gCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgccctg                                                                                                                                        >  2:78652/1‑149 (MQ=255)
     cGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgccctgcc                                                                                                                                      <  1:1050968/149‑1 (MQ=255)
     cGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgccctgcc                                                                                                                                      <  3:902130/149‑1 (MQ=255)
     cGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgccctgcc                                                                                                                                      <  4:276591/149‑1 (MQ=255)
     cGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgccctgcc                                                                                                                                      <  1:433822/149‑1 (MQ=255)
            gATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCTGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgcg                                                                                                                                           <  2:1595707/137‑2 (MQ=255)
            gATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgcc                                                                                                                                           >  1:1595707/1‑137 (MQ=255)
                 ttttACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgccctgcc                                                                                                                                      <  3:1147702/137‑1 (MQ=255)
                 ttttACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgccctgcc                                                                                                                                      >  4:1147702/1‑137 (MQ=255)
                   ttACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatcagttatca                                                                                                                        <  4:739079/149‑1 (MQ=255)
                       tgtgCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgccct                                                                                                                                         <  2:393429/128‑1 (MQ=255)
                       tgtgCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgccct                                                                                                                                         >  1:393429/1‑128 (MQ=255)
                        gtgCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCt                                                                                                                   >  4:915606/1‑149 (MQ=255)
                          gCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCa                                                                                                                 >  1:350758/1‑149 (MQ=255)
                             ggTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCa                                                                                                                 <  2:743982/146‑1 (MQ=255)
                             ggTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCa                                                                                                                 >  1:743982/1‑146 (MQ=255)
                                tttCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTcc                                                                                                           <  1:78652/149‑1 (MQ=255)
                                 ttCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGGGCCACTGGCTCTGCCCttt                                                                                                                                  <  3:608690/125‑2 (MQ=255)
                                 ttCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGGGCCACTGGCTCTGCCCttt                                                                                                                                  >  4:608690/1‑124 (MQ=255)
                                 ttCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCt                                                                                                          <  3:915606/149‑1 (MQ=255)
                                 ttCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCt                                                                                                          <  1:541683/149‑1 (MQ=255)
                                     aGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATc                                                                                                      >  4:1303649/1‑149 (MQ=255)
                                     aGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATc                                                                                                      >  3:1236343/1‑149 (MQ=255)
                                     aGATAAGAAAGACACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATc                                                                                                      >  3:1080035/1‑149 (MQ=255)
                                           gaatGCCACGAAGCGAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatcagtt                                                                                                                            >  3:408293‑M1/5‑121 (MQ=255)
                                           gAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCACAAGGGCAGAGCCAGTGGCCCTGCCCttatcagtt                                                                                                                            >  4:408293/1‑121 (MQ=255)
                                           gAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatcagtt                                                                                                                            <  1:49238/121‑1 (MQ=255)
                                           gAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatcagtt                                                                                                                            >  2:49238/1‑121 (MQ=255)
                                           gAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTAc                                                                                                >  1:1261787/1‑149 (MQ=255)
                                              AGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatcagt                                                                                                                             >  2:855489‑M1/1‑117 (MQ=255)
                                              AGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatcagt                                                                                                                             <  1:855489‑M1/117‑1 (MQ=255)
                                                CCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCgggg                                                                                           >  4:1236343‑M1/1‑149 (MQ=255)
                                                CCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTGACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCgggg                                                                                           >  4:1080035‑M1/1‑149 (MQ=255)
                                                  aCGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGAGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCATGCCTCAGGTTCCTGATCCGGTACCTGGGGt                                                                                         >  3:19048/1‑149 (MQ=255)
                                                     tcgcctctggcttcgtggcttGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTTGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGATcc                                                                                                           >  1:1207500‑M2/22‑128 (MQ=255)
                                                     tcgcctctggcttcgtggcttGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTcc                                                                                                           <  2:1207500‑M2/107‑1 (MQ=255)
                                                       gCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCAc                                                                                    <  1:1159236/149‑1 (MQ=255)
                                                            agGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGc                                                                               <  3:1303649/149‑1 (MQ=255)
                                                                      gcttGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttat                                                                                                                                 <  4:415853‑M2/85‑1 (MQ=255)
                                                                      gcttGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttat                                                                                                                                 >  3:415853‑M2/5‑89 (MQ=255)
                                                                        ttGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgcc                                                                                                                                           >  2:217209‑M2/3‑77 (MQ=255)
                                                                        ttGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgcc                                                                                                                                           <  1:217209‑M2/75‑1 (MQ=255)
                                                                             gagTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCa                                                              <  2:350758/149‑1 (MQ=255)
                                                                                   gTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTg                                                                                                         <  3:927034/100‑1 (MQ=255)
                                                                                   gTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTg                                                                                                         >  4:927034/1‑100 (MQ=255)
                                                                                   gTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTc                                                        <  2:1261787/149‑1 (MQ=255)
                                                                                   gTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTc                                                        >  1:35466/1‑149 (MQ=255)
                                                                                    tttGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCa                                                              >  4:1698837/1‑142 (MQ=255)
                                                                                    tttGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCa                                                              <  3:1698837/142‑1 (MQ=255)
                                                                                           aTCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTgg                                                <  4:19048/149‑1 (MQ=255)
                                                                                                    tCCGCTGATTCGTTATGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATtcat                                       >  2:1692484/1‑149 (MQ=255)
                                                                                                    tCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAg                                                             <  2:524075/127‑1 (MQ=255)
                                                                                                    tCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAg                                                             >  1:524075/1‑127 (MQ=255)
                                                                                                     ccGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATtcatc                                      <  3:1245490/149‑1 (MQ=255)
                                                                                                     ccGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATtcatc                                      <  2:35466/149‑1 (MQ=255)
                                                                                                       gCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGa                                                           >  3:1518802/1‑126 (MQ=255)
                                                                                                       gCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGa                                                           <  4:1518802/126‑1 (MQ=255)
                                                                                                          gATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGc                                                                               <  2:685102/103‑1 (MQ=255)
                                                                                                          gATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGc                                                                               >  1:685102/1‑103 (MQ=255)
                                                                                                            ttCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatcagttatc                                                                                                                         >  4:917799/1‑59 (MQ=255)
                                                                                                            ttCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatcagttatc                                                                                                                         <  3:917799/59‑1 (MQ=255)
                                                                                                            ttCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTg                                                 >  1:943433/1‑131 (MQ=255)
                                                                                                            ttCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTg                                                 <  2:943433/131‑1 (MQ=255)
                                                                                                               gTTATGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCa                            <  1:1692484/149‑1 (MQ=255)
                                                                                                                       gagaATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCggg                                                                                            <  1:1086759/77‑1 (MQ=255)
                                                                                                                       gagaATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCggg                                                                                            >  2:1086759/1‑77 (MQ=255)
                                                                                                                       gagaATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATg                       <  2:579629/146‑1 (MQ=255)
                                                                                                                       gagaATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATg                       >  1:579629/1‑146 (MQ=255)
                                                                                                                         gaATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCgggg                                                                                           >  1:363337/1‑76 (MQ=255)
                                                                                                                         gaATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCgggg                                                                                           <  2:363337/76‑1 (MQ=255)
                                                                                                                              gCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATtcat                                       >  1:506963/1‑123 (MQ=255)
                                                                                                                              gCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATtcat                                       <  2:506963/123‑1 (MQ=255)
                                                                                                                                ataaGGGCAGGGCCACTGGCTCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTc                                >  4:842537/3‑128 (MQ=255)
                                                                                                                                ataaGGGCAGGGCCACTGGCTCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTc                                <  3:842537/126‑1 (MQ=255)
                                                                                                                                aaaaGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTAc                                                                                 >  4:1678816/1‑79 (MQ=255)
                                                                                                                                aaaaGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTAc                                                                                 <  3:1678816/79‑1 (MQ=255)
                                                                                                                                aaaaGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGgcgc                                                   <  1:1161393/109‑1 (MQ=255)
                                                                                                                                aaaaGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGgcgc                                                   >  2:1161393/1‑109 (MQ=255)
                                                                                                                                      gCAGGGCCACTGGCTCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTg                                                 <  3:807641/105‑1 (MQ=255)
                                                                                                                                      gCAGGGCCACTGGCTCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTg                                                 >  4:807641/1‑105 (MQ=255)
                                                                                                                                         gagCCAGTGGCCCCGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGGCGCAtt  >  1:736104/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                              
AAGGCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGGCGCATT  >  NC_000913/2432847‑2433132

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