Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,947,435	A→T	6.1%	S239S (TCT→TCA)	yebC ←	UPF0082 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,947,435	0	A	T	6.1%	206.6 / 3.8	65	S239S (TCT→TCA)	yebC	UPF0082 family protein
Reads supporting (aligned to +/- strand):  ref base A (27/34);  new base T (2/2);  total (29/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 = 4.83e-01

CCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACTTCCGCGCTGTCTGCTTTC  >  NC_000913/1947289‑1947579                                                                                                                                                   |                                                                                                                                                 ccAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCag                                                                                                                                                 <  8:47171/148‑1 (MQ=255)      tAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGAtt                                                                                                                                            >  1:370659/1‑148 (MQ=255)            aGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCg                                                                                                                                      <  8:164940/148‑1 (MQ=255)                 cctacctaCCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATg                                                                                                                                 >  3:120776/1‑148 (MQ=255)                  ctacctacCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATgg                                                                                                                                >  5:362774/1‑148 (MQ=255)                    acctacctGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGt                                                                                                                               >  5:386777/1‑147 (MQ=255)                     cctacctGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTaa                                                                                                                             <  3:337372/148‑1 (MQ=255)                         cctGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAAc                                                                                                                           <  7:137447/146‑1 (MQ=255)                         cctGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACtt                                                                                                                         >  3:290301/1‑148 (MQ=255)                           tGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGt                                                                                                                               >  7:37031/1‑140 (MQ=255)                           tGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGt                                                                                                                               <  8:37031/140‑1 (MQ=255)                           tGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTcc                                                                                                                       >  8:371061/1‑148 (MQ=255)                           tGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGACATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTcc                                                                                                                       >  6:153320/1‑148 (MQ=255)                             gCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCt                                                                                                                      <  3:323242/147‑1 (MQ=255)                              cGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCGGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTg                                                                                                                     >  6:247178/1‑147 (MQ=255)                              cGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGc                                                                                                                    <  7:70777/148‑1 (MQ=255)                                     cGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCAcgtcgt                                                                                                             <  2:137831/148‑1 (MQ=255)                                         gTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCAcgt                                                                                                                >  7:273313/1‑141 (MQ=255)                                         gTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCAcgt                                                                                                                <  8:273313/141‑1 (MQ=255)                                              cGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCa                                                                                                         <  1:276461/143‑1 (MQ=255)                                              cGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCa                                                                                                         >  2:276461/1‑143 (MQ=255)                                                 tCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTc                                                                                                   >  2:231463/1‑146 (MQ=255)                                                 tCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGGCTCCGTTTTTCTGTTTAGCAGGCCTCATCAGATAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCa                                                                                                 <  7:314071/148‑1 (MQ=255)                                                      cgcgAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCata                                                                                            >  8:130610/1‑148 (MQ=255)                                                          aaCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGa                                                                                        <  7:130610/148‑1 (MQ=255)                                                           aCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGAt                                                                                       >  8:175724/1‑148 (MQ=255)                                                             ccGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCAAGTCGTCGCAATCTTCCAGCATATCGAt                                                                                       >  8:76149/1‑146 (MQ=255)                                                                gATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGa                                                                                        <  7:284372/142‑1 (MQ=255)                                                                gATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGa                                                                                        >  8:284372/1‑142 (MQ=255)                                                                           agaATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTc                                                                       <  7:175724/148‑1 (MQ=255)                                                                               taataGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTg                                                                   <  7:371061/148‑1 (MQ=255)                                                                                        tCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTc                                                                                                   <  4:326735/107‑1 (MQ=255)                                                                                        tCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTc                                                                                                   >  3:326735/1‑107 (MQ=255)                                                                                        tCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTg                                                          <  6:362774/148‑1 (MQ=255)                                                                                        tCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTg                                                          >  8:329008/1‑148 (MQ=255)                                                                                          aCGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCa                                                        <  5:26041/148‑1 (MQ=255)                                                                                            gcgTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATc                                                      <  7:329008/148‑1 (MQ=255)                                                                                             cgTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGc                                                         >  8:214541/1‑144 (MQ=255)                                                                                             cgTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGc                                                         <  7:214541/144‑1 (MQ=255)                                                                                                   cGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCa                                               >  8:310878/1‑148 (MQ=255)                                                                                                    gTTTTTCTGGTTAGCAGGCCTCATCAGAGAGGAGCGGCGACCTCATCTGAGAGTGCACCGTTATGGAAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTATCGGTGCGGATTCTGCTTCCa                                                    <  6:369589/142‑1 (MQ=255)                                                                                                    gTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCa                                                                                     >  4:92994/1‑109 (MQ=255)                                                                                                    gTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCa                                                                                     <  3:92994/109‑1 (MQ=255)                                                                                                    gTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAg                                              <  7:310384/148‑1 (MQ=255)                                                                                                     ttttGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCTGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCa                                                    >  5:369588/1‑141 (MQ=255)                                                                                                     ttttGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCTGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCa                                                    <  6:369588/141‑1 (MQ=255)                                                                                                     ttttGCGGTTTAGCAGGCCTCATCAGAGAGTAGCTGCGACCTCATCTGTGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCCGACGCATAATTTTCGGTGCGGTTTCTGCATCCattcgatc                                             >  5:369589/1‑142 (MQ=255)                                                                                                      tttGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCt                                            >  1:231869/1‑148 (MQ=255)                                                                                                         gCTGTTTAGCAGGCCTCATCGGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCttt                                          <  5:247178/147‑1 (MQ=255)                                                                                                          cTGTTTAGCAGTCCTCATCAGAGAGTCGCGGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTg                                         <  3:279035/147‑1 (MQ=255)                                                                                                          cTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTgg                                        <  4:120776/148‑1 (MQ=255)                                                                                                          cTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTg                                         <  7:76149/147‑1 (MQ=255)                                                                                                          cTGTCTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTg                                         <  1:126809/147‑1 (MQ=255)                                                                                                             tttAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAg                                     <  2:231869/148‑1 (MQ=255)                                                                                                               tAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGa                                    <  4:290301/147‑1 (MQ=255)                                                                                                                aGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGAc                                   <  5:153320/147‑1 (MQ=255)                                                                                                                  cAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTTACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCGGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACgg                                 >  8:412582/1‑147 (MQ=255)                                                                                                                         catcaGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATg                          >  1:401570/1‑147 (MQ=255)                                                                                                                         catcaGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATg                          >  1:156467/1‑147 (MQ=255)                                                                                                                         catcaGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGAATCCATATCAGCTTTGGTAGACGGGATCATg                          >  2:104355/1‑147 (MQ=255)                                                                                                                               agagTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt                   <  1:104355/148‑1 (MQ=255)                                                                                                                                      cTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTGCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACTTCCGCGCt            <  1:231463/148‑1 (MQ=255)                                                                                                                                       tGCGACCTCATCAGAGATTTCACCGTTATGGTAAACtt                                                                                                                         <  2:341705/38‑1 (MQ=255)                                                                                                                                       tGCGACCTCATCAGAGATTTCACCGTTATGGTAAACtt                                                                                                                         >  1:341705/1‑38 (MQ=255)                                                                                                                                         cGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACTTCCGCGCTg           <  2:370659/146‑1 (MQ=255)                                                                                                                                             cTCATCAGAGATTTCACCGGTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACTTCCGCGCTGTCTGCt     <  3:236583/148‑1 (MQ=255)                                                                                                                                              tcatcaGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTACATCCATATCAGCTTTGGTAGACGGGATCATGGAAACTTCCGCGCTGTCTGCtt    >  7:234524/1‑148 (MQ=255)                                                                                                                                              tcatcaGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTACATCCATATCAGCTTTGGTAGACGGGATCATGGAAACTTCCGCGCTGTCTGCtt    >  5:23957/1‑148 (MQ=255)                                                                                                                                                 tcaGAGATTTAACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATAGATCAGACGCATCAGTTTCGGTGAGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACTTACGCGCTGTCTGCttt   >  7:203180/1‑146 (MQ=255)                                                                                                                                                  caGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTc                                                            >  7:24245/1‑88 (MQ=255)                                                                                                                                                  caGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTc                                                            <  8:24245/88‑1 (MQ=255)                                                                                                                                                  caGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACggg                                <  3:135738/116‑1 (MQ=255)                                                                                                                                                  caGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACggg                                >  4:135738/1‑116 (MQ=255)                                                                                                                                                  caGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCa                            <  1:81582/120‑1 (MQ=255)                                                                                                                                                  caGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCa                            >  2:81582/1‑120 (MQ=255)                                                                                                                                                  caGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACTTCCGCGCTGTCTGCttt   >  3:195189/1‑145 (MQ=255)                                                                                                                                                   agagATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACTTCCGCGCTGTCTGCTTTc  <  4:195189/145‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 CCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACTTCCGCGCTGTCTGCTTTC  >  NC_000913/1947289‑1947579

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