Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,018,249	G→C	17.6%	S17R (AGC→AGG) ‡	ycbZ ←	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,018,249	0	G	C	17.6%	249.1 / 41.6	90	S17R (AGC→AGG)	ycbZ	putative peptidase
Reads supporting (aligned to +/- strand):  ref base G (35/39);  new base C (10/6);  total (45/45)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGG  >  NC_000913/1018154‑1018349                                                                                                |                                                                                                     aaCTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTAt                                                                                                  >  1:1448741/1‑100 (MQ=255)    tGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCgg                                                                                               >  1:1300737/1‑100 (MQ=255)    tGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCgg                                                                                               <  1:1838247/100‑1 (MQ=255)     gCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGt                                                                                              >  1:10653/1‑100 (MQ=255)     gCTCCAGCGCAAAGTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGt                                                                                              >  1:1749175/1‑100 (MQ=255)      cTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTa                                                                                             >  1:702936/1‑100 (MQ=255)       tCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTAt                                                                                            >  1:2562325/1‑100 (MQ=255)        ccAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATc                                                                                           <  1:187579/100‑1 (MQ=255)        ccAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATc                                                                                           >  2:1485476/1‑100 (MQ=255)          agtgcAAATTGCAGCCGCGGTTGAGTATCCCTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAg                                                                                         <  2:1302344/97‑1 (MQ=255)          aGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAg                                                                                         >  2:1334888/1‑100 (MQ=255)            cgcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGa                                                                                       <  1:987320/100‑1 (MQ=255)            cgcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGa                                                                                       >  1:890631/1‑100 (MQ=255)             gcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAGGaa                                                                                      <  2:1744074/100‑1 (MQ=38)              cAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAGGaaa                                                                                     <  2:1325591/100‑2 (MQ=38)              cAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAGGaaa                                                                                     <  2:2137296/100‑2 (MQ=38)               aaaTTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAAcc                                                                                    <  1:1626982/100‑1 (MQ=255)               aaaTTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAAcc                                                                                    <  2:1949651/100‑1 (MQ=255)                 aTTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAg                                                                                  >  1:2272947/1‑100 (MQ=255)                  ttGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAgg                                                                                 <  2:1082216/100‑1 (MQ=255)                  ttGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAgg                                                                                 <  2:1081897/100‑1 (MQ=255)                  gtGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGTAAATATTTCCTGATAGCTATCGGTATCAGGAACCAgg                                                                                 <  2:1393395/99‑1 (MQ=255)                   tggaGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGt                                                                                <  1:1674391/97‑1 (MQ=255)                    gCAGCCGCGGTTGAATATCACGGAAAAAAAGATCGTTTGCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAGGaaatatttg                                                                               <  2:1916875/100‑8 (MQ=21)                     cAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAGGaaatatttgc                                                                              <  2:2537783/100‑9 (MQ=25)                      aGCCGCGCTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCAc                                                                             >  1:1371875/1‑100 (MQ=255)                       gCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCAc                                                                             >  2:86081/1‑99 (MQ=255)                       gCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACg                                                                            <  2:2487983/100‑1 (MQ=255)                       gCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACg                                                                            >  1:877794/1‑100 (MQ=255)                       gCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACg                                                                            <  1:825821/100‑1 (MQ=255)                           cgGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCAt                                                                        <  1:1435563/100‑1 (MQ=255)                            ggTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATg                                                                       <  2:820781/100‑1 (MQ=255)                                gAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAg                                                                   <  2:167638/100‑1 (MQ=255)                                gAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAg                                                                   <  2:130044/100‑1 (MQ=255)                                 aatATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGt                                                                  >  1:1754245/3‑100 (MQ=255)                                 aGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGt                                                                  <  2:778686/100‑1 (MQ=255)                                   tATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGttt                                                                <  1:726085/100‑1 (MQ=255)                                   tATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGttt                                                                >  1:1219259/1‑100 (MQ=255)                                     tCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCg                                                              <  1:1358373/100‑1 (MQ=255)                                      cACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGGCACGCCATGCAAGTTTCGt                                                             <  2:404917/100‑1 (MQ=255)                                       aCTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTa                                                            >  2:165517/1‑100 (MQ=255)                                       aCTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTa                                                            >  2:711921/1‑100 (MQ=255)                                          gAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATg                                                         >  1:1423414/1‑100 (MQ=255)                                           aaTAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATgg                                                        <  2:1748119/100‑1 (MQ=255)                                           aaTAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTGCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATgg                                                        <  2:2600459/100‑1 (MQ=255)                                           aaTAAAGGATCGTTTTCGTCAAACAAATGTGGCTGAGCAAATATCTCCTGATAGCTAGCGGGAACAGGAACCAGGTCACGCCATGCAAGTGTCGTAATgg                                                        <  2:1178204/100‑1 (MQ=255)                                             tAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTc                                                      >  2:2273491/1‑100 (MQ=255)                                               aaGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCaa                                                    <  2:2364753/100‑1 (MQ=255)                                                    tCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTg                                               >  1:2438231/1‑100 (MQ=255)                                                     cGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGa                                              <  2:1016590/100‑1 (MQ=255)                                                      gTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGAt                                             >  1:888490/1‑100 (MQ=255)                                                       ttttCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATg                                            <  2:2182742/100‑1 (MQ=255)                                                          tCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGttt                                         >  1:2491541/1‑100 (MQ=255)                                                           cGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttt                                        >  1:1881230/1‑100 (MQ=255)                                                           cGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttt                                        >  2:2046726/1‑100 (MQ=255)                                                           cGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttt                                        <  1:2273504/100‑1 (MQ=255)                                                           cGTAAATCAAATGTGGCTGAGCAAAGATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttt                                        <  2:633925/100‑1 (MQ=255)                                                            gTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGttttt                                       <  1:822345/100‑1 (MQ=255)                                                             tcaatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttttt                                      <  2:890626/100‑1 (MQ=255)                                                              caatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTa                                     >  1:1208661/1‑100 (MQ=255)                                                              caatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCAAGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTa                                     >  2:2245736/1‑100 (MQ=255)                                                              caatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTAACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTa                                     >  2:697835/1‑100 (MQ=255)                                                               aatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGTCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAg                                    <  2:917680/100‑1 (MQ=255)                                                               aatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAg                                    <  2:894247/100‑1 (MQ=255)                                                               aatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAg                                    <  1:1670483/100‑1 (MQ=255)                                                               aatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAg                                    <  2:888485/100‑1 (MQ=255)                                                               aatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAg                                    <  1:1257913/100‑1 (MQ=255)                                                                           atGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCaaa                        >  2:764055/2‑100 (MQ=255)                                                                             gAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAg                       <  1:2451360/99‑1 (MQ=255)                                                                              aGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGcc                     >  1:2401225/1‑100 (MQ=255)                                                                                gatctggttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCcgc                   <  2:2528495/93‑1 (MQ=37)                                                                                cAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCcgc                   >  2:501172/1‑100 (MQ=255)                                                                                 acctggttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCcgcg                  >  2:896480/7‑100 (MQ=37)                                                                                 aaaTATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCcgcg                  <  2:2052252/100‑1 (MQ=255)                                                                                   ctggttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCcgcg                  >  1:1819239/5‑98 (MQ=37)                                                                                    tggttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGAtt               >  1:1909717/4‑100 (MQ=38)                                                                                    tctttCCTGATACCTATCGGTATCAGGAACCCGGTCACGACATGCAAGTTTCGTAATGGTCAAGGTGGATGTTTTCTAGTCTGTTTTCAAagcagcaa                 >  2:1749228/3‑96 (MQ=255)                                                                                    tatTGCCTGATAGCTATCGGTATCCGGAAACGGGTCACGCCATGGAAGGTTCGTAATGGTCCAAGTACATGTTTTTTGGACTGGTGTCTAAGCCGCGac                <  2:699831/99‑2 (MQ=255)                                                                                     ggttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCcgcagatc              >  1:2528508/3‑95 (MQ=37)                                                                                      tttCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtat             >  2:1598653/1‑100 (MQ=255)                                                                                      gttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtat             >  2:1816321/2‑100 (MQ=38)                                                                                      gttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtat             >  1:1385636/2‑100 (MQ=38)                                                                                      gttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtat             >  1:19450/2‑100 (MQ=38)                                                                                       ttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtata            >  1:576623/1‑100 (MQ=38)                                                                                        tCCTGATACCGATAGCTATCAGGAACAAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtata            >  2:485489/1‑99 (MQ=38)                                                                                         ccTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATAcc          >  2:682821/1‑100 (MQ=255)                                                                                         ccTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATAcc          >  1:1982218/1‑100 (MQ=255)                                                                                          cTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCg         <  2:1939596/100‑1 (MQ=255)                                                                                             aTAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTaa      >  1:1633114/1‑100 (MQ=255)                                                                                              tAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAAc     >  1:176387/1‑100 (MQ=255)                                                                                                gCTATCGGTATCAGGAACC‑GGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCgg  >  2:297990/1‑100 (MQ=255)                                                                                                |                                                                                                     AACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGG  >  NC_000913/1018154‑1018349

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