Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,201,586	G→A	Q35* (CAA→TAA)	zraP ←	zinc responsive, periplasmic protein with chaperone activity

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,201,586	0	G	A	100.0%	210.0 / NA	62	Q35* (CAA→TAA)	zraP	zinc responsive, periplasmic protein with chaperone activity
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/34);  total (28/34)

GCTAACAGGGCATTGTATTCATAACGCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTGCTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTAATCCTCCAGTGGTTGTCTTACTTCGGGT  >  NC_000913/4201456‑4201721                                                                                                                                   |                                                                                                                                        gCTAACAGGGCATTGTATTCATAACGCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACAt                                                                                                                                >  2:203550/1‑140 (MQ=255)   tAACAGGGCATTGTATTCATAACGCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATAc                                                                                                                              <  2:501765/140‑1 (MQ=255)            aTTGTATTCATAACGCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTccgc                                                                                                                     >  1:197167/1‑140 (MQ=255)                  ttCATAACACTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAg                                                                                                               <  1:506401/140‑1 (MQ=255)                       aaCGCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAAt                                                                                                          >  2:559708/1‑140 (MQ=255)                       aaCGCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAAt                                                                                                          >  1:533002/1‑140 (MQ=255)                        aCGCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTccgccg                                                                                                                   <  1:797221/130‑1 (MQ=255)                        aCGCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTccgccg                                                                                                                   >  2:797221/1‑130 (MQ=255)                        aCGCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATg                                                                                                         >  1:356031/1‑140 (MQ=255)                        aCGCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATg                                                                                                         >  1:352213/1‑140 (MQ=255)                          gCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATg                                                                                                         >  1:729820/1‑138 (MQ=255)                          gCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCa                                                                                                       >  2:78024/1‑140 (MQ=255)                          gCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCa                                                                                                       >  2:394075/1‑140 (MQ=255)                          gCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCa                                                                                                       <  2:533822/140‑1 (MQ=255)                          gCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCa                                                                                                       >  2:632216/1‑140 (MQ=255)                          gCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGACTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATg                                                                                                         <  2:729820/138‑1 (MQ=255)                            ttCGTCACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGa                                                                                                     >  2:367583/1‑140 (MQ=255)                                 cACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCg                                                                                                <  2:352213/140‑1 (MQ=255)                                 cACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCg                                                                                                <  2:150646/140‑1 (MQ=255)                                 cACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCg                                                                                                <  1:394075/140‑1 (MQ=255)                                             ttGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCa                                                                                    <  2:125688/140‑1 (MQ=255)                                                       ctgctTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAg                                                                          >  1:285875/1‑140 (MQ=255)                                                         gctTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGGGGCATTTTACTGCCACATCCCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAgcg                                                                        <  2:286285/140‑1 (MQ=255)                                                         gctTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAgcg                                                                        <  1:817124/140‑1 (MQ=255)                                                         gctTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAgcg                                                                        <  2:711941/140‑1 (MQ=255)                                                         gctTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAgcg                                                                        <  1:559708/140‑1 (MQ=255)                                                         gctTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAgcg                                                                        <  1:512158/140‑1 (MQ=255)                                                         gctTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAgcg                                                                        <  1:78024/140‑1 (MQ=255)                                                         gctTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACACTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAgcg                                                                        <  2:371132/140‑1 (MQ=255)                                                               aGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCcatca                                                                  >  2:53017/1‑140 (MQ=255)                                                               aGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCcatca                                                                  >  2:441192/1‑140 (MQ=255)                                                                    aaaaGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGa                                                                             <  2:404409/124‑1 (MQ=255)                                                                    aaaaGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGa                                                                             >  1:404409/1‑124 (MQ=255)                                                                    aaaaGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAgcgc                                                                       >  1:756964/1‑130 (MQ=255)                                                                    aaaaGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAgcgc                                                                       <  2:756964/130‑1 (MQ=255)                                                                        gTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAggg                                                         <  1:710447/140‑1 (MQ=255)                                                                                      tCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCTCATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTc                                           <  2:70178/140‑1 (MQ=255)                                                                                      tCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTc                                           >  2:735505/1‑140 (MQ=255)                                                                                          ccACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGtt                                       >  1:496241/1‑140 (MQ=255)                                                                                           cACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGttt                                      >  2:656160/1‑140 (MQ=255)                                                                                             cGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCa                                    <  1:771198/140‑1 (MQ=255)                                                                                                 gtctgtTGTTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGt                                <  2:533002/140‑1 (MQ=255)                                                                                                 ctctgtTGTTCGCTGGTCAAAGGAGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTGCCAGGGCAATTTTCGTGTTCCGTTTCAGGGt                                <  2:219334/139‑1 (MQ=255)                                                                                                    tgttgtTCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTaa                             <  1:203550/140‑1 (MQ=255)                                                                                                     gttgttCGCTGGTCAAAGGCGCGGCATTTTACTGCCCCATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTAAt                            <  1:53017/140‑1 (MQ=255)                                                                                                     gttgttCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTAAt                            <  2:534218/140‑1 (MQ=255)                                                                                                     gttgttCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACAAGGGCAATTTTCGTGTTCCGTTTCATGGTTAAt                            <  2:356031/140‑1 (MQ=255)                                                                                                          tCGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTAAtcctcc                       >  2:143612/1‑140 (MQ=255)                                                                                                           cGCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTAATCCTCCa                      >  1:787540/1‑140 (MQ=255)                                                                                                            gCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTAATCCTCCAg                     >  1:735282/1‑140 (MQ=255)                                                                                                            gCTGGTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTAATCCTCCAg                     >  2:570708/1‑140 (MQ=255)                                                                                                               ggTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTAATCCTCCAGTgg                  <  1:632216/140‑1 (MQ=255)                                                                                                               ggTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTAATCCTCCAGTgg                  <  2:747364/140‑1 (MQ=255)                                                                                                               ggTCAAAGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTAATCCTCCAGTgg                  <  1:570708/140‑1 (MQ=255)                                                                                                                    aaGGCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTAATCCTCCAGTGGTTGTc             <  1:425206/140‑1 (MQ=255)                                                                                                                      ggCGCGGCATTTTACTGCCACATCCCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTAATCCTCCAGTGGTTGTCtt           <  2:209467/140‑1 (MQ=255)                                                                                                                      ggCGCGGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTAATCCTCCAGTGGTTGTCtt           <  1:441192/140‑1 (MQ=255)                                                                                                                          cgGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTAATCCTCCAGTgg                  <  2:307330/129‑1 (MQ=255)                                                                                                                          cgGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTAATCCTCCAGTgg                  >  1:307330/1‑129 (MQ=255)                                                                                                                          cgGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTAATCCTCCAGTGGTTGTCTTACtt       >  1:576196/1‑140 (MQ=255)                                                                                                                          cgGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTAATCCTCCAGTGGTTGTCTTACtt       >  1:631102/1‑140 (MQ=255)                                                                                                                          cgGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTAATCCTCCAGTGGTTGTCTTACtt       <  1:671330/140‑1 (MQ=255)                                                                                                                          cgGCATTTTACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTAATCCTCCAGTGGTTGTCTTACtt       >  1:724872/1‑140 (MQ=255)                                                                                                                             cATTTTTCTGCCCCATACCGTAACCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCCCCATCATCACCAGTGCAATCTTCGTGTTCCGTTTCAGTGTTAATCCTCCAGTGGTTGTCTTACTTCgt    <  2:197167/140‑2 (MQ=255)                                                                                                                               ttttACTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTAATCCTCCAGTGGTTGTCTTACTTCGGGt  <  2:721726/140‑1 (MQ=255)                                                                                                                                   |                                                                                                                                        GCTAACAGGGCATTGTATTCATAACGCTTCGTCACCAGTTGCTGTTGCAGTGCGCTGCTTTGAGCGTAAAAGTCATTATGGATTTTCTGCCACGCTGTCTGTTGTTCGCTGGTCAAAGGCGCGGCATTTTGCTGCCACATACCGTGTCCGCCGTGAGCAAATGCAGATGTCGATCCCATCGCCATTGCTGAAAGCGCCATCATTACCAGGGCAATTTTCGTGTTCCGTTTCATGGTTAATCCTCCAGTGGTTGTCTTACTTCGGGT  >  NC_000913/4201456‑4201721

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