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%	243.2 / NA	69	Q165R (CAG→CGG)	purT	phosphoribosylglycinamide formyltransferase 2
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/39);  total (30/39)

GTATCCGTCGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGG  >  NC_000913/1931227‑1931517                                                                                                                                                    |                                                                                                                                                gTATCCGTCGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGAc                                                                                                                                              <  2:49261/151‑1 (MQ=255)   aTCCGTCGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGAc                                                                                                                                              >  2:488053/1‑149 (MQ=255)       gTCGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTAt                                                                                                                                        >  2:583138/1‑151 (MQ=255)       gTCGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTAt                                                                                                                                        <  1:523839/151‑1 (MQ=255)           ccTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGt                                                                                                                                    >  1:5844/1‑151 (MQ=255)            cTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGtt                                                                                                                                   >  2:718557/1‑151 (MQ=255)              ggcggcAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTc                                                                                                                                  <  1:515885/150‑1 (MQ=255)                     gaagaGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGc                                                                                                                          <  1:414974/151‑1 (MQ=255)                          gctgcagctgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACtt                                                                                                                     <  2:754071/151‑1 (MQ=255)                          gctgcagctgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACt                                                                                                                      <  2:740589/150‑1 (MQ=255)                              cagctgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTc                                                                                                                 >  1:588309/1‑151 (MQ=255)                                gctgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAg                                                                                                               <  1:318206/151‑1 (MQ=255)                                 ctgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAgg                                                                                                              >  2:721119/1‑151 (MQ=255)                                   gcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGc                                                                                                             <  1:52308/150‑1 (MQ=255)                                    cccACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGCGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAgg                                                                                                              <  1:17686/148‑1 (MQ=255)                                           ccacttATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGt                                                                                                     <  1:752219/150‑1 (MQ=255)                                                 aTCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTc                                                                                               >  1:429249/1‑150 (MQ=255)                                                     ttttGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCa                                                                                           >  1:13408/1‑150 (MQ=255)                                                     ttttGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCa                                                                                           >  1:337225/1‑150 (MQ=255)                                                     ttttGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCa                                                                                           >  1:2549/1‑150 (MQ=255)                                                         gCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTGCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCg                                                                                      <  1:423399/151‑1 (MQ=255)                                                            gATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTc                                                                                   <  2:269710/151‑1 (MQ=255)                                                             aTAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTcg                                                                                  <  1:686778/151‑1 (MQ=255)                                                                 cGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGcc                                                                              <  1:175633/151‑1 (MQ=255)                                                                      gCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGc                                                                         <  1:593224/151‑1 (MQ=255)                                                                       ccTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCg                                                                        <  2:561024/151‑1 (MQ=255)                                                                         ttttCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCggg                                                                      >  1:581996/1‑151 (MQ=255)                                                                         ttttCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCggg                                                                      >  2:329938/1‑151 (MQ=255)                                                                            tCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGcc                                                                    >  2:348290/1‑150 (MQ=255)                                                                               gcgAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGt                                                                >  1:112532/1‑151 (MQ=255)                                                                               gcgAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGt                                                                <  2:654129/151‑1 (MQ=255)                                                                                 gAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTa                                                               >  1:7506/1‑150 (MQ=255)                                                                                       gTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCGGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAAttgtt                                                         <  1:295229/150‑1 (MQ=255)                                                                                       gTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAAttgttg                                                        <  1:454575/151‑1 (MQ=255)                                                                                        ttGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGAACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAAttgttg                                                        >  1:314691/1‑150 (MQ=255)                                                                                         tGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGa                                                       >  2:628536/1‑150 (MQ=255)                                                                                          gCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGa                                                       >  1:173352/1‑149 (MQ=255)                                                                                          gCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAg                                                     <  2:296721/151‑1 (MQ=255)                                                                                            tGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAg                                                     >  2:195088/1‑149 (MQ=255)                                                                                            tGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGc                                                   >  2:476403/1‑151 (MQ=255)                                                                                              aCATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGATCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGcg                                                  <  2:304936/150‑1 (MQ=255)                                                                                              aCATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGcgt                                                 <  2:169417/151‑1 (MQ=255)                                                                                                aTTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGcgtcg                                               >  1:589184/1‑151 (MQ=255)                                                                                                   ggCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGtt                                             >  1:268403/1‑150 (MQ=255)                                                                                                         cccTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCggg                                                                      <  2:2790/119‑1 (MQ=255)                                                                                                         cccTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCggg                                                                      >  1:2790/1‑119 (MQ=255)                                                                                                          ccTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGa                                     <  1:199204/151‑1 (MQ=255)                                                                                                             gCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACt                                   <  2:192252/150‑1 (MQ=255)                                                                                                             gCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACt                                   <  2:376213/150‑1 (MQ=255)                                                                                                              cATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTc                                 <  1:660429/151‑1 (MQ=255)                                                                                                              cATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTc                                 <  2:113216/151‑1 (MQ=255)                                                                                                                  gTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGaaa                             <  1:718601/151‑1 (MQ=255)                                                                                                                  gTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACCTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGGCGCGTAATTGTTGAAGGCGTCGTTAATTTTGACTGCGaa                              <  2:451882/150‑1 (MQ=255)                                                                                                                    aaaaCCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAAtt                           <  2:504077/151‑1 (MQ=255)                                                                                                                      aaCCGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTAc                         >  2:391059/1‑151 (MQ=255)                                                                                                                         cGGTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCt                      <  2:337211/151‑1 (MQ=255)                                                                                                                          ggTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTg                     <  1:268815/151‑1 (MQ=255)                                                                                                                           gTGATGAGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGc                    <  2:149524/151‑1 (MQ=255)                                                                                                                            tgatgaGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCt                   >  2:673867/1‑151 (MQ=255)                                                                                                                               tgaGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAAc                >  1:736551/1‑151 (MQ=255)                                                                                                                                gaGCTCTTCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAAcc               <  2:308408/151‑1 (MQ=255)                                                                                                                                      ttCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAgc         >  1:477313/1‑151 (MQ=255)                                                                                                                                       tCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCa           <  1:517897/148‑1 (MQ=255)                                                                                                                                       tCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCa           <  1:721163/148‑1 (MQ=255)                                                                                                                                       tCCGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAgcg        >  1:142348/1‑151 (MQ=255)                                                                                                                                        ccGGCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAgcgc       <  1:434839/151‑1 (MQ=255)                                                                                                                                           gCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCgg     >  1:613871/1‑150 (MQ=255)                                                                                                                                           gCAAGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGt    <  2:272486/151‑1 (MQ=255)                                                                                                                                             aaGGGGCGGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTgg  >  2:352386/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                GTATCCGTCGCCTGGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGG  >  NC_000913/1931227‑1931517

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: ‑