Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,683,630	C→T	*663* (TAG→TAA)	yhjK ←	cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,683,630	0	C	T	100.0%	272.1 / NA	83	*663* (TAG→TAA)	yhjK	cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (43/40);  total (43/40)

CCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCA  >  NC_000913/3683494‑3683762                                                                                                                                         |                                                                                                                                      ccTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTa                                                                                                                                     >  1:1083041/1‑139 (MQ=255)                 aaCACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTaa                                                                                                                     <  1:1422397/139‑1 (MQ=255)                  acacCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAAc                                                                                                                    <  2:1393458/139‑1 (MQ=255)                          acaacaTCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCg                                                                                                            <  1:1443568/139‑1 (MQ=255)                            aacaTCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGaa                                                                                                          <  2:47883/139‑1 (MQ=255)                            aacaTCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGaa                                                                                                          <  2:1436491/139‑1 (MQ=255)                              caTCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttcttc                                                                                                                            >  2:982613/1‑119 (MQ=255)                              caTCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttcttc                                                                                                                            <  1:982613/119‑1 (MQ=255)                                tCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGAtt                                                                                                      >  1:1356571/1‑139 (MQ=255)                                 cacaGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACGAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGAttt                                                                                                     <  2:1083041/139‑1 (MQ=255)                                    aGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCaa                                                                                                  >  2:1405583/1‑139 (MQ=255)                                      aaaTGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACtt                                                                                                                                  <  1:1623143/105‑1 (MQ=255)                                      aaaTGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACtt                                                                                                                                  >  2:1623143/1‑105 (MQ=255)                                      aaaTGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCg                                                                                                            >  2:24987/1‑127 (MQ=255)                                      aaaTGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCg                                                                                                            <  1:24987/127‑1 (MQ=255)                                         tGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttcttc                                                                                                                            >  2:191436/1‑108 (MQ=255)                                         tGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttcttc                                                                                                                            <  1:191436/108‑1 (MQ=255)                                         tGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑ggg                                                                                            >  2:1592514/1‑139 (MQ=255)                                             ccTGGATCAAAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTg                                                                                        >  1:122335/1‑139 (MQ=255)                                              cTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTgg                                                                                       >  2:75758/1‑139 (MQ=255)                                              cTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTgg                                                                                       >  1:206921/1‑139 (MQ=255)                                                ggATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGa                                                                                                        <  2:665854/121‑1 (MQ=255)                                                ggATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGa                                                                                                        >  1:665854/1‑121 (MQ=255)                                                ggATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGAttt                                                                                                     <  2:838787/124‑1 (MQ=255)                                                ggATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGAttt                                                                                                     >  1:838787/1‑124 (MQ=255)                                                ggATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGc                                                                                     >  2:1405264/1‑139 (MQ=255)                                                  aTCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                             >  1:1138793/1‑98 (MQ=255)                                                  aTCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                             <  2:1138793/98‑1 (MQ=255)                                                  aTCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                             <  2:265305/98‑1 (MQ=255)                                                  aTCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                             >  1:265305/1‑98 (MQ=255)                                                     aTAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCg                                                                                                            <  2:385738/112‑1 (MQ=255)                                                     aTAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCg                                                                                                            >  1:385738/1‑112 (MQ=255)                                                     aTAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCa                                                                                >  1:1006195/1‑139 (MQ=255)                                                     aTAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCa                                                                                >  2:374489/1‑139 (MQ=255)                                                       aaTTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACtt                                                                                                                                  >  2:171017/1‑88 (MQ=255)                                                       aaTTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACtt                                                                                                                                  <  1:171017/88‑1 (MQ=255)                                                        aTTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCaaaa                                                                             >  2:647092/1‑139 (MQ=255)                                                           gAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAgg                                                                          <  1:1405264/139‑1 (MQ=255)                                                           gAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAgg                                                                          <  1:75758/139‑1 (MQ=255)                                                           gAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAgg                                                                          <  1:647092/139‑1 (MQ=255)                                                            aaaCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACtt                                                                                                                                  <  1:951338/83‑1 (MQ=255)                                                            aaaCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACtt                                                                                                                                  >  2:951338/1‑83 (MQ=255)                                                            aaaCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttc                                                                                                                               <  2:1429801/86‑1 (MQ=255)                                                            aaaCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttc                                                                                                                               >  1:1429801/1‑86 (MQ=255)                                                            aaaCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAAct                                                                                                                   >  2:426197/1‑98 (MQ=255)                                                            aaaCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAAct                                                                                                                   <  1:426197/98‑1 (MQ=255)                                                            aaaCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAAc                                                                                                                    >  1:1079137/1‑97 (MQ=255)                                                            aaaCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAAc                                                                                                                    <  1:1159761/97‑1 (MQ=255)                                                            aaaCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAAc                                                                                                                    >  2:1159761/1‑97 (MQ=255)                                                            aaaCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAAc                                                                                                                    <  2:1079137/97‑1 (MQ=255)                                                            aaaCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGaa                                                                                                          >  1:1582314/1‑107 (MQ=255)                                                            aaaCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGaa                                                                                                          <  2:1582314/107‑1 (MQ=255)                                                            aaaCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAGGa                                                                         >  2:915732/1‑139 (MQ=255)                                                              aCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAGGAAg                                                                       >  1:1585917/1‑139 (MQ=255)                                                               cAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGa                                                                                                        >  2:1350030/1‑106 (MQ=255)                                                               cAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGa                                                                                                        <  1:1350030/106‑1 (MQ=255)                                                                aaaaaGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAGGAAGcc                                                                     <  1:1405583/139‑1 (MQ=255)                                                                      gTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGc                                                               >  2:878520/1‑139 (MQ=255)                                                                                 tGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGc                                                    <  1:915732/139‑1 (MQ=255)                                                                                        gCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAGGAAGCCCTgg                                                                 <  1:1735407/119‑1 (MQ=255)                                                                                        gCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAGGAAGCCCTgg                                                                 >  2:1735407/1‑119 (MQ=255)                                                                                         cTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCtttt                                            >  2:403712/1‑139 (MQ=255)                                                                                                  aCCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                             >  2:113972/1‑50 (MQ=255)                                                                                                  aCCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                             <  2:1591953/50‑1 (MQ=255)                                                                                                  aCCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                             >  1:1591953/1‑50 (MQ=255)                                                                                                  aCCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                             <  1:113972/50‑1 (MQ=255)                                                                                                  aCCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGccagccag                                   >  2:1532966/1‑139 (MQ=255)                                                                                                       cACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTcgcg                              <  2:1006195/139‑1 (MQ=255)                                                                                                          tttttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGaa                                                                                                          >  1:696236/1‑61 (MQ=255)                                                                                                          tttttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGaa                                                                                                          <  2:696236/61‑1 (MQ=255)                                                                                                          tttttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTtg                           <  2:498344/139‑1 (MQ=255)                                                                                                            tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGAtt                                                                                                      <  1:628790/63‑1 (MQ=255)                                                                                                            tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGAtt                                                                                                      >  2:628790/1‑63 (MQ=255)                                                                                                            tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGc                                                                                 <  1:460469/83‑1 (MQ=255)                                                                                                            tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGc                                                                                 >  2:460469/1‑83 (MQ=255)                                                                                                            tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTtgtg                         >  1:610416/1‑139 (MQ=255)                                                                                                                aaGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTc                     <  2:1356571/139‑1 (MQ=255)                                                                                                                 aGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCa                    <  1:962995/139‑1 (MQ=255)                                                                                                                 aGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCa                    <  2:1585917/139‑1 (MQ=255)                                                                                                                      tGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGtct                >  2:707511/1‑139 (MQ=255)                                                                                                                               tttGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTc      >  1:812901/1‑139 (MQ=255)                                                                                                                                  ggggTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGc   >  2:1573079/1‑139 (MQ=255)                                                                                                                                   gggTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCa  <  2:122335/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      CCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATA‑GGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCA  >  NC_000913/3683494‑3683762

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