Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,054,723	Δ1 bp	coding (1081/1317 nt)	shiA →	shikimate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,054,722	0	G	.	97.3%	327.4 / 6.1	77	coding (1080/1317 nt)	shiA	shikimate transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (39/34);  minor base T (1/1);  total (41/36)
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.97e-01

GTTGAGCGCATTTCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGTGGTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGG  >  NC_000913/2054623‑2054817                                                                                                    |                                                                                                gTTGAGCGCATTTCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACAtggtg                                                                                                 >  2:179188/1‑100 (MQ=255) gTTGAGCGCATTTCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACAtggtg                                                                                                 >  2:1664509/1‑100 (MQ=255) gTTGAGCGCATTTCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACAtggtg                                                                                                 >  1:2146172/1‑100 (MQ=255)   tgggcgcATTTCCTTTCTTTATGGCGCTTGAAGAACAATCTATTTTCTGGATAGTTTTCTTCTCCATCATGCTGGCAAACACTGCGCATGACAtggtgtg                                                                                               <  2:1224062/97‑3 (MQ=255)   tGAGCGCATTTCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACAtggtgtg                                                                                               <  2:2519182/100‑3 (MQ=255)     aGCGCATTTCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGCCAtggtgtgtg                                                                                             <  2:1246995/100‑5 (MQ=255)      gcgcATTTCCTTTCTTTATGGCGCTTGAATCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTTGCAAACATTGCGCATGCCATGgt‑gtgtgt                                                                                           >  2:2325716/1‑100 (MQ=255)       cgcATTTCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGt‑gtgtgtg                                                                                          <  1:1372829/100‑1 (MQ=255)         cATTTCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑Gtgtgtgtg                                                                                        <  2:805789/100‑1 (MQ=255)         cATTTCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑Gtgtgtgtg                                                                                        >  1:619036/1‑100 (MQ=255)          aTTTCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGc                                                                                       <  2:2541061/100‑1 (MQ=255)           tttCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGca                                                                                      >  2:1849476/1‑100 (MQ=255)           tttCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGca                                                                                      <  1:712978/100‑1 (MQ=255)            ttCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGcaa                                                                                     >  1:274638/1‑100 (MQ=255)             tCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACAGTGCGCATGACATGG‑AGTGTGTGTGcaac                                                                                    <  1:2323666/100‑1 (MQ=255)              ccTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGcaaca                                                                                   >  2:304403/1‑100 (MQ=255)               ctttctttATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGcaacaa                                                                                  >  2:2505021/1‑100 (MQ=255)               ctttctttATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGcaacaa                                                                                  >  2:697456/1‑100 (MQ=255)                 ttctttATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAAcc                                                                                <  1:2621480/100‑1 (MQ=255)                  tctttATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCg                                                                               >  1:506845/1‑100 (MQ=255)                   ctttATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGa                                                                              >  1:2057707/1‑100 (MQ=255)                   ctttATGGCGCTTGAAGCACAATATATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGAGCATGACATGGT‑GTGTGTGTGCAACAACCGa                                                                              <  2:252640/100‑1 (MQ=255)                    tttATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAGACATTGCGCATGGCATGG‑CGTGTGTGTGCAACAACCGAt                                                                             >  1:2470584/1‑100 (MQ=255)                    tttATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCggt                                                                             >  1:1988617/1‑98 (MQ=255)                     ttATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATg                                                                            >  1:2566230/1‑100 (MQ=255)                     ttATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATg                                                                            >  1:1499088/1‑100 (MQ=255)                         ggCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTTCATAATGCTGGCAATCATTGCGCATCACATGGT‑GTGTGTGTGCAACAACCGATGTTTa                                                                        >  2:370746/1‑100 (MQ=255)                           cgcTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTAcc                                                                      <  1:1925374/100‑1 (MQ=255)                           cgcTTGAAGCACAATCTATTTTCTGGATAATTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTAcc                                                                      <  1:1402334/100‑1 (MQ=255)                            gcTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCg                                                                     >  2:691406/1‑100 (MQ=255)                             cTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGa                                                                    <  1:2162535/100‑1 (MQ=255)                                   gCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGtt                                                              >  1:1938669/1‑100 (MQ=255)                                     acaATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTg                                                            <  2:1224260/100‑1 (MQ=255)                                     acaATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTg                                                            <  1:1726725/100‑1 (MQ=255)                                      caATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCGAACATTGAGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTgg                                                           >  1:2509566/1‑100 (MQ=255)                                       aaTCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGt                                                          <  1:1555658/100‑1 (MQ=255)                                       aaTCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGt                                                          >  2:2057730/1‑100 (MQ=255)                                        aTCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTg                                                         >  2:82135/1‑100 (MQ=255)                                        aTCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTg                                                         <  2:655950/100‑1 (MQ=255)                                           tATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACAGATGTTTACCGAAATGTTTGGTGCCa                                                      >  1:1320981/1‑100 (MQ=255)                                            aTTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAg                                                     >  2:248235/1‑100 (MQ=255)                                            aTTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAg                                                     <  1:476241/100‑1 (MQ=255)                                             ttttCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGt                                                    <  1:1984623/100‑1 (MQ=255)                                                   ggATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGc                                              >  2:1888715/1‑100 (MQ=255)                                                    gATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCt                                             <  2:2547122/100‑1 (MQ=255)                                                    gATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCt                                             <  1:1851826/100‑1 (MQ=255)                                                     aTAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCta                                            <  2:783713/100‑1 (MQ=255)                                                        gTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAg                                         >  1:1811293/1‑100 (MQ=255)                                                        gTTTTCTGCACCACAATGCTGGCAAACATAGCGCATGACATGGT‑GTGTGTGTGCAAGAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAg                                         >  2:2142637/1‑100 (MQ=255)                                                         ttttCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGt                                        >  2:1939363/1‑100 (MQ=255)                                                           ttcttcTCCATAATGCTGGCAACAATTTCGCATGACAT‑GTTGTGTTTATGTAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTgg                                      <  2:1873608/100‑1 (MQ=255)                                                             cttctCCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGgcg                                    >  2:1098151/1‑100 (MQ=255)                                                             cttcaCTATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCCACAATCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGgcg                                    >  2:2444964/1‑100 (MQ=255)                                                               tctcCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCt                                  >  1:1528490/1‑100 (MQ=255)                                                               tctcCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCt                                  >  1:828185/1‑100 (MQ=255)                                                               tctcCATAATGCTGGCAAACATTGCGCATGACAT‑GTTGTGTGTGTGCAACAACCGATGTTCACCGAAATGTTGGGCGCCACGTATCGCTAGAGTGACGCt                                  >  2:2084474/1‑100 (MQ=255)                                                                ctcCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTg                                 <  1:705778/100‑1 (MQ=255)                                                                ctcCACAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTg                                 <  1:2501797/100‑1 (MQ=255)                                                                 tcCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTgg                                <  1:1652649/100‑1 (MQ=255)                                                                 tcCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTgg                                <  2:1913208/100‑1 (MQ=255)                                                                 tcCATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTgg                                <  1:1604807/100‑1 (MQ=255)                                                                  ccATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGa                               >  2:306504/1‑100 (MQ=255)                                                                  ccATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGGTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGa                               <  2:2065527/100‑1 (MQ=255)                                                                   cATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAg                              >  2:912234/1‑100 (MQ=255)                                                                   cATAATGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAg                              <  1:1741332/100‑1 (MQ=255)                                                                     tAATGCTGGCAAACAGTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTc                            >  2:1976489/1‑100 (MQ=255)                                                                      aaTGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCg                           >  2:2106889/1‑100 (MQ=255)                                                                        tGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGt                         <  2:115523/100‑1 (MQ=255)                                                                        tGCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGt                         <  1:198939/100‑1 (MQ=255)                                                                        ggCTGGCAAACATTGCGCAAAACATGGT‑GTGGATAAACAACATCCAATGTTTACCGAAATGTTTAGGGCCAGTTAACGATATAGTGGCGCTGGAGTCGGt                         <  2:207512/99‑1 (MQ=255)                                                                         gCTGGCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGtt                        >  1:2583455/1‑100 (MQ=255)                                                                            ggCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATc                     >  2:138173/1‑100 (MQ=255)                                                                            ggCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATc                     >  2:748784/1‑100 (MQ=255)                                                                            ggCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATc                     >  1:2294683/1‑100 (MQ=255)                                                                             gCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCa                    <  2:2234725/100‑1 (MQ=255)                                                                             gCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCa                    <  1:189802/100‑1 (MQ=255)                                                                             gCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCa                    <  1:345178/100‑1 (MQ=255)                                                                             gCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCa                    <  2:751575/100‑1 (MQ=255)                                                                             gCAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCa                    <  2:702241/100‑1 (MQ=255)                                                                              cAAACATTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAg                   >  2:479485/1‑100 (MQ=255)                                                                                   aTTGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGc              >  1:2342778/1‑100 (MQ=255)                                                                                    ttGCGCATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGcc             <  2:2353726/100‑1 (MQ=255)                                                                                    gtgagcATGACATGGT‑GCGTGCGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCCGTTCTCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGcc             <  2:56259/96‑1 (MQ=255)                                                                                     ngcgcATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAACTGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGCGTCGGTTATCAGGTTGCCa            >  1:885913/2‑100 (MQ=255)                                                                                       cgcATGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAgt          >  1:1781850/1‑100 (MQ=255)                                                                                           tGACATGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTgg      >  1:452458/1‑100 (MQ=255)                                                                                               aTGGT‑GTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTgg  <  1:1285919/100‑1 (MQ=255)                                                                                                    |                                                                                                GTTGAGCGCATTTCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGTGGTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGG  >  NC_000913/2054623‑2054817

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: ‑