Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,910,662	A→T	D288E (GAT→GAA)	pstS ←	phosphate ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,910,662	0	A	T	100.0%	366.0 / NA	91	D288E (GAT→GAA)	pstS	phosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (54/37);  total (54/37)

ATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGT  >  NC_000913/3910516‑3910806                                                                                                                                                   |                                                                                                                                                  atatTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCtt                                                                                                                                               >  1:212288/1‑151 (MQ=255) atatTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCtt                                                                                                                                               >  2:227821/1‑151 (MQ=255) atatTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCtt                                                                                                                                               >  1:829845/1‑151 (MQ=255)    ttGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgt                                                                                                                                            >  2:72338/1‑151 (MQ=255)       gTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGa                                                                                                                                         >  1:347347/1‑151 (MQ=255)         cTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGAt                                                                                                                                        >  2:173631/1‑150 (MQ=255)         cTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATc                                                                                                                                       <  2:1190978/151‑1 (MQ=255)              aCGCGGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAAt                                                                                                                                  <  2:1413799/151‑1 (MQ=255)                   gcgcGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACg                                                                                                                             >  2:1596772/1‑151 (MQ=255)                   gcgcGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACg                                                                                                                             >  2:1593449/1‑151 (MQ=255)                   gcgcGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACg                                                                                                                             >  2:1316112/1‑151 (MQ=255)                   gcgcGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACg                                                                                                                             >  1:584277/1‑151 (MQ=255)                   gcgcGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACg                                                                                                                             >  2:1433639/1‑151 (MQ=255)                    cgcgAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAAc                                                                                                                              <  2:1509198/149‑1 (MQ=255)                    cgcgAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACg                                                                                                                             <  1:8562/150‑1 (MQ=255)                        aaCCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGt                                                                                                                         <  1:983444/150‑1 (MQ=255)                        aaCCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTa                                                                                                                        >  1:217378/1‑151 (MQ=255)                         aCCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTag                                                                                                                       <  2:491134/151‑1 (MQ=255)                         aCCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTag                                                                                                                       <  1:898377/151‑1 (MQ=255)                          ccTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTaga                                                                                                                      >  2:46010/1‑151 (MQ=255)                           cTGCTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTAGGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTagag                                                                                                                     >  1:566040/1‑151 (MQ=255)                            tGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAgg                                                                                                                    <  1:1315031/151‑1 (MQ=255)                            tGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAgg                                                                                                                    <  2:1356596/151‑1 (MQ=255)                            tGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAgg                                                                                                                    <  2:959043/151‑1 (MQ=255)                             gTTCAACTACACTATCTGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGt                                                                                                                   >  2:438067/1‑151 (MQ=255)                             gTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGt                                                                                                                   >  2:197963/1‑151 (MQ=255)                             gTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGt                                                                                                                   >  2:214457/1‑151 (MQ=255)                             gTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGt                                                                                                                   >  1:1724910/1‑151 (MQ=255)                              ttCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGt                                                                                                                   <  2:874520/150‑1 (MQ=255)                              ttCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGt                                                                                                                   <  2:761385/150‑1 (MQ=255)                                cAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAAt                                                                                                                >  2:1012083/1‑151 (MQ=255)                                cAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAAt                                                                                                                >  2:38411/1‑151 (MQ=255)                                cAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAAt                                                                                                                <  2:584500/151‑1 (MQ=255)                                cAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTCCAGCACTTCTGTGCCTTGTTCGGGTTTCTTCCGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAAt                                                                                                                >  2:1452734/1‑151 (MQ=255)                                   cTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAgg                                                                                                             >  2:1026376/1‑151 (MQ=255)                                   cTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAgg                                                                                                             <  2:1009274/151‑1 (MQ=255)                                   cTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAgg                                                                                                             >  2:1652845/1‑151 (MQ=255)                                    tACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGc                                                                                                            <  2:261498/151‑1 (MQ=255)                                    tACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGc                                                                                                            <  1:25415/151‑1 (MQ=255)                                    tACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGc                                                                                                            <  2:67337/151‑1 (MQ=255)                                     acacTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGc                                                                                                            >  2:1286122/1‑150 (MQ=255)                                     acacTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGc                                                                                                            >  1:1016030/1‑150 (MQ=255)                                          aTCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCa                                                                                                      >  1:438841/1‑151 (MQ=255)                                          aTCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCa                                                                                                      >  2:1220753/1‑151 (MQ=255)                                          aTCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCa                                                                                                      >  2:199426/1‑151 (MQ=255)                                          aTCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCa                                                                                                      <  2:954673/151‑1 (MQ=255)                                            ccGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATc                                                                                                    <  2:1241576/151‑1 (MQ=255)                                                   gCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGcc                                                                                              <  2:1473679/150‑1 (MQ=255)                                                      ggCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCttt                                                                                           >  2:1656955/1‑150 (MQ=255)                                                            aTCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggtt‑g                                                                                   >  2:903812/1‑151 (MQ=255)                                                            aTCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggtt‑g                                                                                   >  2:959044/1‑151 (MQ=255)                                                            aTCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggtt‑g                                                                                   >  1:565659/1‑151 (MQ=255)                                                                  gTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATCCATCTTCGCCTTTCTGGTT‑GGTc                                                                                <  2:135266/148‑1 (MQ=255)                                                                          ccTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAg                                                                      <  2:817454/150‑1 (MQ=255)                                                                          ccTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAg                                                                      <  2:181189/150‑1 (MQ=255)                                                                             gTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCgca                                                                  >  1:1446293/1‑149 (MQ=255)                                                                                        gTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGGTCAGATCCTGAGCGAAGGTTTTGCt                                                         >  1:1557453/1‑150 (MQ=255)                                                                                        gTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTc                                                        >  1:1516564/1‑150 (MQ=255)                                                                                        gTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTc                                                        <  2:575212/150‑1 (MQ=255)                                                                                        gTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTc                                                        >  2:1043302/1‑150 (MQ=255)                                                                                           ttGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGt                                                    >  1:1103865/1‑151 (MQ=255)                                                                                            tGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGt                                                    >  1:726647/1‑150 (MQ=255)                                                                                            tGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTc                                                   <  2:646571/151‑1 (MQ=255)                                                                                            tGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTc                                                   <  1:1194250/151‑1 (MQ=255)                                                                                            tGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTc                                                   <  1:45994/151‑1 (MQ=255)                                                                                            tGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTc                                                   <  1:1006669/151‑1 (MQ=255)                                                                                                  cccAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCAcc                                             <  1:19765/151‑1 (MQ=255)                                                                                                          aagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTtgctg                                     <  2:466191/151‑1 (MQ=255)                                                                                                          aagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTtgctg                                     <  2:462924/151‑1 (MQ=255)                                                                                                            gaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCa                                   >  1:1009499/1‑151 (MQ=255)                                                                                                            gaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCa                                   >  2:624841/1‑151 (MQ=255)                                                                                                               tttCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTa                                >  2:1353048/1‑151 (MQ=255)                                                                                                                     cACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCg                             <  1:1381811/148‑1 (MQ=255)                                                                                                                     cACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCg                             <  2:492346/148‑1 (MQ=255)                                                                                                                       cTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGtt                        >  2:1442874/1‑151 (MQ=255)                                                                                                                       cTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGtt                        >  2:917763/1‑151 (MQ=255)                                                                                                                       cTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGtt                        >  1:652169/1‑151 (MQ=255)                                                                                                                       cTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGtt                        >  1:1500287/1‑151 (MQ=255)                                                                                                                         tCTGTGCCTTGTTTTGGATTCTTCTGTTCTTTGTGGATTAGAATGAACGTGGTAGCGGCCATAGGCCATGCATCTTCGCCTTTCTCGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGACGgttt                      >  1:860772/1‑151 (MQ=255)                                                                                                                             tgCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGGGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCACCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTttcttc                  <  2:633075/151‑1 (MQ=255)                                                                                                                              gCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGAGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTtcg                 >  2:468017/1‑151 (MQ=255)                                                                                                                                 ttGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTtcggt               >  1:1532549/1‑150 (MQ=255)                                                                                                                                   gTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCATCTTTTGCTGCATTAGCGAAGCTTGCTtcggtcgg            >  1:1121105/1‑151 (MQ=255)                                                                                                                                   gTTCCGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTtcggtcgg            >  1:864915/1‑151 (MQ=255)                                                                                                                                      cTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACt         <  1:131737/151‑1 (MQ=255)                                                                                                                                        ggTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTaa       >  1:498894/1‑151 (MQ=255)                                                                                                                                         gTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTaa       >  1:1184614/1‑150 (MQ=255)                                                                                                                                          tttCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAAcc     <  2:212352/151‑1 (MQ=255)                                                                                                                                            tcttctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCgg   <  2:739579/151‑1 (MQ=255)                                                                                                                                             cttctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGt  >  1:513143/1‑151 (MQ=255)                                                                                                                                             cttctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGt  >  2:1058144/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                  ATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTT‑GGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGT  >  NC_000913/3910516‑3910806

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