Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,398,201	G→A	P141L (CCG→CTG)	uspE ←	stress‑induced protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,398,201	0	G	A	100.0%	358.0 / NA	89	P141L (CCG→CTG)	uspE	stress‑induced protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (41/48);  total (41/48)

GTATGGTTGACTTGTTCTGCCAGTTCGATCGTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTCCGCCGCTGATCACTTCCTGAATGATGGCTTCGAAAGGACGGTTATG  >  NC_000913/1398054‑1398346                                                                                                                                                    |                                                                                                                                                  gTATGGTTGACTTGTTCTGCCAGTTCGATCGTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGc                                                                                                                                                <  2:143454/151‑1 (MQ=255)  tATGGTTGACTTGTTCTGCCAGTTCGATCGTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCt                                                                                                                                               >  2:538732/1‑151 (MQ=255)   aTGGTTGACTTGTTCTGCCAGTTCGATCGTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCtt                                                                                                                                              >  2:1056244/1‑151 (MQ=255)    tGGTTGACTTGTTCTGCCAGTTCGATCGTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTg                                                                                                                                             >  1:1038653/1‑151 (MQ=255)       ttGACTTGTTCTGCCAGTTCGATCGTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCATGCTTggg                                                                                                                                           <  2:1111050/150‑1 (MQ=255)         gACTTGTTCTGCCAGTTCGATCGTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGc                                                                                                                                          <  2:1161212/149‑1 (MQ=255)         gACTTGTTCTGCCAGTTCGATCGTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCAt                                                                                                                                        <  1:622400/151‑1 (MQ=255)                 cTGCCAGTTCGATCGTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGt                                                                                                                                  >  2:561883/1‑149 (MQ=255)                  tGCCAGTTCGATCGTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTg                                                                                                                                             >  2:1262135/1‑137 (MQ=255)                  tGCCAGTTCGATCGTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTg                                                                                                                                             <  1:1262236/137‑1 (MQ=255)                  tGCCAGTTCGATCGTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAAc                                                                                                                               >  2:959922/1‑151 (MQ=255)                    ccAGTTCGATCTTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCTCTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAg                                                                                                                             <  2:770064/151‑1 (MQ=255)                         tCGATCGTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGcc                                                                                                                        >  2:587188/1‑151 (MQ=255)                          cGATCGTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGcaa                                                                                                                       <  1:253508/151‑3 (MQ=255)                          cGATCGTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCa                                                                                                                       <  1:1328189/151‑1 (MQ=255)                            aTCGTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGt                                                                                                                     >  1:1473168/1‑151 (MQ=255)                            aTCGTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGt                                                                                                                     <  2:818178/151‑1 (MQ=255)                              cGTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTcc                                                                                                                   >  2:1428166/1‑151 (MQ=255)                              cGTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTc                                                                                                                    <  2:1383628/150‑1 (MQ=255)                               gTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCg                                                                                                                  <  1:1343082/151‑1 (MQ=255)                               gTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCg                                                                                                                  >  1:1025970/1‑151 (MQ=255)                                 ctctTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGtt                                                                                                                >  2:654074/1‑151 (MQ=255)                                     ttGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTgg                                                                                                              <  2:72819/149‑1 (MQ=255)                                     ttGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTg                                                                                                               <  1:602984/148‑1 (MQ=255)                                     ttGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCg                                                                                                            >  2:1084968/1‑151 (MQ=255)                                        aCCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTaa                                                                                                         <  2:599563/151‑1 (MQ=255)                                          cAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTaaaa                                                                                                       <  2:1151800/151‑1 (MQ=255)                                            gTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAAt                                                                                                     >  2:108/1‑151 (MQ=255)                                              ttttCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATc                                                                                                    <  2:603649/150‑1 (MQ=255)                                               tttCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCAc                                                                                                  >  2:1110911/1‑151 (MQ=255)                                                 tCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCAcc                                                                                                 <  1:1130818/150‑1 (MQ=255)                                                  cATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCg                                                                                                <  1:186752/150‑1 (MQ=255)                                                  cATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGc                                                                                               >  2:791286/1‑151 (MQ=255)                                                  cATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGc                                                                                               <  1:818247/151‑1 (MQ=255)                                                     tGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTc                                                                                            >  1:1369055/1‑151 (MQ=255)                                                        gcgcATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCa                                                                                          <  2:1028555/150‑1 (MQ=255)                                                            aTTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACgg                                                                                     >  2:1292740/1‑151 (MQ=255)                                                              tATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTc                                                                                   >  1:1220975/1‑151 (MQ=255)                                                               aTGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTc                                                                                   <  1:774936/150‑1 (MQ=255)                                                               aTGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTc                                                                                   <  2:876092/150‑1 (MQ=255)                                                               aTGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTc                                                                                   >  1:650456/1‑150 (MQ=255)                                                                    aCGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGtt                                                                              <  1:666950/150‑1 (MQ=255)                                                                    aCGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTg                                                                             <  1:210374/151‑1 (MQ=255)                                                                     cGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCCGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTgg                                                                            >  2:767074/1‑151 (MQ=255)                                                                           ttCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGcc                                                                      >  2:754446/1‑151 (MQ=255)                                                                            tCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGcc                                                                      >  2:732801/1‑150 (MQ=255)                                                                             cACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCAt                                                                    >  1:1456091/1‑151 (MQ=255)                                                                                tGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCAtttt                                                                 >  2:962169/1‑151 (MQ=255)                                                                                tGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCAtttt                                                                 >  1:195804/1‑151 (MQ=255)                                                                                 ggCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCAttttt                                                                <  1:947065/151‑1 (MQ=255)                                                                                 ggCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCAttt                                                                  <  1:601961/149‑1 (MQ=255)                                                                                 ggCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCAttt                                                                  <  1:303444/149‑1 (MQ=255)                                                                                  gCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTa                                                               <  2:264878/151‑1 (MQ=255)                                                                                    gagaTTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAg                                                              <  1:1191291/150‑1 (MQ=255)                                                                                    gagaTTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAg                                                              <  1:1192085/150‑1 (MQ=255)                                                                                      gaTTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCAc                                                           <  1:237453/151‑1 (MQ=255)                                                                                        ttCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCAcc                                                          >  1:359896/1‑150 (MQ=255)                                                                                        ttCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCAcc                                                          >  2:738824/1‑150 (MQ=255)                                                                                         tCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCa                                                         >  2:639892/1‑150 (MQ=255)                                                                                          cACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCGCCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCaaa                                                       <  1:1392611/151‑1 (MQ=255)                                                                                           aCCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAAt                                                      >  2:610773/1‑151 (MQ=255)                                                                                            ccGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATc                                                     <  1:1056340/151‑1 (MQ=255)                                                                                            ccGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATc                                                     <  1:561948/151‑1 (MQ=255)                                                                                             cGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCg                                                    >  2:810623/1‑151 (MQ=255)                                                                                              gCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCgt                                                   >  1:992401/1‑151 (MQ=255)                                                                                              gCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCgt                                                   >  1:476972/1‑151 (MQ=255)                                                                                                 accaGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTc                                                <  1:1057075/151‑1 (MQ=255)                                                                                                  ccaGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTcc                                               <  1:1137285/151‑1 (MQ=255)                                                                                                   caGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTccg                                              <  1:1115250/151‑1 (MQ=255)                                                                                                      cgcCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTccgccg                                           <  1:829719/151‑1 (MQ=255)                                                                                                         cTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTCCGCCGCTg                                        >  1:659861/1‑151 (MQ=255)                                                                                                                  ttCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTCCGCCGCTGATCACTTcc                               >  1:3188/1‑151 (MQ=255)                                                                                                                   tCCGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTCCGCCGCTGATCACTTc                                <  1:587257/149‑1 (MQ=255)                                                                                                                     cGGCCACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTCCGCCGCTGATCACTTCCTg                             >  1:1280474/1‑150 (MQ=255)                                                                                                                        ccACGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTCCGCCGCTGATCACTTCCTGAatg                         <  1:538794/151‑1 (MQ=255)                                                                                                                          aCGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTCCGCCGCTGATCACTTCCTGAatgat                       <  2:1328752/151‑1 (MQ=255)                                                                                                                           cGGCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTCCGCCGCTGATCACTTCCTGAatgatg                      <  2:705569/151‑1 (MQ=255)                                                                                                                            ggCTGGTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTCCGCCGCTGATCACTTCCTGAATGATgg                     >  1:1388247/1‑151 (MQ=255)                                                                                                                                ggTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTCCGCCGCTGATCACTTCCTGAATGATGGCt                   <  1:140824/149‑1 (MQ=255)                                                                                                                                ggTCTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTCCGCCGCTGATCACTTCCTGAATGATGGCt                   <  1:1229146/149‑1 (MQ=255)                                                                                                                                   cTTTCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTCCGCCGCTGATCACTTCCTGAATGATGGCTTCGa               >  1:1354911/1‑150 (MQ=255)                                                                                                                                     ttCACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTCCGCCGCTGATCACTTCCTGAATGATGGCTTCGaa              >  2:631878/1‑149 (MQ=255)                                                                                                                                       cACCATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTCCGCCGCTGATCACTTCCTGAATGATGGCTTCGAAAgt           <  2:702529/150‑2 (MQ=255)                                                                                                                                         ccATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTCCGCCGCTGATCACTTCCTGAATGATGGCTTCGAAAGGACg        <  1:1085068/151‑1 (MQ=255)                                                                                                                                          cATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTCCGCCGCTGATCACTTCCTGAATGATGGCTTCGAAAGGACgg       >  1:176889/1‑151 (MQ=255)                                                                                                                                          cATCCACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTCCGCCGCTGATCACTTCCTGAATGATGGCTTCGAAAGGACgg       >  1:1381203/1‑151 (MQ=255)                                                                                                                                             ccACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTCCGCCGCTGATCACTTCCTGAATGATGGCTTCGAAAGGACGGtt     >  2:1169327/1‑150 (MQ=255)                                                                                                                                             ccACACCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTCCGCCGCTGATCACTTCCTGAATGATGGCTTCGAAAGGACGGTTa    <  1:654149/151‑1 (MQ=255)                                                                                                                                                cacCAGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTCCGCCGCTGATCACTTCCTGAATGATGGCTTCGAAAGGACGGTTATg  <  2:465214/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GTATGGTTGACTTGTTCTGCCAGTTCGATCGTCTCTTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGTTGGTGCGCCATTTTTAGCACCAAATCGTGTCCGCCGCTGATCACTTCCTGAATGATGGCTTCGAAAGGACGGTTATG  >  NC_000913/1398054‑1398346

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