Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,095,330	G→A	6.5%	P218L (CCT→CTT)	rhaA ←	L‑rhamnose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,095,330	0	G	A	6.5%	194.5 / 2.5	61	P218L (CCT→CTT)	rhaA	L‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base G (28/29);  new base A (2/2);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-01

TCAGTCGGGTGGAAGTGCCCGGCGTCCAGGCACAGCGCAGTCTGGCGGCTGGTGGCATACCCCATGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCGATGGTGTGCCGAGTTGCTCGCCAAAATAG  >  NC_000913/4095183‑4095473                                                                                                                                                    |                                                                                                                                                tCAGTCGGGTGGAAGTGCCCGGCGTCCAGGCACAGCGCAGTCTGGCGGCTGGTGGCATACCCCATGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAg                                                                                                                                                 >  4:133857/1‑148 (MQ=255)        ggTGGAAGTGCCCGGCGTCCAGGCACAGCGCAGTCTGGCGGCTGGTGGCATACCCCATGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAg                                                                                                                                          <  3:332140/148‑1 (MQ=255)        ggTGGAAGTGCCCGGCGTCCAGGCACAGCGCAGTCTGGCGGCTGGTGGCATACCCCATGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAg                                                                                                                                          <  2:200876/148‑1 (MQ=255)         gTGGAAGTGCCCGGCGTCCAGGCACAGCGCAGTCTGGCGGCTGGTGGCATACCCCATGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGc                                                                                                                                         >  4:746/1‑148 (MQ=255)             aaGTGCCCGGCGTCCAGGCACAGCGCAGTCTGGCGGCTGGTGGCATACCCCATGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTtct                                                                                                                                     <  3:133857/148‑1 (MQ=255)              agggCCCGGCGTCCAGGCACAGCGCAGGCTGCCGGCTGGTGGCATACCCCATGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTtctc                                                                                                                                    <  7:106183/145‑1 (MQ=255)                    cggcggCCAGGCACAGCGCAGTCTGGCGGCTGGGGGCATACCCCATGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGAt                                                                                                                              <  2:79940/148‑1 (MQ=255)                     ggCGTCCAGGCACAGCGCAGTCTGGCGGCTGGTGGCATACCCCATGTAAACCTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGAAATGGTGCGCAGGGTTTAGCTTCTCGCTGATc                                                                                                                             <  5:71609/148‑1 (MQ=255)                       cGTCCAGGCACAGCGCAGTCTGGCGGCTGGTGGCATACCCCATGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCAc                                                                                                                           >  1:45905/1‑148 (MQ=255)                            aGGCACAGCGCAGTCTGGCGGCTGGTGGCATACCCCATGTAAAACGCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCCAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGGTTAGCGTCTCGCTGATCACCTCAt                                                                                                                      >  6:112460/1‑148 (MQ=255)                               cacaGCGCAGTCTGGCGGCTGGTGGCATACCCCATGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATc                                                                                                                     >  1:55715/1‑146 (MQ=255)                                  aGCGCAGTCTGGCGGCTGGTGGCATACCCCATGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGtg                                                                                                                >  4:3699/1‑148 (MQ=255)                                         tCTGGCGGCTGGTGGCATACCCCATGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATTGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCcag                                                                                                         <  2:25655/148‑1 (MQ=255)                                          cTGGCGGCTGGTGGCATACCCCATGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCcag                                                                                                         >  6:101043/1‑147 (MQ=255)                                                 gtggtggCATACCCCATCTAAATCTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCACACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGTCGCt                                                                                                 <  5:305395/147‑1 (MQ=255)                                                    gtggCATACCCCATGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTg                                                                                                <  5:137534/146‑1 (MQ=255)                                                            ccccATGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCg                                                                                      <  3:746/148‑1 (MQ=255)                                                             cccATGTAAAACTCATTGGAGCCAACCGTGTAGATCTCTGCGCCAATGCCAATCAATTTGCTCTCAACGGCTTCCATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGATGCCCGCAGACGCTGACGCGGGGc                                                                                       >  2:79627/1‑146 (MQ=255)                                                              ccATGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCgag                                                                                    <  2:45905/148‑1 (MQ=255)                                                              ccATGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCCGCCAGCAGCCGCTGACGCGGGGc                                                                                       >  4:89003/1‑145 (MQ=255)                                                                  gTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCcagc                                                                                                        >  2:299648/1‑124 (MQ=255)                                                                  gTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCcagc                                                                                                        <  1:299648/124‑1 (MQ=255)                                                                  gTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACgg                                                                                >  8:125813/1‑148 (MQ=255)                                                                  gTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACgg                                                                                >  1:4570/1‑148 (MQ=255)                                                                  gTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACgg                                                                                >  2:17374/1‑148 (MQ=255)                                                                  gTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACCGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACgg                                                                                >  1:98501/1‑148 (MQ=255)                                                                    aaaaCTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGt                                                                               >  8:26667/1‑147 (MQ=255)                                                                       aCTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAAc                                                                           >  8:230662/1‑148 (MQ=255)                                                                             tGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCACGCTGATCGCCTCATCCAGTGCTGCCAGCGGACGCTGACGCGGGGCGAGACGGTCAACGGTGa                                                                      >  6:59510/1‑147 (MQ=255)                                                                             tGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAAGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGa                                                                      >  5:332837/1‑147 (MQ=255)                                                                             tGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAAGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGa                                                                      >  5:338367/1‑147 (MQ=255)                                                                               gAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGatat                                                                   <  7:125813/148‑1 (MQ=255)                                                                                  ccAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAAGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCtt                                                                <  6:338367/148‑1 (MQ=255)                                                                                  ccAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAAGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCtt                                                                <  6:332837/148‑1 (MQ=255)                                                                                   cAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGTGGGGCGAGACGGTCAACGGTGATATCttt                                                               <  3:3699/148‑1 (MQ=255)                                                                                   cAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCttt                                                               <  1:144466/148‑1 (MQ=255)                                                                                               tctctGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCgg                                                   >  5:182083/1‑148 (MQ=255)                                                                                                 tctGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGa                                                 >  8:147443/1‑148 (MQ=255)                                                                                                  ctGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGc                                                                                                                                  <  6:44690/66‑1 (MQ=255)                                                                                                  ctGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGc                                                                                                                                  >  5:44690/1‑66 (MQ=255)                                                                                                  ctGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGAt                                                >  2:373016/1‑148 (MQ=255)                                                                                                   tGCGCCAATGCCAAACAATTTGCTCTCAAATGCGTCTATATCGTGCGCAGGGTTTAGATTCTCGCTGATAACCTCATCAAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGCTGATATCTTTCATAACATCCGGGAt                                                >  5:147787/1‑147 (MQ=255)                                                                                                      gcCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAg                                            <  1:17374/148‑1 (MQ=255)                                                                                                         aaTGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGa                                                 <  8:327218/140‑1 (MQ=255)                                                                                                         aaTGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGa                                                 >  7:327218/1‑140 (MQ=255)                                                                                                          aTGCCAAACAATTTGCTCTCAACGGCGGCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGt                                        <  5:106307/148‑1 (MQ=255)                                                                                                             ccAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTtca                                     <  7:230662/148‑1 (MQ=255)                                                                                                              cAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTtcat                                    >  6:2205/1‑148 (MQ=255)                                                                                                               aaaCAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGCTCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTtcat                                    >  2:338067/1‑147 (MQ=255)                                                                                                               aaaCAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGAt                                                >  2:203385/1‑135 (MQ=255)                                                                                                               aaaCAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGAt                                                <  1:203385/135‑1 (MQ=255)                                                                                                               aaaCAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTtcatc                                   <  1:373016/148‑1 (MQ=255)                                                                                                                   aaTTTGCTCTCAACTGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGTGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCg                               <  6:291214/148‑1 (MQ=255)                                                                                                                        gCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCGATgg                           <  2:4570/147‑1 (MQ=255)                                                                                                                         ctctCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCGATGgt                          <  5:2205/147‑1 (MQ=255)                                                                                                                               aCGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCGATGGTGTGcc                     >  2:125026/1‑146 (MQ=255)                                                                                                                                      cGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCGATGGTGTGCCGAGTTGCTc            >  3:309035/1‑148 (MQ=255)                                                                                                                                             gTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCGATGGTGTGCCGAGTTGCTCGCCaaa      <  7:147443/147‑1 (MQ=255)                                                                                                                                             gTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCGATGGTGTGCCGAGTTGCTCGCCaaa      <  4:309035/147‑1 (MQ=255)                                                                                                                                             gTGCGAAGGGTTTAGCTTGTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCTTGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCGATGGTGTGCCGAGTTGCTCGCCaaa      <  6:147787/147‑1 (MQ=255)                                                                                                                                              tGCGCAGGGTTTAGCTTCTCGCTGTTCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCCGGGCGAGACGGTGAACGGTGATATCTTTCATACCATCCGGGATGGAGATGTTCATCAGCGATGGTGTGCAGAGTTGCTCGCCcaaa     >  6:344512/1‑147 (MQ=255)                                                                                                                                                cgcAGGGTTTAGCTTCTCGCTGATCACCTCTTCCAGTGCTGCCAGCAGACGCTGACGCGGTGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCGATGGTGTGCCGAGTTGCTCGCCAAAATAg  <  1:79627/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                TCAGTCGGGTGGAAGTGCCCGGCGTCCAGGCACAGCGCAGTCTGGCGGCTGGTGGCATACCCCATGTAAAACTCATTGGAGCCAACCGTGTAGCTCTCTGCGCCAATGCCAAACAATTTGCTCTCAACGGCGTCGATATGGTGCGCAGGGTTTAGCTTCTCGCTGATCACCTCATCCAGTGCTGCCAGCAGACGCTGACGCGGGGCGAGACGGTCAACGGTGATATCTTTCATACCATCCGGGATCCAGATGTTCATCACCGATGGTGTGCCGAGTTGCTCGCCAAAATAG  >  NC_000913/4095183‑4095473

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑