Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,432,995 C→T 8.8% intergenic (‑53/+17) 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,9950CT 8.8% 202.7 / 9.9 68intergenic (‑53/+17)folC/accDbifunctional folylpolyglutamate synthase/ dihydrofolate synthase/acetyl‑CoA carboxylase, beta (carboxyltransferase) subunit
Reads supporting (aligned to +/- strand):  ref base C (32/30);  new base T (3/3);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.20e-01

AAGGCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGGCGCATTTCCGGACGACG  >  NC_000913/2432847‑2433143
                                                                                                                                                    |                                                                                                                                                    
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)
            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                                                                                                                                                    >  1:393429/1‑128 (MQ=255)
                       tgtgCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgccct                                                                                                                                                    <  2:393429/128‑1 (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                                                                                                                                             >  4:608690/1‑124 (MQ=255)
                                 ttCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGGGCCACTGGCTCTGCCCttt                                                                                                                                             <  3:608690/125‑2 (MQ=255)
                                 ttCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCt                                                                                                                     <  1:541683/149‑1 (MQ=255)
                                 ttCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCt                                                                                                                     <  3:915606/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                                                                                                                                       >  2:49238/1‑121 (MQ=255)
                                           gAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatcagtt                                                                                                                                       <  1:49238/121‑1 (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                                                                                                                                            >  3:415853‑M2/5‑89 (MQ=255)
                                                                      gcttGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttat                                                                                                                                            <  4:415853‑M2/85‑1 (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                                                                         <  3:1698837/142‑1 (MQ=255)
                                                                                    tttGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCa                                                                         >  4:1698837/1‑142 (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                                                                      <  4:1518802/126‑1 (MQ=255)
                                                                                                       gCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGa                                                                      >  3:1518802/1‑126 (MQ=255)
                                                                                                          gATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGc                                                                                          <  2:685102/103‑1 (MQ=255)
                                                                                                          gATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGc                                                                                          >  1:685102/1‑103 (MQ=255)
                                                                                                            ttCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatcagttatc                                                                                                                                    <  3:917799/59‑1 (MQ=255)
                                                                                                            ttCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatcagttatc                                                                                                                                    >  4:917799/1‑59 (MQ=255)
                                                                                                            ttCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTg                                                            >  1:943433/1‑131 (MQ=255)
                                                                                                            ttCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTg                                                            <  2:943433/131‑1 (MQ=255)
                                                                                                               gTTATGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCa                                       <  1:1692484/149‑1 (MQ=255)
                                                                                                                       gagaATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCggg                                                                                                       >  2:1086759/1‑77 (MQ=255)
                                                                                                                       gagaATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCggg                                                                                                       <  1:1086759/77‑1 (MQ=255)
                                                                                                                       gagaATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATg                                  >  1:579629/1‑146 (MQ=255)
                                                                                                                       gagaATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATg                                  <  2:579629/146‑1 (MQ=255)
                                                                                                                         gaATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCgggg                                                                                                      <  2:363337/76‑1 (MQ=255)
                                                                                                                         gaATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCgggg                                                                                                      >  1:363337/1‑76 (MQ=255)
                                                                                                                              gCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATtcat                                                  <  2:506963/123‑1 (MQ=255)
                                                                                                                              gCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATtcat                                                  >  1:506963/1‑123 (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)
                                                                                                                                                   ccctgccctTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGaa                                    <  2:1410240/116‑1 (MQ=255)
                                                                                                                                                   ccctgccctTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGaa                                    >  1:1410240/1‑116 (MQ=255)
                                                                                                                                                   ccccgccctTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGGCGCATTTCcggacgac   <  2:736104/149‑1 (MQ=255)
                                                                                                                                                    cctgccctTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGCTTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGGCGCATTTCcggacgacg  >  2:1708808/1‑149 (MQ=255)
                                                                                                                                                    cctgccctTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGaa                                    >  3:1652004/1‑115 (MQ=255)
                                                                                                                                                    cctgccctTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGaa                                    <  4:1652004/115‑1 (MQ=255)
                                                                                                                                                    cctgccctTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGGCGCATTTCcggacgacg  >  3:908109/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAGGCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGGCGCATTTCCGGACGACG  >  NC_000913/2432847‑2433143

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: