Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,565,758	C→T	*461* (TAG→TAA)	dosC ←	diguanylate cyclase, cold‑ and stationary phase‑induced oxygen‑dependent biofilm regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,565,758	0	C	T	100.0%	281.1 / NA	87	*461* (TAG→TAA)	dosC	diguanylate cyclase, cold‑ and stationary phase‑induced oxygen‑dependent biofilm regulator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (40/47);  total (40/47)

CTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAG  >  NC_000913/1565634‑1565893                                                                                                                              |                                                                                                                                        cTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCttt                                                                                                                           >  1:73613/1‑139 (MQ=255)   tcatcaTTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTATCCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTc                                                                                                                          <  1:1123945/139‑1 (MQ=255)      tcatTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCaga                                                                                                                      <  2:73613/139‑1 (MQ=255)         ttttCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGtt                                                                                                                   >  1:1192992/1‑139 (MQ=255)              gTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAaca                                                                                                              <  2:1192992/139‑1 (MQ=255)              gTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAaca                                                                                                              >  2:181763/1‑139 (MQ=255)              gTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAaca                                                                                                              <  2:252102/139‑1 (MQ=255)               ttaattaaCACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTTCAGAGTTCAacac                                                                                                             <  1:60952/139‑1 (MQ=255)                       caccacACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTa                                                                                                     >  1:854988/1‑139 (MQ=255)                         ccGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTAcc                                                                                                   <  2:849176/139‑1 (MQ=255)                            cACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCaga                                                                                                                      <  2:216802/117‑1 (MQ=255)                            cACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCaga                                                                                                                      >  1:216802/1‑117 (MQ=255)                            cACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCg                                                                                                >  1:1237571/1‑139 (MQ=255)                              cccATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTc                                                                                              >  1:269002/1‑139 (MQ=255)                                 atcatATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCttt                                                                                           <  1:1331444/139‑1 (MQ=255)                                 atcatATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCttt                                                                                           >  1:779208/1‑139 (MQ=255)                                 atcatATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCttt                                                                                           >  2:1340924/1‑139 (MQ=255)                                 atcatATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCttt                                                                                           <  2:1055442/139‑1 (MQ=255)                                  tcatATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCtttt                                                                                          >  2:961010/1‑139 (MQ=255)                                  tcatATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCtttt                                                                                          >  1:706201/1‑139 (MQ=255)                                   catATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTg                                                                                         <  1:1186523/139‑1 (MQ=255)                                   catATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCAGCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTg                                                                                         <  2:1180969/139‑1 (MQ=255)                                     tatTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGc                                                                                       <  1:1105271/139‑1 (MQ=255)                                      atTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCg                                                                                      <  2:854988/139‑1 (MQ=255)                                      atTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCg                                                                                      <  2:824911/139‑1 (MQ=255)                                         ttGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGata                                                                                   >  2:1137968/1‑139 (MQ=255)                                         ttGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGata                                                                                   >  1:1621557/1‑139 (MQ=255)                                         ttGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGata                                                                                   >  1:88111/1‑139 (MQ=255)                                         ttGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGata                                                                                   >  2:1581050/1‑139 (MQ=255)                                           gCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAaca                                                                                                              >  2:1151267/1‑110 (MQ=255)                                           gCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAaca                                                                                                              <  1:1151267/110‑1 (MQ=255)                                           gCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTg                                                                                         >  1:987279/1‑131 (MQ=255)                                           gCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTg                                                                                         <  2:987279/131‑1 (MQ=255)                                           gCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGAtata                                                                                 <  1:986672/139‑1 (MQ=255)                                            cTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGAc                                                                                                                                  <  2:1160462/89‑1 (MQ=255)                                            cTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGAc                                                                                                                                  >  1:1160462/1‑89 (MQ=255)                                            cTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAac                                                                                                               >  2:212244/1‑108 (MQ=255)                                            cTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAac                                                                                                               <  1:212244/108‑1 (MQ=255)                                               aaGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTg                                                                                         >  1:1292697/1‑127 (MQ=255)                                               aaGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTg                                                                                         <  2:1292697/127‑1 (MQ=255)                                                 gggCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTAcc                                                                                                   <  1:234076/115‑1 (MQ=255)                                                 gggCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTAcc                                                                                                   >  2:234076/1‑115 (MQ=255)                                                     ggggAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCa                                                                       <  1:1137968/139‑1 (MQ=255)                                                      gggAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCAt                                                                      >  1:1205064/1‑139 (MQ=255)                                                      gggAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCAt                                                                      >  2:494033/1‑139 (MQ=255)                                                       ggAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATc                                                                     >  2:631390/1‑139 (MQ=255)                                                        gAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCg                                                                    <  1:181763/139‑1 (MQ=255)                                                        gAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCg                                                                    <  2:1205064/139‑1 (MQ=255)                                                         aaaaaaaTGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCgg                                                                   <  1:1340924/139‑1 (MQ=255)                                                               aTGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTAttt                                                             <  1:1581050/139‑1 (MQ=255)                                                                 gCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGAc                                                                                                                                  >  2:432566/1‑68 (MQ=255)                                                                 gCCATCGGCAGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGAc                                                                                                                                  <  1:432566/68‑1 (MQ=255)                                                                   cATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAAt                                                         <  2:77915/139‑1 (MQ=255)                                                                       ggcggcATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGTCTATTTGAATGAgg                                                     <  1:804792/139‑1 (MQ=255)                                                                       ggcggcATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAgg                                                     <  1:961010/139‑1 (MQ=255)                                                                       ggcggcATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAgg                                                     <  1:631390/139‑1 (MQ=255)                                                                           gcATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCt                                                 <  2:779208/139‑1 (MQ=255)                                                                                       aTCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGa                                                       <  1:1660342/122‑1 (MQ=255)                                                                                       aTCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGa                                                       >  2:1660342/1‑122 (MQ=255)                                                                                         cGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATg                                    >  2:602551/1‑139 (MQ=255)                                                                                          ggTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGa                                   <  1:350033/139‑1 (MQ=255)                                                                                           gTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTgg                                                                                        <  2:1372097/85‑1 (MQ=255)                                                                                           gTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTgg                                                                                        >  1:1372097/1‑85 (MQ=255)                                                                                           gTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGAc                                  >  2:1365857/1‑139 (MQ=255)                                                                                            ttAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTTCAGAGTTCAACACGGTTTCTACCTCGTCTTTTg                                                                                         <  1:531727/83‑1 (MQ=255)                                                                                            ttAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTg                                                                                         >  2:531727/1‑83 (MQ=255)                                                                                                   cATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACaga                                                                             <  1:1590954/88‑1 (MQ=255)                                                                                                   cATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACaga                                                                             >  2:1590954/1‑88 (MQ=255)                                                                                                     tGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACAt                        >  1:1399926/1‑139 (MQ=255)                                                                                                        ttACCTCTGCATCCTGGTGCATTTAGAGACACACTTTCCAGAGATCAGCCTGGTTTCTACCTCGTATATTGGCGTTATACAGAGCGTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATggc                     <  2:269002/139‑1 (MQ=255)                                                                                                            ctctGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCa                 <  2:1399926/139‑1 (MQ=255)                                                                                                              ctGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCAcc               <  1:602551/139‑1 (MQ=255)                                                                                                               tGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCg                                                                                      >  1:1619633/1‑67 (MQ=255)                                                                                                               tGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCg                                                                                      <  2:1619633/67‑1 (MQ=255)                                                                                                               tGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCg              >  2:1453215/1‑139 (MQ=255)                                                                                                                 cATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGAt            >  2:1468956/1‑139 (MQ=255)                                                                                                                  aTCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATg           <  2:1374068/139‑1 (MQ=255)                                                                                                                   tCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCg                                                                                      <  2:1216424/63‑1 (MQ=255)                                                                                                                   tCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCg                                                                                      >  1:1216424/1‑63 (MQ=255)                                                                                                                    ccTGGCGCATTTAAAGACTGGATTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGa         <  2:1237571/139‑1 (MQ=255)                                                                                                                      tGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGaaa       >  1:1476480/1‑139 (MQ=255)                                                                                                                         cgcATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAgg                                                     >  2:1720/1‑90 (MQ=255)                                                                                                                         cgcATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAgg                                                     <  1:1720/90‑1 (MQ=255)                                                                                                                         cgcATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTg    >  1:1181629/1‑139 (MQ=255)                                                                                                                         cgcATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTg    >  1:218025/1‑139 (MQ=255)                                                                                                                          gcATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTga   <  1:1365857/139‑1 (MQ=255)                                                                                                                           cATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTgag  <  2:357577/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                        CTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTAT‑CCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAG  >  NC_000913/1565634‑1565893

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