Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,656,129	C→T	G77D (GGT→GAT)	pepB ←	aminopeptidase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,656,129	0	C	T	95.7%	141.1 / ‑1.9	46	G77D (GGT→GAT)	pepB	aminopeptidase B
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (13/31);  total (14/32)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

CCCCATATAACCTTGCTCACGCAGATCTTCGCCTTTGGTGATCCGATAAGTCACACGATCGCCCGCGACGTTGCTGATCAGATCAACAGCACGCTGTGCCAGTTGCGATGGTCCCAATTCTTCTGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAACCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTA  >  NC_000913/2655882‑2656362                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                          ccccATATAACCTTGCTCACGCAGATCTTCGCCTTTGGTGATCCGATAAGTCACACGATCGCCCGCGACGTTGCTGATCAGATCAACAGCACGCTGTGCCAGTTGCGATGGTCCCAATTCTTCTGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCtt                                                                                                                                                                                                                                        <  1:364429/251‑1 (MQ=255)              tGCTCACGCAGATCTTCGCCTTTGGTGATCCGATAAGTCACACGATCGCCCGCGACGTTGCTGATCAGATCAACAGCACGCTGTGCCAGTTGCGATGGTCCCAATTCTTCTGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGccc                                                                                                                                                                                                                            >  1:240389/1‑250 (MQ=255)               gCTCACGCAGATCTTCGCCTTTGGTGATCCGATAAGTCACACGATCGCCCGCGACGTTGCTGATCAGATCAACAGCACGCTGTGCCAGTTGCGATGGTCCCAATTCTTCTGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGccc                                                                                                                                                                                                                            >  1:439966/1‑249 (MQ=255)                     gCAGATCTTCGCCTTTGGTGATCCGATAAGTCACACGATCGCCCGCGACGTTGCTGATCAGATCAACAGCTCGCTGTGCCAGTTGCGATGGTCCCAATTCTTCTGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTAGCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCGCGTTACGCGTGCCTGTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGa                                                                                                                                                                                                                    <  2:411914/251‑1 (MQ=255)                       aGATCTTCGCCTTTGGTGATCCGATAAGTCACACGATCGCCCGCGACGTTGCTGATCAGATCAACAGCACGCTGTGCCAGTTGCGATGGTCCCAATTCTTCTGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCGGCGATc                                                                                                                                                                                                                  <  1:119376/251‑1 (MQ=255)                             tCGCCTTTGGTGATCCGATAAGTCACACGATCGCCCGCGACGTTGCTGATCAGATCAACAGCACGCTGTGCCAGTTGCGATGGTCCCAATTCTTCTGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATcc                                                                                                                                                                                                                 <  1:110618/246‑1 (MQ=255)                                                aaGTCACCCGATCGCCCGCGACGTTGCTGATCAGATCAACAGCACGCTGTGCCAGTTGCGATGGTCCCAATTCTTCTGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTaa                                                                                                                                                                                          <  2:3706/250‑1 (MQ=255)                                                      cacGATCGCCCGAGCAGTTGCTGATCAGATCACCAGAACGCTGTGCCAGTTGCGATGGCCCCAATTCTTCTGCCGGAGCGTTGATGCCGTAACGACCCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTACAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCCGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAAc                                                                                                                                                                                    >  2:43430/1‑250 (MQ=255)                                                        cgctCGCCCGCGACGTTGCTGATCAGATCAACCGCACGCTGTGCCAGTTGCGCTGGTCCCAATTCTTCTGCCGGTGCGTTGATGGAGTAACGCAACCAGTCGGTGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACAACACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACAtt                                                                                                                                                                                 >  2:308276/4‑250 (MQ=255)                                                             gCCCGCGACGTTGCTGATCATATCACCAGCACGCTGTGCCAGTTGCGATGGTCCCAATTCTTCTGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGa                                                                                                                                                                            <  2:241726/251‑1 (MQ=255)                                                              cccGCGACGTTGCTGATCAGATCAACAGCACGCTGTGCCAGTTGCGATGGTCCCAATTCTTCTGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGAt                                                                                                                                                                           <  1:401684/251‑1 (MQ=255)                                                                  cgACGTTGCTGATCAGATCAACAGCACGCTGTGCCAGTTGCGATGGTCCCAATTCTTCTGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGccc                                                                                                                                                                       <  1:63115/251‑1 (MQ=255)                                                                      gTTGCTGATCAGATCAACAGCACGCTGTGCCAGTTGCGATGGTCCCAATTCTTCTGCCGGTGCGTTGATGGTCTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCCCACCACTTTACCCGTGCCTTTCGGGGCTTTGAAATTTTGCCAGAATGCCCGGCAGCGATCGGCATCCAAGCCTTCACCGCTAAACGAAACATCCTTGATGccaggac                                                                                                                                                                   >  1:399443/1‑246 (MQ=255)                                                                         gCTGATCAGATCAACAGCACGCTGTGCCAGTTGCGATGGTCCCAATTCTTCTGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGt                                                                                                                                                                >  1:362879/1‑251 (MQ=255)                                                                            gatcagatcaACAGCACGCTGTGCCAGTTGCGATGGTCCCAATTCTTCTGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAACCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCa                                                                                                                                                              >  1:196664/1‑250 (MQ=255)                                                                            gatcagatcaACAGCACGCTGTGCCAGTTGCGATGGTCCCAATTCTTCTGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCAGCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCa                                                                                                                                                              <  1:174534/250‑1 (MQ=255)                                                                               cagatcaACAGCACGCTGTGCCAGTTGCGATGGTCCCAATTCTTCTGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCTCTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCt                                                                                                                                                          <  1:151654/251‑1 (MQ=255)                                                                               cagatcaACAGCACGCTGTGCCAGTTGCGATGGTCCCAATTCTTCTGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATc                                                                                                                                                           <  1:133275/250‑1 (MQ=255)                                                                                    caACAGCACGCTGTGCCAGTTGCGATGGTCCCAATTCTTCTGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGcc                                                                                                                                                                        >  2:397267/1‑232 (MQ=255)                                                                                    caACAGCACGCTGTGCCAGTTGCGATGGTCCCAATTCTTCTGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGcc                                                                                                                                                                        <  1:397267/232‑1 (MQ=255)                                                                                                           gATGGTCCCAATTCTTCTGCCGGGGCGTTGATGGTGTCACGCACCCAGTCGATGATCAGCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACGTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTGGATGCCCAGGCCGGCAATCTTGCGCGCCGCACGCTGGTTCAGCCCCAg                                                                                                                              <  2:45198/251‑1 (MQ=255)                                                                                                              ggTCCCAATTCTTCTGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATACGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCg                                                                                                                                                                                                                     >  2:372905/1‑161 (MQ=255)                                                                                                              ggTCCCAATTCTTCTGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGAGCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCg                                                                                                                                                                                                                     <  1:372905/161‑1 (MQ=255)                                                                                                                          tctGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCg                                                                                                                 <  1:176785/249‑1 (MQ=255)                                                                                                                             gCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGtt                                                                                                               <  1:44937/248‑1 (MQ=255)                                                                                                                              ccGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCaaa                                                                                                           <  1:161271/251‑1 (MQ=255)                                                                                                                                   gCGTTGATGTTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCTCTGTGCATCGTCCAGATCCGGCCCCACCACGTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTGCACCGCTTAACTGAACATGCTTGATGGCCAGACCTTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCa                                                                                                       <  2:79644/250‑1 (MQ=255)                                                                                                                                         aTGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAAt                                                                                                <  1:153238/251‑1 (MQ=255)                                                                                                                                              gtCACGCACCCTGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCTTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTACATCTTGCGCGCCGCACGCTGGATCAGCCCCAGCTCGACTGCCCCGTTCAAATGCAGGGTACTGCCGt                                                                                           >  2:284955/1‑251 (MQ=255)                                                                                                                                                  cGCACCCAGTCGATGATCTGCAGGCTGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCATACCTTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCatt                                                                                       <  2:206970/251‑1 (MQ=255)                                                                                                                                                          gTCGAGGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGTATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTGGATGCCCAGACCGTCAATCTGGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGc                                                                               <  2:141739/251‑1 (MQ=255)                                                                                                                                                                gatCATCAGGCGTTTATCCATTTCCTGCCGCTTGGCATCGTCCAGATCCAGTCGCACCGCTCTACGCGTGCCTTTCGGGGCTTTGTTATCTTGCCAGAATTCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTACCTGAACATGCTTGATGCCCAGACCGTCAATTTCGCGCGCCGCACGCTGTATCAGCCCCAGATCGTCTGCTCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAg                                                                         <  2:306779/251‑1 (MQ=255)                                                                                                                                                                   catcaGGCGGATATCCAGTTCCTGGAGCTGGGCATCGTCCAGATCCGGCACCACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCCGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGGCTGCCCCGTTCAAATGCAGGGTAATGCCGTACTTATTAATGCTGTAAGTTg                                                                      >  2:59235/1‑251 (MQ=255)                                                                                                                                                                              tATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTcc                                                            <  1:329941/250‑1 (MQ=255)                                                                                                                                                                                 ccAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGAGCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAg                                                        <  1:444141/251‑1 (MQ=255)                                                                                                                                                                                      tCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCACCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAgcgcg                                                    <  2:434249/250‑1 (MQ=255)                                                                                                                                                                                       ccTGGCGCTGGGCATCGTCCAGATCCGGCCACTCTACTGTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGTCTGCACCGCTTAACTGAACATGCTTGATGCCCAGACCGGCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGcgcgcg                                                  <  1:360363/251‑1 (MQ=255)                                                                                                                                                                                               tGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAgg                                          >  1:280040/1‑251 (MQ=255)                                                                                                                                                                                                 ggCATCGTCCAGATCCGGCCACACCACCTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGtt                                        <  2:236759/251‑1 (MQ=255)                                                                                                                                                                                                                         cacTTTACGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTTTTCTTCATCg                <  1:97527/251‑1 (MQ=255)                                                                                                                                                                                                                           cTTTACGCGTGCCCTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCt              >  1:157316/1‑251 (MQ=255)                                                                                                                                                                                                                               aCGCGTGCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGGTTTTTCTCCCCAGCGCGCGTCGGCAGGGTGGGTAGAGAGGGTAATCTTCATCGCTTCt           <  1:43430/250‑1 (MQ=255)                                                                                                                                                                                                                               aCGCGTGCCTTTCGGGGCTTTGTAACCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTTCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTg          <  2:196664/251‑1 (MQ=255)                                                                                                                                                                                                                                     tccTTTCTGTGCTTTGTAATCTGGCCAGAATGCCCAGCAGCGATCCGCATCCCATCCTTCACCGCTTAACTGAACATGCTTGATTCCCAGGCCCTCAATCTTGCGCGCCGCACGCTGGATCAGCCTCAGATCGTCTTCCCCTTTCAAATGCATGGTAATGCCGTCATTAGTAATTCTGTAAGTTGCTTTTTCTCCCCTGCGCTCTTCGGCAGGTTGGGTATGGAAGCTAACCTGCATCTCTTCTGTCAtt     <  2:19365/249‑1 (MQ=255)                                                                                                                                                                                                                                     gCCTTTCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCAttt    >  2:352324/1‑251 (MQ=255)                                                                                                                                                                                                                                         ttCGGGGCTTTGTAATCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTa  <  1:305194/249‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                          CCCCATATAACCTTGCTCACGCAGATCTTCGCCTTTGGTGATCCGATAAGTCACACGATCGCCCGCGACGTTGCTGATCAGATCAACAGCACGCTGTGCCAGTTGCGATGGTCCCAATTCTTCTGCCGGTGCGTTGATGGTGTCACGCACCCAGTCGATGATCATCAGGCGGTTATCCAGTTCCTGGCGCTGGGCATCGTCCAGATCCGGCCACACCACTTTACGCGTGCCTTTCGGGGCTTTGTAACCTTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTA  >  NC_000913/2655882‑2656362

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