Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,908,340	A→G	S9S (AGT→AGC)	pstB ←	phosphate ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,908,340	0	A	G	100.0%	323.0 / NA	90	S9S (AGT→AGC)	pstB	phosphate ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (38/52);  total (38/52)

TACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGT  >  NC_000913/3908204‑3908473                                                                                                                                         |                                                                                                                                      tACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTc                                                                                                                                    <  1:1209184/140‑1 (MQ=255)  aCGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCg                                                                                                                                   >  1:1024471/1‑140 (MQ=255)   cGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCgg                                                                                                                                  >  1:160264/1‑140 (MQ=255)   cGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCgg                                                                                                                                  >  2:908947/1‑140 (MQ=255)    gcagcaGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCggg                                                                                                                                 >  2:347739/1‑140 (MQ=255)     cagcagCGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCgggg                                                                                                                                <  2:1374839/140‑1 (MQ=255)     cagcagCGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCgggg                                                                                                                                <  2:160264/140‑1 (MQ=255)     cagcagCGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCgggg                                                                                                                                <  2:300828/140‑1 (MQ=255)       gcagcGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCggg                                                                                                                                 >  1:457765/1‑137 (MQ=255)       gcagcGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCggg                                                                                                                                 <  2:457765/137‑1 (MQ=255)        cagcGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAg                                                                                                                             >  1:1068184/1‑140 (MQ=255)         agcGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGt                                                                                                                            <  1:1209628/140‑1 (MQ=255)         agcGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGt                                                                                                                            <  2:1175670/140‑1 (MQ=255)         agcGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGt                                                                                                                            >  2:994697/1‑140 (MQ=255)          gcGACGATTTTCCGCAGGCGGCCGGCCGGATAAACGCCGTTGCCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGACCCTGAATTTTGCTCGGGGCAGtt                                                                                                                           <  2:134979/140‑1 (MQ=255)           cGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGttt                                                                                                                          >  2:1462127/1‑140 (MQ=255)           cGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGttt                                                                                                                          <  1:483889/140‑1 (MQ=255)           cGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGttt                                                                                                                          <  1:802608/140‑1 (MQ=255)           cGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGttt                                                                                                                          <  2:1024471/140‑1 (MQ=255)           cGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGttt                                                                                                                          >  2:1583844/1‑140 (MQ=255)           cGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGttt                                                                                                                          >  1:1447554/1‑140 (MQ=255)           cGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGttt                                                                                                                          <  2:1608150/140‑1 (MQ=255)           cGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGttt                                                                                                                          <  2:821523/140‑1 (MQ=255)           cGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGttt                                                                                                                          >  2:802608/1‑140 (MQ=255)              cGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCaa                                                                                                                       >  1:1051878/1‑140 (MQ=255)                aTTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAAcc                                                                                                                     >  1:675749/1‑140 (MQ=255)                aTTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAAcc                                                                                                                     >  1:1028487/1‑140 (MQ=255)                          gCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTg                                                                                                           <  2:1447554/140‑1 (MQ=255)                          gCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTg                                                                                                           <  1:1462127/140‑1 (MQ=255)                          gCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTg                                                                                                           <  1:1583844/140‑1 (MQ=255)                            cGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCa                                                                                                         <  1:347739/140‑1 (MQ=255)                            cGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCa                                                                                                         <  2:1051878/140‑1 (MQ=255)                                   cccGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCtttt                                                                                                  >  2:1640933/1‑140 (MQ=255)                                    ccGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTc                                                                                                 <  2:513013/140‑1 (MQ=255)                                    ccGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTc                                                                                                 <  1:1320815/140‑1 (MQ=255)                                     cGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTcc                                                                                                <  2:1028487/140‑1 (MQ=255)                                         aaaCGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCAt                                                                                            <  1:824680/140‑1 (MQ=255)                                         aaaCGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCAt                                                                                            <  1:908947/140‑1 (MQ=255)                                          aaCGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATc                                                                                           <  1:607694/140‑1 (MQ=255)                                          aaCGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATc                                                                                           <  1:994697/140‑1 (MQ=255)                                           aCGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCg                                                                                          <  2:718269/140‑1 (MQ=255)                                             gCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTg                                                                                        >  2:1119015/1‑140 (MQ=255)                                               cGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGcc                                                                                      >  1:826970/1‑140 (MQ=255)                                                 ttACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTATCCTCATCGTGCCTg                                                                                    >  2:473983/1‑140 (MQ=255)                                                  tACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGa                                                                                   <  1:1119015/140‑1 (MQ=255)                                                    ccTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATg                                                                                 <  2:546972/140‑1 (MQ=255)                                                          ttttAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTc                                                                           <  1:1641935/140‑1 (MQ=255)                                                                cGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAAt                                                                     <  2:534836/140‑1 (MQ=255)                                                                cGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTATCCTCATCGTGCCTGATGCACTTCGCTAAt                                                                     <  1:473983/140‑1 (MQ=255)                                                                     tCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGc                                                                <  2:675749/140‑1 (MQ=255)                                                                     tCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGc                                                                <  1:1640933/140‑1 (MQ=255)                                                                      ccAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCt                                                               <  1:969396/140‑1 (MQ=255)                                                                      ccAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCt                                                               <  2:1629600/140‑1 (MQ=255)                                                                      ccAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCAAACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCt                                                               <  1:1327920/140‑1 (MQ=255)                                                                          gTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTAcc                                                           >  2:36205/1‑140 (MQ=255)                                                                           ttGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCg                                                          <  2:826970/140‑1 (MQ=255)                                                                             gATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTa                                                        >  1:41071/1‑140 (MQ=255)                                                                              aTGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAAt                                                                     <  1:286267/126‑1 (MQ=255)                                                                              aTGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAAt                                                                     >  2:286267/1‑126 (MQ=255)                                                                              aTGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAg                                                       >  2:1184181/1‑140 (MQ=255)                                                                               tGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAgg                                                      >  1:758724/1‑140 (MQ=255)                                                                               tGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAgg                                                      >  2:617875/1‑140 (MQ=255)                                                                               tGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTGCGCTAATCAGGCTTACCGTAgg                                                      >  1:1177250/1‑140 (MQ=255)                                                                                gTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGc                                                     <  1:36205/140‑1 (MQ=255)                                                                                    ttCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGa                                                 >  2:1201595/1‑140 (MQ=255)                                                                                     tCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGAt                                                <  2:1177250/140‑1 (MQ=255)                                                                                         ggCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAgg                                            >  1:276556/1‑140 (MQ=255)                                                                                         ggCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAgg                                            >  1:363457/1‑140 (MQ=255)                                                                                             tGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGtt                                        <  1:1184181/140‑1 (MQ=255)                                                                                                 aTTTGCCGTAGTAGAAGTTCAAATTACGAGCCTGAATTTTGCTCGGGGCAGTTTCAACCAAGCTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTACTCAGGCTTACCGTAGGCCGGACCAGGCTTTCAcg                                    <  2:1068184/140‑1 (MQ=255)                                                                                                         tagtagAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATcc                            >  1:1076786/1‑140 (MQ=255)                                                                                                            tagAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGc                         <  2:363457/140‑1 (MQ=255)                                                                                                              gAAGTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCaa                       <  1:617875/140‑1 (MQ=255)                                                                                                                 gTTCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAgt                    >  2:93339/1‑140 (MQ=255)                                                                                                                  ttCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAgtg                   >  1:851022/1‑140 (MQ=255)                                                                                                                   tCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAgtgt                  <  1:1201595/140‑1 (MQ=255)                                                                                                                   tCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAgtgt                  <  2:1076786/140‑1 (MQ=255)                                                                                                                   tCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAgtgt                  <  1:1364061/140‑1 (MQ=255)                                                                                                                   tCAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAgtgt                  >  2:1284531/1‑140 (MQ=255)                                                                                                                    cAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGtt                 >  2:617958/1‑140 (MQ=255)                                                                                                                    cAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGtt                 >  1:638687/1‑140 (MQ=255)                                                                                                                    cAAATTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGtt                 >  2:77838/1‑140 (MQ=255)                                                                                                                     aaaTTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTa                <  1:617958/140‑1 (MQ=255)                                                                                                                     aaaCTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTa                <  1:95970/140‑1 (MQ=255)                                                                                                                      aaTTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAgtgt                  >  2:643523/1‑137 (MQ=255)                                                                                                                      aaTTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAgtgt                  <  1:643523/137‑1 (MQ=255)                                                                                                                      aaTTACGAACCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTaa               >  1:874311/1‑140 (MQ=255)                                                                                                                             aaCCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTa        <  1:77838/140‑1 (MQ=255)                                                                                                                             aaCCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTa        <  2:601489/140‑1 (MQ=255)                                                                                                                              aCCTGAATTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTAc       <  2:638687/140‑1 (MQ=255)                                                                                                                                   aaTTTTGCTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGt  >  2:1370988/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                      TACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGT  >  NC_000913/3908204‑3908473

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