Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,559,342	G→A	5.5%	P62L (CCA→CTA)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,559,342	0	G	A	5.5%	231.3 / 2.2	72	P62L (CCA→CTA)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (35/33);  new base A (2/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.90e-01

GCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTG  >  NC_000913/4559194‑4559490                                                                                                                                                     |                                                                                                                                                     gCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTg                                                                                                                                                      <  5:193507/149‑1 (MQ=255) gCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTg                                                                                                                                                      <  5:188742/149‑1 (MQ=255)      agaTAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGTCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCag                                                                                                                                                 <  2:227821/149‑1 (MQ=255)      agaTAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTggg                                                                                                                                                    >  8:92034/1‑146 (MQ=255)      agaTAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTggg                                                                                                                                                    >  5:314322/1‑146 (MQ=255)       gaTAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCaga                                                                                                                                                <  2:22590/149‑1 (MQ=255)          agagTCGGTGCCCAGCAGACCAACCACTGACGTGACTACCGATTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCAGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGa                                                                                                                                             >  3:42194/1‑149 (MQ=255)           gagTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGAt                                                                                                                                            <  1:260409/149‑1 (MQ=255)             gTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCt                                                                                                                                          >  6:157532/1‑149 (MQ=255)                  tGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCACCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCa                                                                                                                                     >  1:281500/1‑149 (MQ=255)                    cccATCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCTTCGTGGGACCTGCTTCGCGGCCACCGCCAATCAATTGGACGTGTTGATCAATAATACCTGGGCTGAGGTTCTGccgtct                                                                                                                                   <  8:204320/149‑5 (MQ=255)                       agcagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTgg                                                                                                                                                     >  8:57593/1‑128 (MQ=255)                       agcagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACt                                                                                                                                   >  1:96211/1‑146 (MQ=255)                        gcagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACt                                                                                                                                   >  6:336091/1‑145 (MQ=255)                         cagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGAt                                                                                                                              <  2:311871/149‑1 (MQ=255)                         cagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGAt                                                                                                                              <  6:314322/149‑1 (MQ=255)                          agACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATc                                                                                                                             <  6:348748/149‑1 (MQ=255)                                 ccaCTGACGTGACACCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAAcc                                                                                                                      >  5:265104/1‑149 (MQ=255)                                  caCTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCg                                                                                                                     >  8:11612/1‑149 (MQ=255)                                   aCTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCACCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGt                                                                                                                    <  2:281500/149‑1 (MQ=255)                                    cTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTg                                                                                                                   >  4:68053/1‑149 (MQ=255)                                     tgacgtgacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTg                                                                                                                   >  2:225612/1‑148 (MQ=255)                                       acgtgacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCa                                                                                                                 <  8:168693/148‑1 (MQ=255)                                              gCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCg                                                                                                            >  4:125873/1‑146 (MQ=255)                                              gCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCg                                                                                                            <  3:125873/146‑1 (MQ=255)                                              gCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTAGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTa                                                                                                         >  8:326854/1‑149 (MQ=255)                                              gCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTAGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTa                                                                                                         >  8:321916/1‑149 (MQ=255)                                                    ttCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGACGCCACCGCCAATAAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCTGTACAATGt                                                                                                   >  2:81760/1‑149 (MQ=255)                                                        gTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAg                                                                                                >  8:272310/1‑148 (MQ=255)                                                            ggCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAggg                                                                                           <  1:279290/149‑1 (MQ=255)                                                             gCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGa                                                                                          <  4:110282/149‑1 (MQ=255)                                                               gACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGata                                                                                        >  4:326393/1‑149 (MQ=255)                                                                 cTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATAtt                                                                                      <  2:96211/149‑1 (MQ=255)                                                                  ttAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTg                                                                                     >  2:214598/1‑149 (MQ=255)                                                                  ttAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTg                                                                                     >  6:67271/1‑149 (MQ=255)                                                                    aGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCt                                                                                   >  2:42354/1‑149 (MQ=255)                                                                     gcgcCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTAGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCt                                                                                   <  7:321916/148‑1 (MQ=255)                                                                     gcgcCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTAGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCt                                                                                   <  7:326854/148‑1 (MQ=255)                                                                         cACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCaa                                                                              >  8:335430/1‑149 (MQ=255)                                                                           cTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACg                                                                            >  4:343087/1‑149 (MQ=255)                                                                            ttCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACgg                                                                           >  5:7201/1‑149 (MQ=255)                                                                               cGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGa                                                                        <  7:335430/149‑1 (MQ=255)                                                                               cGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGa                                                                        <  3:343087/149‑1 (MQ=255)                                                                                   gTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGAAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAAt                                                                    <  5:208633/149‑1 (MQ=255)                                                                                         gtcgtGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGcc                                                              <  7:161925/149‑1 (MQ=255)                                                                                         gtcgtGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGcc                                                              <  5:67271/149‑1 (MQ=255)                                                                                               ggACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGc                                                        <  3:326393/149‑1 (MQ=255)                                                                                                   cTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGcgac                                                     >  7:143736/1‑148 (MQ=255)                                                                                                      cTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAg                                                  <  7:92034/148‑1 (MQ=255)                                                                                                       ttCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGa                                                >  6:116713/1‑149 (MQ=255)                                                                                                           ccgccaccgccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATc                                            <  5:157532/149‑1 (MQ=255)                                                                                                           ccgccaccgccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATc                                            <  5:336091/149‑1 (MQ=255)                                                                                                             gccaccgccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGc                                          <  6:265104/149‑1 (MQ=255)                                                                                                                  cgccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAAtt                                     >  7:189348/1‑149 (MQ=255)                                                                                                                   gccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAg                                                  <  5:215215/135‑1 (MQ=255)                                                                                                                   gccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAg                                                  >  6:215215/1‑135 (MQ=255)                                                                                                                   gccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTc                                    >  6:325205/1‑149 (MQ=255)                                                                                                                     caATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGcgacga                                                   <  8:196461/132‑1 (MQ=255)                                                                                                                     caATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGcgacga                                                   >  7:196461/1‑132 (MQ=255)                                                                                                                     caATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTccc                                  >  7:70916/1‑149 (MQ=255)                                                                                                                      aaTCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTcccc                                 >  7:110277/1‑149 (MQ=255)                                                                                                                                  gtgtTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTccc                                  >  2:170513/1‑136 (MQ=255)                                                                                                                                  gtgtTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTccc                                  <  1:170513/136‑1 (MQ=255)                                                                                                                                  gtgtTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGgcg                     <  4:198258/149‑1 (MQ=255)                                                                                                                                   tgtTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGgcgc                    <  1:42354/149‑1 (MQ=255)                                                                                                                                     ttGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCa                   >  6:66451/1‑148 (MQ=255)                                                                                                                                     ttGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCa                   >  1:52525/1‑148 (MQ=255)                                                                                                                                           aaTAAATCCTGGGCAGAGTATCTGCCCACTGAGATCGACAACCGTGAAGTTCGGTACAATGTCAGAAGTGATAATGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAAt            <  4:42194/149‑1 (MQ=255)                                                                                                                                           aaTAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAAt            <  7:274320/149‑1 (MQ=255)                                                                                                                                           aaTAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAAt            >  3:208549/1‑149 (MQ=255)                                                                                                                                            aTAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATg           >  3:114792/1‑149 (MQ=255)                                                                                                                                             tAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATgc          >  2:366359/1‑149 (MQ=255)                                                                                                                                             atctcagtGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATgc          >  8:301246/8‑149 (MQ=255)                                                                                                                                               aaaCCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATgcgc        <  3:68053/149‑1 (MQ=255)                                                                                                                                                aaCCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCt       <  7:272310/149‑1 (MQ=255)                                                                                                                                                  ccTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTcc     >  7:227602/1‑149 (MQ=255)                                                                                                                                                  ccTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTcc     <  7:11612/149‑1 (MQ=255)                                                                                                                                                     gggCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTg  >  6:191367/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTG  >  NC_000913/4559194‑4559490

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