Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,506,492	C→T	22.9%	A47A (GCG→GCA)	ypdF ←	Xaa‑Pro aminopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,506,492	0	C	T	22.9%	178.0 / 50.0	96	A47A (GCG→GCA)	ypdF	Xaa‑Pro aminopeptidase
Reads supporting (aligned to +/- strand):  ref base C (41/33);  new base T (8/14);  total (49/47)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

ATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGCACTGCATCCAGTT  >  NC_000913/2506394‑2506589                                                                                                   |                                                                                                  aTAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCGGAATGTgcg                                                                                                  >  2:1153970/1‑100 (MQ=255)   aGTGGTAAGCGTGTTCGTCGCGTCAAGCATATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACATAAACGTAATAGCGCGAATCCACCAGAATGTGCGCa                                                                                                <  1:1394529/100‑1 (MQ=255)   aGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCa                                                                                                <  1:1100414/100‑1 (MQ=255)     tGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACt                                                                                              >  1:859146/1‑100 (MQ=255)      ggTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACtt                                                                                             >  2:27365/1‑100 (MQ=255)        tAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCCCGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTc                                                                                           >  2:1746214/1‑100 (MQ=255)         aaGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCa                                                                                          >  2:2220518/1‑100 (MQ=255)          aGCGTGTTCGTCGCGTCAAACAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCAAATCCACCAGAATGTGCGCACTTTCAc                                                                                         >  2:890362/1‑100 (MQ=255)             gtgtTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGTGCACTTTCACGGc                                                                                      <  1:1662507/100‑1 (MQ=255)             gtgtTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGc                                                                                      <  2:1671760/100‑1 (MQ=255)              tgtTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCt                                                                                     <  1:159570/100‑1 (MQ=255)                  cgtcgCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGTGCACTTTCACGGCTAATc                                                                                 <  2:51356/100‑1 (MQ=255)                  cgtcgCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCAAATCCACCAGAATGTGTGCACTTTCACAGCTAATc                                                                                 <  2:1536744/100‑1 (MQ=255)                   gtcgCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGTGCACTTTCACGGCTAATca                                                                                <  1:1593995/100‑1 (MQ=255)                   gtcgCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATca                                                                                >  2:2047636/1‑100 (MQ=255)                   gtcgCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCAACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATca                                                                                <  1:926017/100‑1 (MQ=255)                      gcgTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATcacca                                                                             <  1:549599/100‑1 (MQ=255)                      gcgTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATcacc                                                                              >  2:1849733/1‑99 (MQ=255)                      gcgTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCAACAGAATGTGCGCACTTTCACGGCTAATcacca                                                                             >  2:916096/1‑100 (MQ=255)                       cgTCAAGCAACGGCAGCTGGTCGTCGTGCCGACGGGCCTCCTCCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATcac                                                                               <  2:2095308/97‑1 (MQ=255)                        gTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATATGCACACTTTCACGGCTAATCacaaca                                                                           >  1:798637/1‑100 (MQ=255)                          cAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATa                                                                         >  2:1274419/1‑100 (MQ=255)                           aaGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGTGCACTTTCACGGCTAATCACCACATaa                                                                        >  1:1239883/1‑100 (MQ=255)                            aGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAAc                                                                       <  2:1493158/100‑1 (MQ=255)                            aGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAAc                                                                       >  2:430670/1‑100 (MQ=255)                              cAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATaa                                                                        >  1:1951397/1‑97 (MQ=255)                              cAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCt                                                                     >  1:816288/1‑100 (MQ=255)                                    cAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCg                                                               >  1:586440/1‑100 (MQ=255)                                      gCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGt                                                             >  1:1878847/1‑100 (MQ=255)                                       cTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTg                                                            <  2:385210/100‑1 (MQ=255)                                       cTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTg                                                            <  2:385724/100‑1 (MQ=255)                                        tGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTgg                                                           <  1:526317/100‑1 (MQ=255)                                         ggTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTAATCCGATGga                                                          >  1:216839/1‑100 (MQ=255)                                         ggTAGCCTTGCGCACGGGCTTCCACATCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTACTCACCACATAACCTGCTCCGGTGga                                                          >  2:918905/1‑100 (MQ=255)                                          gTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGTGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGgag                                                         <  1:167490/100‑1 (MQ=255)                                          gTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGgag                                                         <  1:27365/100‑1 (MQ=255)                                          gTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGgag                                                         <  2:843071/100‑1 (MQ=255)                                           tAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTCTCACGGCTAATCACCACATAACCTGATCCGGTGgaga                                                        <  2:1191947/100‑1 (MQ=255)                                            aGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGTGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGAt                                                       <  1:1118294/100‑1 (MQ=255)                                             gCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGTGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATc                                                      >  1:877713/1‑100 (MQ=255)                                             gCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATc                                                      >  1:1126808/1‑100 (MQ=255)                                             gCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCAGTGGAGAt                                                       >  1:751660/1‑99 (MQ=255)                                              ccTTGCGTACGGTCTTCCACCTCAACGTAATAGCGCGAAACCACCAGAATGTACGCATTTTTACGGTTATTTACCACATAACCTAAGCCGATg                                                            <  2:1097167/93‑1 (MQ=255)                                               cTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGCAGATcc                                                     >  1:659940/1‑99 (MQ=255)                                                 tccgcACGCGCTTCCACTTCAACTTAATACCGCTAATCCACCAGAATGTGCGCACTATCACTGCTAATCACCACATAACCTGATCCGGTGgag                                                         <  2:1183199/91‑1 (MQ=255)                                                 tGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGCTCCCCa                                                  >  1:340212/1‑100 (MQ=255)                                                  gcgcACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAg                                                 <  2:670449/100‑1 (MQ=255)                                                   cgcACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGTGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGa                                                >  2:1160730/1‑100 (MQ=255)                                                   cgcACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGa                                                >  1:117975/1‑100 (MQ=255)                                                    gcACGGGCTTCCACCTCAACGTAATAGCGCGAAACCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGAGCCGGTGGAGATCCCCAGa                                                >  2:42731/1‑99 (MQ=255)                                                     cACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATg                                              >  2:1116814/1‑100 (MQ=255)                                                       cGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCg                                            >  2:435446/1‑100 (MQ=255)                                                        gggCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGTGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCgg                                           <  1:1949163/100‑1 (MQ=255)                                                         ggCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGc                                          >  1:513508/1‑100 (MQ=255)                                                         ggCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGc                                          >  2:1688603/1‑100 (MQ=255)                                                          gCTTCCCCCTCAACGTAATAGCACGAATCCACCAGAATGCGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGc                                          >  2:887822/1‑99 (MQ=255)                                                          gCTTCCACCTCAACGTAATAGTGCGAATCCACCAGAATGTGTGCCCTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCt                                         <  2:773873/100‑1 (MQ=255)                                                           cTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCtg                                        >  2:49852/1‑100 (MQ=255)                                                              ccACCTCAACGTAATAGCGCGAATCCAACAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCtgttt                                     >  2:248980/1‑100 (MQ=255)                                                               cACCTCAACGTAATAGCGCGAATCCACCAGAATGTGTGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCtgtttg                                    <  1:1165946/100‑1 (MQ=255)                                                               cACCTCAACGTAATAGCGCGAATCCACCAGAATGTGTGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCtgtttg                                    >  1:506109/1‑100 (MQ=255)                                                               cACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCtgtttg                                    >  2:777779/1‑100 (MQ=255)                                                               cACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCtgtttg                                    <  2:1286492/100‑1 (MQ=255)                                                                aCCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCtgtttgt                                   >  2:438912/1‑100 (MQ=255)                                                                 ccTCAACGTAATAGCGCGAATCCACCAGAATGTGTGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCtgtttgtt                                  >  2:790008/1‑100 (MQ=255)                                                                 ccTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCtgtttgtt                                  >  2:643367/1‑100 (MQ=255)                                                                  cTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTc                                 <  1:925639/100‑1 (MQ=255)                                                                   tCAACGTAATAGCGCGAATCCACCATAATGTGCGCACTTTCACGGCTAATCGCCACATAACCTGATCCGGTGGAGTTCCCCAGATGCGGCTGTTTGTTc                                 <  2:816319/99‑1 (MQ=255)                                                                   tCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCtgtttgttn                                 <  2:1770869/99‑2 (MQ=255)                                                                   tCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCt                                <  2:1951980/100‑1 (MQ=255)                                                                      aCGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGcc                             >  1:2130170/1‑100 (MQ=255)                                                                       cGTAATAGCGCGAATCCACCAGAATGTGTGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGcc                             >  2:653057/1‑99 (MQ=255)                                                                       cGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCg                            <  2:521642/100‑1 (MQ=255)                                                                         tAATTGCGCGAATCCACCAGAATGTGTGCACTTTCACGGCTAAGCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTg                          <  1:71622/100‑1 (MQ=255)                                                                           aTAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTgag                        <  2:1985382/100‑1 (MQ=255)                                                                           aTAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTgag                        <  1:1209861/100‑1 (MQ=255)                                                                            tAGCGCGAATCCACCAGAATGTGTGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTgagg                       <  2:1256440/100‑1 (MQ=255)                                                                            tAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTgagg                       <  1:951234/100‑1 (MQ=255)                                                                             aGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTgagga                      <  1:269979/100‑1 (MQ=255)                                                                              gcgcGAATCCACCAGAATGTGTGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTgagg                       <  2:907772/98‑1 (MQ=255)                                                                              gcgcGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTgaggag                     <  1:2220533/100‑1 (MQ=255)                                                                               cgcgAATCCACCAGAATGTGTGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGgaga                    <  1:1035431/100‑1 (MQ=255)                                                                                 cgAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGgagaga                  <  1:1587720/100‑1 (MQ=255)                                                                                 cgAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGgagaga                  <  2:513509/100‑1 (MQ=255)                                                                                 cgAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGgagaga                  <  1:1248014/100‑1 (MQ=255)                                                                                      ccaccaGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGCAc             >  1:1718099/1‑100 (MQ=255)                                                                                       caccaGAATGTGTGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGCAc             >  2:401984/1‑99 (MQ=255)                                                                                       caccaGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGCACt            >  2:546646/1‑100 (MQ=255)                                                                                       caccaGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGCACt            >  1:493305/1‑100 (MQ=255)                                                                                         ccaGAATGTGTGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGCACTGc          <  1:2090936/100‑1 (MQ=255)                                                                                         ccaGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGCACTGc          >  1:1245389/1‑100 (MQ=255)                                                                                          caGAATGTGTGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGCACTGCa         >  2:624968/1‑100 (MQ=255)                                                                                           aGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGCAc             <  1:1354611/95‑1 (MQ=255)                                                                                            gAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGCACTGCATc       >  1:1741492/1‑100 (MQ=255)                                                                                               tgtgCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGCACTGCATCCAg    <  1:1977343/100‑1 (MQ=255)                                                                                                 tgCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGCACTGCATCCAgtt  >  2:1908777/1‑100 (MQ=255)                                                                                                 tgCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGCACTGCATCCAgtt  >  1:475334/1‑100 (MQ=255)                                                                                                   |                                                                                                  ATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGCACTGCATCCAGTT  >  NC_000913/2506394‑2506589

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