Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,504,799	T→A	6.1%	K250M (AAG→ATG) ‡	ypdE ←	aminopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,504,799	0	T	A	6.1%	190.0 / 2.6	66	K250M (AAG→ATG) ‡	ypdE	aminopeptidase
Reads supporting (aligned to +/- strand):  ref base T (29/33);  new base A (2/2);  total (31/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 = 6.82e-01

CCACTGTGGGCACGCCGGTGCCGGTTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATC  >  NC_000913/2504653‑2504946                                                                                                                                                   |                                                                                                                                                    ccacTGTGGGCACGCCGGTGCCGGTTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCGCGGGGGTAAGTTTGGGCGGCGCAATCAGCGACTTg                                                                                                                                                   >  2:654070/1‑149 (MQ=255)     tgtgGGCACGCCGGTGCCGGTTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGc                                                                                                                                               <  2:679004/149‑1 (MQ=255)       tgGGCACGCCGGTGGCGGTTAAGTGCACCGCCACGCCGTCCGTGCCGCCGTTGCTGAACTTATCTGCCTGCAACGGCACGCCTATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTTTCGCtt                                                                                                                                             <  2:634108/149‑1 (MQ=255)       ggggCACGCCGGTGCCGGTTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCtt                                                                                                                                             <  2:31727/148‑1 (MQ=255)        gggCACGCCGGTGCCGGTTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGATCATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACTGTTTCGACCCAGGCGGGAAGTTTTGGCGGCGCAATCAGCGACAAGTCGCtga                                                                                                                                            <  1:645366/149‑3 (MQ=255)           cACGCCGGTGCCGGTTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAaca                                                                                                                                         <  2:283497/149‑1 (MQ=255)            aCGCCGGTGCCGGTTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAacac                                                                                                                                        <  2:397926/149‑1 (MQ=255)                      cggtTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGc                                                                                                                              <  1:678465/149‑1 (MQ=255)                       ggtTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGcc                                                                                                                             >  2:1028934/1‑149 (MQ=255)                           aaGTGCACCGCCACGCCGTCAGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGACCGtt                                                                                                                         >  1:469460/1‑149 (MQ=255)                            aGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTa                                                                                                                        <  2:422992/149‑1 (MQ=255)                            aGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTa                                                                                                                        <  1:353028/149‑1 (MQ=255)                                cACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCaa                                                                                                                    >  1:984098/1‑149 (MQ=255)                                    gccccgccGTCCGTGCCGCCGTTGCTGATCATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCGCGGTTTCGACCCCGGCGGTAAGTTTTGGCGGCGCAATCAGCGCCTTGTCGCTTAACCCCAGCATCGGCCCGTTCCCAATCTg                                                                                                                >  1:662330/1‑149 (MQ=255)                                        cgccgTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGGTTCGACCCAGGCGGGATGTTTTGGCGGCGCAATCAGCGACTGGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGa                                                                                                            >  1:941606/1‑149 (MQ=255)                                             tccgtGCCGCCGTTGCTGAACACATCTGCCTGCAACGGCACGCCAATTTCTGCCGACACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGtt                                                                                                       <  2:93528/149‑1 (MQ=255)                                                gtGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGc                                                                                                    >  2:1034754/1‑149 (MQ=255)                                                 tGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGcc                                                                                                   <  2:377093/149‑1 (MQ=255)                                                   ccgccgTGGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCcgc                                                                                                 <  2:596594/149‑1 (MQ=255)                                                    cgccgTTGCTGAACATATCTGCCTGCAACGGCACGCCAAATTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGACGCTTAACCCCAGCATCGGCCCGTTACCAATCAGGCGATGGTAAGCcgcg                                                                                                >  2:590612/1‑149 (MQ=255)                                                     gccgTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCgcgc                                                                                               <  2:442800/149‑1 (MQ=255)                                                            cTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATc                                                                                        >  1:1025594/1‑149 (MQ=255)                                                             tGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCTTTACCAATCTGGCTATGGTTAGCCGCTCCATAATCa                                                                                       >  1:773735/1‑149 (MQ=255)                                                             tGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCa                                                                                       <  2:1019984/149‑1 (MQ=255)                                                               aaCATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGTGATGGTTAGCCGCGCCATAATCaaa                                                                                     <  2:585880/149‑1 (MQ=255)                                                               aaCATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCaaa                                                                                     >  1:986942/1‑149 (MQ=255)                                                                      ctgcctgcAAAGGCACGCCAATGTCCGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTc                                                                              <  2:1140900/149‑1 (MQ=255)                                                                       tgcccgcAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCg                                                                             <  1:640297/149‑1 (MQ=255)                                                                        gcctgcAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGttt                                                                                 >  1:516978/1‑144 (MQ=255)                                                                        gcctgcAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGttt                                                                                 <  2:516978/144‑1 (MQ=255)                                                                           tgcAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCca                                                                         >  1:426278/1‑149 (MQ=255)                                                                             cAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCcagc                                                                       >  1:93773/1‑149 (MQ=255)                                                                                  gCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGTCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCg                                                                  <  2:266515/149‑1 (MQ=255)                                                                                   cACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGATACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCgg                                                                 >  1:289472/1‑149 (MQ=255)                                                                                      gCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCTTCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTa                                                               <  2:355740/148‑1 (MQ=255)                                                                                      gCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCTTCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTa                                                               >  1:355740/1‑148 (MQ=255)                                                                                       ccAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATc                                                             <  1:413247/149‑1 (MQ=255)                                                                                       ccAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATc                                                             >  2:722929/1‑149 (MQ=255)                                                                                        cAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTATCACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCa                                                            >  1:331218/1‑149 (MQ=255)                                                                                         aaTTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCTCAATCAGCTAATTGTCGCTTAACACCAGCATCGGCCCGATACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCaa                                                           <  2:611448/149‑1 (MQ=255)                                                                                         aaTTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCAAAAACAAAGTTTTTCGCCCAGCAGGCGGTATCaa                                                           >  2:465442/1‑149 (MQ=255)                                                                                          aTTTGTGCCGCCACGGTTTCGACCCAGGCGGTAAGGTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCGCGTTACCAATCTGGCGAGGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAg                                                          <  1:692268/149‑1 (MQ=255)                                                                                          aTTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAg                                                          <  2:482039/149‑1 (MQ=255)                                                                                          aTTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAg                                                          <  1:1004397/149‑1 (MQ=255)                                                                                          aTTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAg                                                          <  2:146611/149‑1 (MQ=255)                                                                                               tGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACaa                                                     >  2:791603/1‑149 (MQ=255)                                                                                                gccgccACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAAt                                                    <  1:595223/149‑1 (MQ=255)                                                                                                 ccgccACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATg                                                   <  2:241104/149‑1 (MQ=255)                                                                                                    ccACGGTTTCGCCCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGCCTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCg                                                >  1:959835/1‑149 (MQ=255)                                                                                                      aCGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGAc                                              >  1:436915/1‑149 (MQ=255)                                                                                                        ggTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGCCGATGGTTAGCCGCGCCATAATCAAAGGTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGt                                            >  2:566643/1‑149 (MQ=255)                                                                                                           ttCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCg                                         >  2:440947/1‑149 (MQ=255)                                                                                                                cccAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGac                                    <  1:917958/149‑1 (MQ=255)                                                                                                                     gCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACcgc                               >  1:1006055/1‑149 (MQ=255)                                                                                                                        gTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACcgcgcg                            >  1:768305/1‑149 (MQ=255)                                                                                                                         tacgTTGTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCgg                           <  1:452522/146‑1 (MQ=255)                                                                                                                         tAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCgg                           <  1:457756/149‑1 (MQ=255)                                                                                                                             ttttGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACCCCGCGCGGGTgg                       >  1:303094/1‑149 (MQ=255)                                                                                                                                   cggcgCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGttt                 <  1:791603/149‑1 (MQ=255)                                                                                                                                    ggcgCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTATATCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTg                <  1:598130/149‑1 (MQ=255)                                                                                                                                     gcgcAATCAGCGACTTGTCGCTTAACACCAGCATCGGCGCGTTACCAATCTGGCGCTGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCTGTTTGc               >  2:651180/1‑149 (MQ=255)                                                                                                                                      cgcAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGcc              <  2:1119114/149‑1 (MQ=255)                                                                                                                                       gtgttaaGCGACAAGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGccc             >  1:797750/7‑149 (MQ=255)                                                                                                                                        tgttaaGCGACAAGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCccg            >  2:701725/6‑149 (MQ=255)                                                                                                                                            aaGCGACAAGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCcgccgc        <  2:797750/148‑1 (MQ=255)                                                                                                                                                  cTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCACCCGCGTAATc  >  1:109874/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CCACTGTGGGCACGCCGGTGCCGGTTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATC  >  NC_000913/2504653‑2504946

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