Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,137,229	A→G	S143S (AGT→AGC)	yijF ←	DUF1287 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,137,229	0	A	G	100.0%	294.2 / NA	82	S143S (AGT→AGC)	yijF	DUF1287 family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (38/44);  total (38/44)

TTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTACTCAGCAAAATTC  >  NC_000913/4137097‑4137361                                                                                                                                     |                                                                                                                                     ttcCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTgggg                                                                                                                                <  2:368075/139‑1 (MQ=255)   ccTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGtt                                                                                                                              <  1:139302/139‑1 (MQ=255)     tGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCt                                                                                                                            <  2:33371/139‑1 (MQ=255)      gcgcgcC‑GCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTg                                                                                                                          >  2:1408952/1‑139 (MQ=255)       cgcgcCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCt                                                                                                                                    >  1:1152472/1‑129 (MQ=255)       cgcgcCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCt                                                                                                                                    <  2:1152472/129‑1 (MQ=255)       cgcgcCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGt                                                                                                                               <  1:328254/134‑1 (MQ=255)       cgcgcCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGt                                                                                                                               >  2:328254/1‑134 (MQ=255)       cgcgcCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTg                                                                                                                          >  2:1436846/1‑139 (MQ=255)       cgcgcCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTg                                                                                                                          >  2:1163290/1‑139 (MQ=255)         cgcCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCt                                                                                                                        >  2:1040204/1‑139 (MQ=255)            cTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCtggt                                                                                                                     <  1:1436846/139‑1 (MQ=255)               cgcCAATGTTATGTGTCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATCTGCCAGCCCGTTAGCCAGTCGCCAGGAGACGATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCtggtggt                                                                                                                  <  2:1262881/139‑2 (MQ=255)                gcCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTgggg                                                                                                                 <  1:1580981/139‑1 (MQ=255)                 ccAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                <  1:535751/139‑1 (MQ=255)                 ccAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                <  1:453063/139‑1 (MQ=255)                  cAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCg                                                                                                               <  2:204621/139‑1 (MQ=255)                         aTGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCt                                                                                                                                    <  1:53467/111‑1 (MQ=255)                         aTGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTa                                                                                                        >  2:751744/1‑139 (MQ=255)                               cACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATgg                                                                                                  >  2:1233005/1‑139 (MQ=255)                               cACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATgg                                                                                                  >  1:596288/1‑139 (MQ=255)                                 ccAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCgtg                                                                                                             >  1:1396466/1‑126 (MQ=255)                                 ccAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCgtg                                                                                                             <  2:1396466/126‑1 (MQ=255)                                   aGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggcggc                                                                                              <  1:1040204/139‑1 (MQ=255)                                    gcggcgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTg                                                                                                                                   <  2:900480/101‑1 (MQ=255)                                    gcggcgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTg                                                                                                                                   >  1:900480/1‑101 (MQ=255)                                    gcggcgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCt                                                                                             >  1:1017565/1‑139 (MQ=255)                                    gcggcgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCt                                                                                             >  1:1118108/1‑139 (MQ=255)                                        cgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggc                                                                                                 <  1:1381053/131‑1 (MQ=255)                                        cgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAc                                                                                         <  2:596288/139‑1 (MQ=255)                                        cgTACCGTCACGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTGTTATCATggc                                                                                                 >  2:1381053/1‑131 (MQ=255)                                            ccGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGtt                                                                                                                              >  1:1977850/1‑98 (MQ=255)                                            ccGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGtt                                                                                                                              <  2:1977850/98‑1 (MQ=255)                                                  cgGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTTGCTGGGGTTCTTGCTGGTggg                                                                                                                  <  2:933188/104‑1 (MQ=255)                                                  cgGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTTGCTGGGGTTCTTGCTGGTggg                                                                                                                  >  1:933188/1‑104 (MQ=255)                                                        aaGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAg                                                                         <  2:1017565/139‑1 (MQ=255)                                                        aaGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGACTAAACCAGGTTTCCAGATTAg                                                                         <  1:1408952/139‑1 (MQ=255)                                                         aGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAgg                                                                        <  2:1409224/139‑1 (MQ=255)                                                             aTCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCAcc                                                                    <  1:751744/139‑1 (MQ=255)                                                               cTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCCTGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCg                                                                  >  2:249524/1‑139 (MQ=255)                                                                   aaCCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGATAAACCAGGTTTCCAGATTAGGCACCCGACgg                                                              >  1:1120181/1‑139 (MQ=255)                                                                     ccaccCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTg                                                            <  1:1163290/139‑1 (MQ=255)                                                                     ccaccCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTg                                                            >  2:1934565/1‑139 (MQ=255)                                                                        ccccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTg                                                                                                                                   >  1:778839/1‑65 (MQ=255)                                                                        ccccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTg                                                                                                                                   <  2:778839/65‑1 (MQ=255)                                                                        ccccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTgggg                                                                                                                 <  2:1939598/83‑1 (MQ=255)                                                                        ccccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTgggg                                                                                                                 >  1:1939598/1‑83 (MQ=255)                                                                         cccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCg                                                        >  2:1840250/1‑139 (MQ=255)                                                                         cccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCg                                                        >  2:1494970/1‑139 (MQ=255)                                                                            aaTATGTGCCAGCCCGTTATCTAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATg                                                     >  2:163090/1‑139 (MQ=255)                                                                            aaTATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATg                                                     >  2:491613/1‑139 (MQ=255)                                                                               atGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTg                                                  <  1:186288/139‑1 (MQ=255)                                                                                      aGCCCGTTATCTAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCg                                           <  1:163090/139‑1 (MQ=255)                                                                                        cccGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTgggg                                                                                                                 <  1:1632613/67‑1 (MQ=255)                                                                                        cccGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTgggg                                                                                                                 >  2:1632613/1‑67 (MQ=255)                                                                                        cccGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATg                                                                                                   >  2:1534490/1‑81 (MQ=255)                                                                                        cccGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATg                                                                                                   <  1:1534490/81‑1 (MQ=255)                                                                                             tATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCAACGGTg                                                            <  1:423369/115‑1 (MQ=255)                                                                                             tATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCAACGGTg                                                            >  2:423369/1‑115 (MQ=255)                                                                                              aTCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCttt                                   <  2:1118108/139‑1 (MQ=255)                                                                                                  aGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTgggg                                                                                                                 <  2:581692/57‑1 (MQ=255)                                                                                                  aGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTgggg                                                                                                                 >  1:581692/1‑57 (MQ=255)                                                                                                    tCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGgc                                        >  1:411616/1‑128 (MQ=255)                                                                                                    tCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGgc                                        <  2:411616/128‑1 (MQ=255)                                                                                                         aGGAGACAATATCGCCCGCTTGATTGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCAttt                        >  1:1393866/1‑139 (MQ=255)                                                                                                         aGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCggg                                         <  2:511376/122‑1 (MQ=255)                                                                                                         aGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCggg                                         >  1:511376/1‑122 (MQ=255)                                                                                                          ggAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGATAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCAtttt                       <  2:1120181/139‑1 (MQ=255)                                                                                                           gagaCAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATg                                                                                                   >  1:997366/1‑62 (MQ=255)                                                                                                           gagaCAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATg                                                                                                   <  2:997366/62‑1 (MQ=255)                                                                                                                aaTATCGCCCGCTTGATAGTGGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATg                                                                                                   >  2:1396233/1‑57 (MQ=39)                                                                                                                aaTATCGCCCGCTTGATAGTCGCTGGGGTTCTTGTTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGt                                                    <  1:747135/104‑1 (MQ=255)                                                                                                                aaTATCGCCCGCTTGATAGTCGCTGGGGTTCTTGTTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGt                                                    >  2:747135/1‑104 (MQ=255)                                                                                                                aaTATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATg                                                                                                   <  1:1396233/57‑1 (MQ=255)                                                                                                                  tatCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTa               <  1:1934565/139‑1 (MQ=255)                                                                                                                    tCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTACt             <  1:1494970/139‑1 (MQ=255)                                                                                                                     cGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTcc                                                                                <  1:1327365/71‑1 (MQ=255)                                                                                                                     cGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTcc                                                                                >  2:1327365/1‑71 (MQ=255)                                                                                                                       cccGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGa                                                                 >  2:1535540/1‑84 (MQ=255)                                                                                                                       cccGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGa                                                                 <  1:1535540/84‑1 (MQ=255)                                                                                                                       cccGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTACTCAg          >  2:890267/1‑139 (MQ=255)                                                                                                                               aTAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTACTCAGCAAAAttc  <  1:1840250/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                     TTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTACTCAGCAAAATTC  >  NC_000913/4137097‑4137361

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