Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,103,566	T→A	100%	M185K (ATG→AAG)	mutY →	adenine DNA glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,103,566	0	T	A	100.0%	293.0 / NA	80	M185K (ATG→AAG)	mutY	adenine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (41/39);  total (41/39)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGCCGGAGCGCACCGGC  >  NC_000913/3103421‑3103711                                                                                                                                                  |                                                                                                                                                  ctcGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGa                                                                                                                                                >  4:679576/1‑149 (MQ=255) ctcGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTACGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGa                                                                                                                                                >  4:867581/1‑149 (MQ=255)  tcGACGGTAACGTCAAACGCGTGCTGGCTCGCTCCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGGGGTCGAGAATAAATTATGGAGTCTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGag                                                                                                                                               >  4:893763/1‑148 (MQ=255)   cGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGAtt                                                                                                                                              <  2:196680/149‑1 (MQ=255)      cGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTgg                                                                                                                                           <  4:307181/149‑1 (MQ=255)      cGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTgg                                                                                                                                           <  1:706043/149‑1 (MQ=255)         tAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTg                                                                                                                                        >  1:191265/1‑149 (MQ=255)         tAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCCGGGAAAAAAGAGGTCGAGAATAAATTCTGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTg                                                                                                                                        >  4:31215/1‑149 (MQ=255)            tgTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCga                                                                                                                                     <  1:303186/148‑1 (MQ=255)             gTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCgat                                                                                                                                    >  4:172857/1‑149 (MQ=255)             gTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCgat                                                                                                                                    <  4:274985/149‑1 (MQ=255)               cAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCgatga                                                                                                                                  <  2:229187/149‑1 (MQ=255)                aaaCGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCgatgat                                                                                                                                 <  3:31215/149‑1 (MQ=255)                aaaCGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCgatgat                                                                                                                                 <  4:173180/149‑1 (MQ=255)                aaaCGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCgatgat                                                                                                                                 <  4:335052/149‑1 (MQ=255)                aaaCGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCgatgat                                                                                                                                 <  4:570043/149‑1 (MQ=255)                aaaCGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCgatgat                                                                                                                                 <  4:811151/149‑1 (MQ=255)                 aaCGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGAtt                                                                                                                                <  1:620375/149‑1 (MQ=255)                  aCGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGAttt                                                                                                                               <  3:89149/149‑1 (MQ=255)                       tGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTAcg                                                                                                                          >  1:682053/1‑149 (MQ=255)                                  tgctATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGaaa                                                                                                                 >  2:425310/1‑147 (MQ=255)                                  tgctATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGaaa                                                                                                                 <  1:425310/147‑1 (MQ=255)                                     tATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGaa                                                                                                            >  4:491952/1‑149 (MQ=255)                                      aTGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGaaa                                                                                                           >  3:920569/1‑149 (MQ=255)                                      aTGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGaaa                                                                                                           <  2:59784/149‑1 (MQ=255)                                      aTGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGaaa                                                                                                           <  4:232663/149‑1 (MQ=255)                                      aTGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGaaa                                                                                                           <  2:386109/149‑1 (MQ=255)                                       tGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAAt                                                                                                          <  4:517440/149‑1 (MQ=255)                                          tGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGtt                                                                                                       >  1:11238/1‑149 (MQ=255)                                             aaGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGc                                                                                                    <  2:730537/149‑1 (MQ=255)                                                cGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAACCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGctct                                                                                                 <  4:588669/149‑1 (MQ=255)                                                 ggCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTg                                                                                                >  3:63117/1‑149 (MQ=255)                                                   cTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTc                                                                                              <  1:735655/149‑1 (MQ=255)                                                      gCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGc                                                                                           >  3:35945/1‑149 (MQ=255)                                                        cTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTa                                                                                         >  4:51275/1‑149 (MQ=255)                                                         tGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTg                                                                                                                              <  1:871755/111‑1 (MQ=255)                                                         tGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTg                                                                                                                              >  2:871755/1‑111 (MQ=255)                                                          gggAAAAAAGGGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTTATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACa                                                                                       >  4:482755/1‑149 (MQ=255)                                                           ggAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGaaa                                                                                                                 >  1:156694/1‑122 (MQ=255)                                                           ggAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGaaa                                                                                                                 <  2:156694/122‑1 (MQ=255)                                                                aaaGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACgg                                                                                 >  3:70457/1‑149 (MQ=255)                                                                 aaGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGa                                                                                <  1:254360/149‑1 (MQ=255)                                                                  agagGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGAt                                                                               <  3:482755/149‑1 (MQ=255)                                                                  agagGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGAt                                                                               <  4:914223/149‑1 (MQ=255)                                                                  agagGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGAt                                                                               <  2:11238/149‑1 (MQ=255)                                                                  agagGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGAt                                                                               >  2:340941/1‑149 (MQ=255)                                                                       tCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTAtt                                                                          >  4:378200/1‑149 (MQ=255)                                                                        cGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTg                                                                         <  4:124640/149‑1 (MQ=255)                                                                         gagaATAAATTTTGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTgc                                                                        >  2:922647/1‑149 (MQ=255)                                                                                     tGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCcaacaa                                                            <  2:386750/149‑1 (MQ=255)                                                                                       gAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATa                                                          >  4:509493/1‑149 (MQ=255)                                                                                         gTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGc                                                        >  1:471873/1‑149 (MQ=255)                                                                                          tttGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCt                                                       >  4:66540/1‑149 (MQ=255)                                                                                           ttGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTg                                                      <  4:35945/149‑1 (MQ=255)                                                                                            tGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTgg                                                     <  4:70457/149‑1 (MQ=255)                                                                                                       gTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATcc                                          >  1:420237/1‑149 (MQ=255)                                                                                                         gACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCgg                                        <  2:119056/149‑1 (MQ=255)                                                                                                           cGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGc                                      <  1:232320/149‑1 (MQ=255)                                                                                                              ccGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCaaa                                   >  1:744199/1‑149 (MQ=255)                                                                                                               cgcgGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTgccgcc                                                                    >  4:200629/1‑115 (MQ=255)                                                                                                               cgcgGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTgccgcc                                                                    <  3:200629/115‑1 (MQ=255)                                                                                                                 cgGTTGGCGTGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCaaaaaa                                <  2:632654/149‑1 (MQ=255)                                                                                                                          tGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGa                            >  4:675675/1‑144 (MQ=255)                                                                                                                          tGGAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGa                            <  3:675675/144‑1 (MQ=255)                                                                                                                           ggAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGa                      <  2:905948/149‑1 (MQ=255)                                                                                                                            gAACGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGATACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGAc                     >  4:335013/1‑149 (MQ=255)                                                                                                                             aaCGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGTAAAAAACCGAAACAGACg                    >  3:583717/1‑149 (MQ=255)                                                                                                                             aaCGGTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACg                    >  2:243396/1‑149 (MQ=255)                                                                                                                                 gTTTAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGc                >  3:166549/1‑149 (MQ=255)                                                                                                                                  tttAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGcc               >  4:345440/1‑149 (MQ=255)                                                                                                                                  tttAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGcc               >  2:379030/1‑149 (MQ=255)                                                                                                                                  tttAATCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGcc               >  3:661325/1‑149 (MQ=255)                                                                                                                                     aaTCAGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGCCGGa            <  1:922647/149‑1 (MQ=255)                                                                                                                                         aGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGaaa                          >  2:573176/1‑131 (MQ=255)                                                                                                                                         aGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGaaa                          <  1:573176/131‑1 (MQ=255)                                                                                                                                         aGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGATACAGACGCTGCCGGAgcgc        >  4:396198/1‑149 (MQ=255)                                                                                                                                         aGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGCCGGAgcgc        >  2:736583/1‑149 (MQ=255)                                                                                                                                         aGGCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCCGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGCCGGAgcgc        >  2:246780/1‑149 (MQ=255)                                                                                                                                           gCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAGCAGACGCTGCCGGAGCGCAc      >  1:424953/1‑149 (MQ=255)                                                                                                                                           gCGATGAAGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGCCGGAGCGCAc      >  3:174253/1‑149 (MQ=255)                                                                                                                                             gatgaaGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGCCGGAGCGCACCg    >  2:9782/1‑149 (MQ=255)                                                                                                                                             gatgaaGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCACCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGCCGGAGCGCACCg    >  4:712867/1‑149 (MQ=255)                                                                                                                                               tgaagGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGCCGGAGCGCACCGGc  <  1:682676/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  CTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGCCGGAGCGCACCGGC  >  NC_000913/3103421‑3103711

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: ‑