Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,931,374	A→G	Q165R (CAG→CGG)	purT →	phosphoribosylglycinamide formyltransferase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,931,374	0	A	G	100.0%	246.8 / NA	69	Q165R (CAG→CGG)	purT	phosphoribosylglycinamide formyltransferase 2
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/36);  total (33/36)

TATCCGTCGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGG  >  NC_000913/1931228‑1931521                                                                                                                                                   |                                                                                                                                                    tATCCGTCGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACg                                                                                                                                                 >  1:385368/1‑151 (MQ=255)  aTCCGTCGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACg                                                                                                                                                 <  1:880552/150‑1 (MQ=255)  aTCCGTCGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGt                                                                                                                                                >  1:749318/1‑151 (MQ=255)     cGTCGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGttt                                                                                                                                              >  1:645752/1‑150 (MQ=255)      gTCGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTAt                                                                                                                                            <  1:720595/151‑1 (MQ=255)      gTCGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTAt                                                                                                                                            >  1:615063/1‑151 (MQ=255)        cGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGGGATGAGCCCTTCCGGCAAGGGGCGGACGTTTAtt                                                                                                                                           >  2:623533/1‑150 (MQ=255)         gCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCg                                                                                                                                         <  1:114253/151‑1 (MQ=255)               cggcAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTg                                                                                                                                    >  2:946816/1‑150 (MQ=255)                   agaagaGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCa                                                                                                                                  >  2:812349/1‑148 (MQ=255)                     aagaGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCa                                                                                                                             <  2:114612/151‑1 (MQ=255)                     aagaGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCa                                                                                                                             <  1:935715/151‑1 (MQ=255)                       gagCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAAc                                                                                                                           <  1:882411/151‑1 (MQ=255)                        agctgcagctgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACt                                                                                                                          <  1:586206/151‑1 (MQ=255)                            gcagctgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCt                                                                                                                      <  2:483166/151‑1 (MQ=255)                               gctgcCTACTTCCGCTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGATTGATGACATTGGCTATCCCTACATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCa                                                                                                                    <  2:3432/150‑1 (MQ=255)                               gctgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGGTGAGCTCTTCCAGCAAGGGGCGGACGTTTATTCGGTCTGCAGAGCAACTTGc                                                                                                                       >  2:1029724/1‑147 (MQ=255)                                ctgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAgg                                                                                                                  >  2:720845/1‑151 (MQ=255)                                  gcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGc                                                                                                                 <  1:387652/150‑1 (MQ=255)                                   cccACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCAt                                                                                                               <  1:314840/151‑1 (MQ=255)                                     cacttccacttATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATgg                                                                                                             <  2:956842/151‑1 (MQ=255)                                        ttccacttATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGa                                                                                                            >  1:302944/1‑149 (MQ=255)                                         tccacttATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAg                                                                                                          <  1:752328/150‑1 (MQ=255)                                         tccacttATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGt                                                                                                         >  1:846657/1‑151 (MQ=255)                                                   gTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGc                                                                                                >  1:233272/1‑150 (MQ=255)                                                   gTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGc                                                                                                >  1:319248/1‑150 (MQ=255)                                                   gTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGGAAAACCGGGGATGAGCTCTTCCGGCAAGGGGGGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAg                                                                                                 >  2:950138/1‑149 (MQ=255)                                                     tttGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCaa                                                                                              <  2:481761/150‑1 (MQ=255)                                                           gATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTc                                                                                       <  1:547324/151‑1 (MQ=255)                                                           gATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTc                                                                                       <  2:479540/151‑1 (MQ=255)                                                           gATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTc                                                                                       <  1:376264/151‑1 (MQ=255)                                                            aTAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTcg                                                                                      >  1:565404/1‑151 (MQ=255)                                                                 gAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGcc                                                                                  >  2:907370/1‑150 (MQ=255)                                                                   aaGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGa                                                                               <  1:983016/151‑1 (MQ=255)                                                                     gCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGc                                                                             <  2:359723/151‑1 (MQ=255)                                                                            ccGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGcc                                                                        <  1:64031/149‑1 (MQ=255)                                                                            ccGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCcgc                                                                      <  2:811134/151‑1 (MQ=255)                                                                             cgcgAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCcgcg                                                                     >  1:834093/1‑151 (MQ=255)                                                                                    cGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAAttgt                                                              >  2:504861/1‑151 (MQ=255)                                                                                     ggTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAg                                                                                                 >  1:1030719/1‑115 (MQ=255)                                                                                     ggTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAg                                                                                                 <  2:1030677/115‑1 (MQ=255)                                                                                        tGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGa                                                           >  1:438650/1‑150 (MQ=255)                                                                                         gCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAg                                                         <  2:104514/151‑1 (MQ=255)                                                                                           tGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGc                                                       >  2:291356/1‑151 (MQ=255)                                                                                            gACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGcg                                                      >  1:963016/1‑151 (MQ=255)                                                                                                 tGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGcgtcgt                                                  <  1:105757/150‑1 (MQ=255)                                                                                                  ggCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTa                                                <  2:767795/151‑1 (MQ=255)                                                                                                   gCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTa                                                >  2:302129/1‑150 (MQ=255)                                                                                                    cTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTa                                                >  2:1053273/1‑149 (MQ=255)                                                                                                      aTCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGt                                             <  2:350623/150‑1 (MQ=255)                                                                                                        cccTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTg                                          >  2:345280/1‑151 (MQ=255)                                                                                                            gCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACtt                                      <  2:232769/151‑1 (MQ=255)                                                                                                              aTTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCg                                    <  1:399626/151‑1 (MQ=255)                                                                                                              aTTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCg                                    <  1:675244/151‑1 (MQ=255)                                                                                                               ttGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCg                                    >  1:930926/1‑150 (MQ=255)                                                                                                                tGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGa                                   <  2:291292/150‑1 (MQ=255)                                                                                                                   aaaaCCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAAtt                               >  2:288474/1‑151 (MQ=255)                                                                                                                             atgaGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTaa                     <  2:969677/151‑1 (MQ=255)                                                                                                                             atgaGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTaa                     <  2:155908/151‑1 (MQ=255)                                                                                                                               gaGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAAcc                   <  1:499181/151‑1 (MQ=255)                                                                                                                               gaGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAAcc                   >  1:356027/1‑151 (MQ=255)                                                                                                                                aGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCg                  <  1:998224/151‑1 (MQ=255)                                                                                                                                    ctTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTc                >  1:144937/1‑149 (MQ=255)                                                                                                                                     ttCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAgc             >  2:194036/1‑151 (MQ=255)                                                                                                                                       ccGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAgcg            <  1:905502/150‑1 (MQ=255)                                                                                                                                        cGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGcgcg          >  2:921926/1‑151 (MQ=255)                                                                                                                                             aGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAAATACCCTGCTAACCGTCAGCGCGGTgg      >  1:181711/1‑150 (MQ=255)                                                                                                                                              ggggCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGAt    <  1:55859/151‑1 (MQ=255)                                                                                                                                                ggCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATgg  >  1:636807/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                    TATCCGTCGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGG  >  NC_000913/1931228‑1931521

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