Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,956,003 G→A D225N (GAT→AAT)  ilvA → l‑threonine dehydratase, biosynthetic; also known as threonine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,956,0030GA96.3% 1340.9 / 25.2 406D225N (GAT→AAT) ilvAl‑threonine dehydratase, biosynthetic; also known as threonine deaminase
Reads supporting (aligned to +/- strand):  ref base G (8/7);  new base A (206/185);  total (214/192)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

Only 100 of 407 total aligned reads displayed.

GCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATG  >  NC_000913/3955855‑3956152
                                                                                                                                                    |                                                                                                                                                     
  agGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGc                                                                                                                                                   <  8:198832/151‑1 (MQ=255)
  agGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGc                                                                                                                                                   >  4:1166354/1‑151 (MQ=255)
    gACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCgg                                                                                                                                                 <  6:746716/151‑1 (MQ=255)
    gACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCgg                                                                                                                                                 <  3:515350/151‑1 (MQ=255)
     aCGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCggg                                                                                                                                                <  8:1460502/151‑1 (MQ=255)
       gCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTc                                                                                                                                              <  6:894953/151‑1 (MQ=255)
            tctcGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCg                                                                                                                                         >  8:792554/1‑151 (MQ=255)
              tcGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGt                                                                                                                                       >  6:1170522/1‑151 (MQ=255)
               cGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCgg                                                                                                                                        <  3:1043419/149‑1 (MQ=255)
               cGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGtt                                                                                                                                      >  3:962115/1‑151 (MQ=255)
                gACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTg                                                                                                                                     >  5:887338/1‑151 (MQ=255)
                 aCCGCGTATTTGTGCCAGTCGGCGGCGGCGGGCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGAGCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGa                                                                                                                                    <  6:135972/151‑1 (MQ=255)
                   cgcgTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATc                                                                                                                                  >  1:910069/1‑151 (MQ=255)
                    gcgTATTTGTGCCAGTCGGCGGCTGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCTTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCt                                                                                                                                 <  7:1229853/151‑1 (MQ=255)
                     cgTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTg                                                                                                                                >  1:798107/1‑151 (MQ=255)
                     cgTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTg                                                                                                                                <  4:1069649/151‑1 (MQ=255)
                      gTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGCTCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGGCATCCGGTTGCTCTGc                                                                                                                               >  2:1427876/1‑151 (MQ=255)
                      gTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGc                                                                                                                               <  6:1051124/151‑1 (MQ=255)
                        aTTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCcg                                                                                                                             <  8:561548/151‑1 (MQ=255)
                         tttGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCcgc                                                                                                                            >  7:870133/1‑151 (MQ=255)
                         tttGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCcgc                                                                                                                            >  5:240696/1‑151 (MQ=255)
                             tgCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGt                                                                                                                        >  1:1024016/1‑151 (MQ=255)
                              gCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTa                                                                                                                       <  3:863753/151‑1 (MQ=255)
                               ccAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCAGCCGGTTGATCTGCCGCGCGt                                                                                                                        >  3:1164945/1‑149 (MQ=255)
                                cAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGCCTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAgg                                                                                                                     >  6:614637/1‑151 (MQ=255)
                                cAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAgg                                                                                                                     <  5:444368/151‑1 (MQ=255)
                                cAGTCGGCGGCGGCGGTCTGGCTGCTGGAGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAgg                                                                                                                     >  8:1403356/1‑151 (MQ=255)
                                  gTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGc                                                                                                                   >  5:185395/1‑151 (MQ=255)
                                   gcggcggcggcggGCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCt                                                                                                                  <  2:1587933/150‑1 (MQ=255)
                                    cggcggcggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCt                                                                                                                  <  3:706121/150‑1 (MQ=255)
                                       cggcggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTAttt                                                                                                              <  6:317973/151‑1 (MQ=255)
                                        ggcggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTg                                                                                                             <  5:1626280/151‑1 (MQ=255)
                                        ggcggcggTATCCCTGCTGCCGTGGAGGAGCTCATCAAACAACTGATCCCGCAAATCAAAGTGATCGACGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTAGATCTGCCGCGAGTAGGGCTAtttc                                                                                                             >  6:1656568/1‑150 (MQ=255)
                                         gcggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGc                                                                                                            >  8:1556048/1‑151 (MQ=255)
                                           ggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTg                                                                                                          <  2:314676/151‑1 (MQ=255)
                                            gcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTg                                                                                                          >  8:233685/1‑150 (MQ=255)
                                              ggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGa                                                                                                         >  3:182632/1‑149 (MQ=255)
                                                 cTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCg                                                                                                    >  2:1074486/1‑151 (MQ=255)
                                                   ggCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGc                                                                                                     <  7:1080694/148‑1 (MQ=255)
                                                   ggCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTa                                                                                                  <  3:1521712/151‑1 (MQ=255)
                                                   ggCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTa                                                                                                  <  7:1457329/151‑1 (MQ=255)
                                                     ctgctgGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCTGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGc                                                                                                >  5:619933/1‑151 (MQ=255)
                                                     ctgctgGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGc                                                                                                >  4:1269772/1‑151 (MQ=255)
                                                      tgctgGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCg                                                                                               >  1:344081/1‑151 (MQ=255)
                                                       gcttGCGTGGCGGTGGTGTTGTTTCATGTGATGCGGCAAATCAATGTGATCGCCTTAGAAGATGAAGACTGCGCCTGCCTGAAAGCAGCGCTGAATGAGGGTGATGGGGTTGTTCTGCGGCGCGTAGGGCTATTTGCTGAAGGGGTAggg                                                                                               <  6:1698424/150‑3 (MQ=255)
                                                         tggcgtggcgGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTACCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGt                                                                                             >  3:1362978/1‑150 (MQ=255)
                                                             gtggcgGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAAc                                                                                        <  5:1524870/151‑1 (MQ=255)
                                                               ggcgGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGc                                                                                      <  2:93750/151‑1 (MQ=255)
                                                                   gTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCg                                                                                  >  2:375080/1‑151 (MQ=255)
                                                                    tGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCgg                                                                                 >  5:1205840/1‑151 (MQ=255)
                                                                    tGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCgg                                                                                 >  1:408420/1‑151 (MQ=255)
                                                                    tGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCgg                                                                                 >  7:389954/1‑151 (MQ=255)
                                                                    ggCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCgg                                                                                 <  1:1408167/150‑1 (MQ=255)
                                                                      cTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTg                                                                               >  8:11719/1‑151 (MQ=255)
                                                                        gATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGAc                                                                             <  3:1236591/151‑1 (MQ=255)
                                                                        gATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGAc                                                                             <  5:329758/151‑1 (MQ=255)
                                                                        gATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGAc                                                                             >  4:1675003/1‑151 (MQ=255)
                                                                           cAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGa                                                                           <  6:377980/150‑1 (MQ=255)
                                                                             aacaacTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCTGAAGACTCCGCCTGCATGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAAc                                                                        <  1:124638/151‑1 (MQ=255)
                                                                             aacaacTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAAc                                                                        >  4:359999/1‑151 (MQ=255)
                                                                             aacaacTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAAc                                                                        <  4:804421/151‑1 (MQ=255)
                                                                              acaacTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAAcc                                                                       >  6:1475234/1‑151 (MQ=255)
                                                                                 acTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTc                                                                    <  7:775061/151‑1 (MQ=255)
                                                                                        ccGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTAt                                                             <  8:294243/151‑1 (MQ=255)
                                                                                                cAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAgg                                                       >  1:627443/1‑149 (MQ=255)
                                                                                                cAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAg                                                     <  6:82709/151‑1 (MQ=255)
                                                                                                 aaaGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCAGTTTATGCCAGGAGt                                                    >  5:1239006/1‑151 (MQ=255)
                                                                                                   agggATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGGAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTAt                                                  <  6:945963/148‑1 (MQ=255)
                                                                                                      gATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTAtctc                                               <  3:338/151‑1 (MQ=255)
                                                                                                      gATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTAtctc                                               >  1:1496844/1‑151 (MQ=255)
                                                                                                      gATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTAtctc                                               >  2:440185/1‑151 (MQ=255)
                                                                                                      gATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTAtct                                                <  7:705137/150‑1 (MQ=255)
                                                                                                        tCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGGAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTcg                                              <  8:1233274/150‑1 (MQ=255)
                                                                                                           ccGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTcgacga                                          <  8:1025347/151‑1 (MQ=255)
                                                                                                           ccGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTcgacg                                           >  5:747364/1‑150 (MQ=255)
                                                                                                           ccGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTcgacg                                           >  5:1069200/1‑150 (MQ=255)
                                                                                                              tAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGAcat                                       >  6:475736/1‑151 (MQ=255)
                                                                                                                 aaGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGAcatcat                                    >  1:1064851/1‑151 (MQ=255)
                                                                                                                   gCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACatcatc                                   >  2:1227397/1‑150 (MQ=255)
                                                                                                                   ccGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACAtcatca                                  <  4:261773/150‑1 (MQ=255)
                                                                                                                     ggAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCAcc                                >  4:627502/1‑151 (MQ=255)
                                                                                                                     ggAAGACTCCGACTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCAcc                                >  4:1509764/1‑151 (MQ=255)
                                                                                                                        aGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGCCATCATCACCGTc                             >  1:766675/1‑151 (MQ=255)
                                                                                                                            tCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGAt                          >  5:550050/1‑150 (MQ=255)
                                                                                                                             ccGCCTGCCTGAAAGCAGCGCTGAATGCGTGTCATCCGGTTGATCTGCCGCGCGTAGTGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGACAGGAGTATCTCGACGACATCATCACCGTCGAt                          >  4:29497/1‑149 (MQ=255)
                                                                                                                             ccGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGAt                          >  1:129759/1‑149 (MQ=255)
                                                                                                                             ccGCCTGCCGGAAAGCAGCGCGGAATGCGGGTCATCCGGGTGATCTGCCGCGCGTAGGGCTATTTGCGGAAGGCGTAGCGGGAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATa                         <  1:976207/150‑1 (MQ=255)
                                                                                                                              cGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgc                       >  7:347630/1‑151 (MQ=255)
                                                                                                                              cGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgc                       >  2:1311075/1‑151 (MQ=255)
                                                                                                                                 ctgcctgAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgcgat                    >  4:528322/1‑151 (MQ=255)
                                                                                                                                   gcctgAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgcgatgc                  >  8:90746/1‑151 (MQ=255)
                                                                                                                                     ctgAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgcgatgcg                 >  3:451114/1‑150 (MQ=255)
                                                                                                                                       gAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgcgatgcga                <  7:214813/149‑1 (MQ=255)
                                                                                                                                           gcagcGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCATCGTCGATAGCGATGCGATCtgtg          <  8:863874/151‑1 (MQ=255)
                                                                                                                                           gcagcGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCtgtg          >  6:1020281/1‑151 (MQ=255)
                                                                                                                                            cagcGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTgc         >  7:463321/1‑151 (MQ=255)
                                                                                                                                            cagcGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTgc         <  1:1683998/151‑1 (MQ=255)
                                                                                                                                              gcgcTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGCCGACATCATCACCGTCGATAGCGATGCGATCTGTgcgg       >  7:614212/1‑151 (MQ=255)
                                                                                                                                                gcTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTgcggcg     >  7:1595912/1‑151 (MQ=255)
                                                                                                                                                gcTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTgcggcg     <  2:805936/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATG  >  NC_000913/3955855‑3956152

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