Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,485,626	A→G	V109A (GTC→GCC)	pepA ←	aminopeptidase A/I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,485,626	0	A	G	97.1%	237.6 / ‑2.3	69	V109A (GTC→GCC)	pepA	aminopeptidase A/I
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (30/37);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

AGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGT  >  NC_000913/4485491‑4485764                                                                                                                                        |                                                                                                                                           aGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGc                                                                                                                                        >  1:279803/1‑140 (MQ=255) aGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTATTGCCGTTTCGACAGCCTGACGCACTTTCCTGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGTCCGc                                                                                                                                        >  2:76119/1‑140 (MQ=255)  gCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCt                                                                                                                                       <  2:279803/140‑1 (MQ=255)  gCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCt                                                                                                                                       <  2:319922/140‑1 (MQ=255)   cGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCtt                                                                                                                                      >  1:767936/1‑140 (MQ=255)     gacgacGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTcc                                                                                                                                    >  1:608835/1‑140 (MQ=255)      acgacgCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCa                                                                                                                                   <  1:862174/140‑1 (MQ=255)       cgacgCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCAt                                                                                                                                  <  2:545021/140‑1 (MQ=255)        gacgCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCAtt                                                                                                                                 <  2:1024910/140‑1 (MQ=255)        gacgCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCAtt                                                                                                                                 <  1:485015/140‑1 (MQ=255)          cgcgGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGa                                                                                                                               <  2:558644/140‑1 (MQ=255)            cgGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGc                                                                                                                             <  1:100308/140‑1 (MQ=255)            cgGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGc                                                                                                                             <  2:1063174/140‑1 (MQ=255)            cgGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGc                                                                                                                             <  1:920211/140‑1 (MQ=255)            cgGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGc                                                                                                                             <  1:869371/140‑1 (MQ=255)            cgGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTACAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGc                                                                                                                             <  1:1152960/140‑1 (MQ=255)             ggTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGcc                                                                                                                            >  2:1044322/1‑140 (MQ=255)                    tctTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCAtt                                                                                                                                 >  2:1102996/1‑128 (MQ=255)                    tctTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCAtt                                                                                                                                 <  1:1102996/128‑1 (MQ=255)                        gTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCa                                                                                                                 <  2:608835/140‑1 (MQ=255)                         ttCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCAtt                                                                                                                   <  2:782093/137‑1 (MQ=255)                         ttCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCAtt                                                                                                                   >  1:782093/1‑137 (MQ=255)                         ttCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAg                                                                                                                >  1:411845/1‑140 (MQ=255)                             ttttCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTa                                                                                                            >  2:277251/1‑140 (MQ=255)                                tCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTa                                                                                                         <  1:967192/140‑1 (MQ=255)                                        tCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGtttt                                                                                                 <  1:615679/140‑1 (MQ=255)                                                gTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATaa                                                                                         <  1:536007/140‑1 (MQ=255)                                                gTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATaa                                                                                         <  1:1044322/140‑1 (MQ=255)                                                 tAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAAc                                                                                        >  2:929483/1‑140 (MQ=255)                                                   gagCGTCTCTTTTGACGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCCGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCt                                                                                      >  2:597645/1‑140 (MQ=255)                                                     gCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGTGTATTAATGGTTTTCTGAATAACCTGc                                                                                    <  1:350407/140‑1 (MQ=255)                                                     gCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGc                                                                                    <  1:769168/140‑1 (MQ=255)                                                      cGTCTCTTTTGCCGTCTCGACAGTCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCt                                                                                   >  2:48684/1‑140 (MQ=255)                                                      cGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATaa                                                                                         <  1:409174/134‑1 (MQ=255)                                                      cGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATaa                                                                                         >  2:409174/1‑134 (MQ=255)                                                      cGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCt                                                                                   >  2:165716/1‑140 (MQ=255)                                                      cGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGaa                                                                                            >  2:564341/1‑131 (MQ=255)                                                      cGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGaa                                                                                            <  1:564341/131‑1 (MQ=255)                                                         ctctTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGt                                                                                >  1:815181/1‑140 (MQ=255)                                                          tctTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTa                                                                               >  2:55701/1‑140 (MQ=255)                                                           ctTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTAc                                                                              >  2:596878/1‑140 (MQ=255)                                                           ctTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTAc                                                                              >  1:542700/1‑140 (MQ=255)                                                            ttttGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACt                                                                             <  1:165716/140‑1 (MQ=255)                                                              ttGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGa                                                                           >  1:275655/1‑140 (MQ=255)                                                                     ctcGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCAt                                                                    <  2:815181/140‑1 (MQ=255)                                                                       cGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATcc                                                                  <  1:76119/140‑1 (MQ=255)                                                                         aCAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAg                                                                <  1:100496/140‑1 (MQ=255)                                                                             ccTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCa                                                            <  1:597645/140‑1 (MQ=255)                                                                                gACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGt                                                         >  1:186857/1‑140 (MQ=255)                                                                                        ttCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGcc                                                 >  2:239648/1‑140 (MQ=255)                                                                                            agtagtTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCaa                                             <  2:275655/140‑1 (MQ=255)                                                                                              tagtTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAAcc                                           <  2:411845/140‑1 (MQ=255)                                                                                              tagtTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAAcc                                           <  2:542700/140‑1 (MQ=255)                                                                                              tagtTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGATTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAAcc                                           >  2:340999/1‑140 (MQ=255)                                                                                                 ttgttACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAAt                                        >  2:505030/1‑140 (MQ=255)                                                                                                   gttACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATaa                                      >  2:342492/1‑140 (MQ=255)                                                                                                   gttACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATaa                                      >  2:21608/1‑140 (MQ=255)                                                                                                   gttACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATaa                                      >  1:24288/1‑140 (MQ=255)                                                                                                    ttACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAg                                     <  1:929483/140‑1 (MQ=255)                                                                                                        ggCCTTTAACGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGgaga                                 >  1:897191/1‑140 (MQ=255)                                                                                                               aaCGTGCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCt                          >  2:969589/1‑140 (MQ=255)                                                                                                                    gCAGCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGaa                     >  2:53978/1‑140 (MQ=255)                                                                                                                       gCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGt                  <  1:565975/140‑1 (MQ=255)                                                                                                                       gCTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGt                  <  2:897191/140‑1 (MQ=255)                                                                                                                        cTCAGTCAGAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTa                 >  1:1081957/1‑140 (MQ=255)                                                                                                                             tcagAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTc            >  2:62945/1‑140 (MQ=255)                                                                                                                              cagAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCg           <  1:55701/140‑1 (MQ=255)                                                                                                                                gAAAGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGa         <  1:48684/140‑1 (MQ=255)                                                                                                                                 aaaGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGaa        <  2:1081957/140‑1 (MQ=255)                                                                                                                                 aaaGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGaa        <  1:596878/140‑1 (MQ=255)                                                                                                                                   aGCAGGCCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACa      >  2:300976/1‑140 (MQ=255)                                                                                                                                       gACCGCTTCCATTGAACCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGt  >  2:202928/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                           AGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGT  >  NC_000913/4485491‑4485764

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