Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,501,912	T→A	100%	I117I (ATT→ATA)	tehB →	tellurite, selenium methyltransferase, SAM‑dependent; tellurite, selenium resistance protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,501,912	0	T	A	97.9%	305.1 / ‑1.7	93	I117I (ATT→ATA)	tehB	tellurite, selenium methyltransferase, SAM‑dependent; tellurite, selenium resistance protein
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base A (42/49);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

TTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGT  >  NC_000913/1501815‑1502006                                                                                                  |                                                                                               ttGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTCAAACCATCCCCGGGTTGATAGc                                                                                              <  2:1754818/100‑1 (MQ=255)    aTCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCaa                                                                                           >  1:237320/1‑100 (MQ=255)    aTCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACAATCCCCGGGTTGATAGCCaa                                                                                           >  1:34984/1‑100 (MQ=255)     tCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAat                                                                                          <  1:1000999/100‑1 (MQ=255)      cTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAata                                                                                         >  2:2446271/1‑100 (MQ=255)      cTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAata                                                                                         <  1:145546/100‑1 (MQ=255)       tGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAatat                                                                                        >  2:19655/1‑100 (MQ=255)       tGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAata                                                                                         <  1:2423982/99‑1 (MQ=255)        gAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATg                                                                                       >  1:1874251/1‑100 (MQ=255)        gAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTAAGGATAAAACCATCCCCGGGTTGATAGCCAatat                                                                                        >  1:1043761/1‑99 (MQ=255)         aaTAACCTCACATTTGATAGACAGTACGACTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGCCTAAAACCATCCCCGGGTTGATAGCCAATATGc                                                                                      <  2:1824255/100‑1 (MQ=255)             aCCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGAGTTGATAGCCAATATGCAAc                                                                                   >  1:1497172/1‑99 (MQ=255)               cTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGtt                                                                                >  2:1141242/1‑100 (MQ=255)               cTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGtt                                                                                >  2:794278/1‑100 (MQ=255)               cTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGtt                                                                                >  2:2332232/1‑100 (MQ=255)                tCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGttn                                                                               <  2:756329/100‑2 (MQ=255)                tCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGtt                                                                                >  2:2001130/1‑99 (MQ=255)                tCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTg                                                                               <  1:543116/100‑1 (MQ=255)                 cccaTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGc                                                                              >  2:2239409/3‑100 (MQ=255)                 cacaTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGTTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGc                                                                              <  1:25695/100‑1 (MQ=255)                 cacaTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGtt                                                                                >  2:66796/1‑98 (MQ=255)                 cacaTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGc                                                                              <  2:546084/100‑1 (MQ=255)                   caTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCAc                                                                            <  1:1472071/100‑1 (MQ=255)                   caTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCAc                                                                            <  2:2231508/100‑1 (MQ=255)                   caTTTGATAGACAGTACGATTTTATGCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCAc                                                                            <  2:365620/100‑1 (MQ=255)                    aTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCAc                                                                            >  2:708684/1‑99 (MQ=255)                       tGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTaaa                                                                        >  2:800285/1‑100 (MQ=255)                        gATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAAc                                                                       <  2:1178984/100‑1 (MQ=255)                        gATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAAc                                                                       <  2:1178873/100‑1 (MQ=255)                        gATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAAc                                                                       >  1:46962/1‑100 (MQ=255)                        gATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAAc                                                                       <  1:1252871/100‑1 (MQ=255)                        gATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAAc                                                                       >  2:1840061/1‑100 (MQ=255)                         aTAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATACCCGGGTTGATAGCCAATATGCAACGTTACACTAAAcc                                                                      >  2:1926894/1‑100 (MQ=255)                            gACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtgg                                                                   >  2:66988/1‑100 (MQ=255)                            gACAGTACGATTTTATTCTTTCGACTGTGGTGATGATGTTCCTTGAGGCTAAAACCATCCTCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtgg                                                                   >  2:2483415/1‑100 (MQ=255)                             aCAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtggt                                                                  <  2:3019555/100‑1 (MQ=255)                              cAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtggtg                                                                 <  2:1519873/100‑1 (MQ=255)                              cAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtggtg                                                                 <  1:350131/100‑1 (MQ=255)                              cAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtggtg                                                                 <  1:1541916/100‑1 (MQ=255)                                gTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTggtggt                                                               >  1:895608/1‑100 (MQ=255)                                  aCGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTa                                                             <  1:1652862/100‑1 (MQ=255)                                   cGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTAc                                                            >  1:2229923/1‑100 (MQ=255)                                     aTTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGTGTTTATAGCCAATATACAACGTTGCACTAAACCTGGTGGTTACaa                                                          >  1:2353012/1‑100 (MQ=255)                                      ttttATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAc                                                         <  1:1241527/100‑1 (MQ=255)                                          aTTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGa                                                     <  2:1615539/100‑1 (MQ=255)                                           ttCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGAt                                                    <  1:66988/100‑1 (MQ=255)                                           ttCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGAt                                                    >  1:1397376/1‑100 (MQ=255)                                           ttCTTTCGACTGTAGTGCTAATATTCCTTGAGACTAAAACCATGCCCGGGTTGATAGCCAATATGCAACGTTGAACTAAACCTGGTGGTTACAACCTGAt                                                    >  2:695155/1‑100 (MQ=255)                                            tCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGAt                                                    >  2:845135/1‑99 (MQ=255)                                             cTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATtg                                                  >  2:852488/1‑100 (MQ=255)                                              tttCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATtgt                                                 >  1:437441/1‑100 (MQ=255)                                                 cGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggc                                              <  1:1840083/100‑1 (MQ=255)                                                 cGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTATggc                                              >  1:150373/1‑100 (MQ=255)                                                  gACTGTGGTGCTTATTTTCCTTAATTATTAAACCTTCCTCTAGTTGATAGCCAATATGCTACGTTTTACGAATCCTTTTGGTTACAACCTGATTGTggc                                              <  2:1317624/99‑1 (MQ=255)                                                  gACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCAGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggcg                                             >  2:1188866/1‑100 (MQ=255)                                                     tgtgGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGgcggcg                                          <  1:1970579/100‑1 (MQ=255)                                                          tgCTGATGTTCCTTGAGGCTAAAACCATCCTCGGGTTGATAGGCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATg                                       <  1:2827134/98‑1 (MQ=255)                                                          tgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGa                                     <  2:2567036/100‑1 (MQ=255)                                                             tGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATAc                                  <  1:1899041/100‑1 (MQ=255)                                                             tGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATAc                                  >  2:820030/1‑100 (MQ=255)                                                              gctGTTCCTCGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGAAACGTTGTACTAAATCTTGTGGTTACTACCTTATTGTGGCGGCGATGGGTAcc                                 <  2:245370/98‑1 (MQ=255)                                                               aTGTTCCTTGAGGCTAAAACCAACCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATAcc                                 <  2:2988696/99‑1 (MQ=255)                                                                tGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGc                               <  2:1508222/100‑1 (MQ=255)                                                                tGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGc                               >  2:1135959/1‑100 (MQ=255)                                                                  ttacTTGAGTCTAAAACCATCCCCGGGTTGATAGCGAATATGCAACGTTGCACTAATCCTGGTGGTTACAACCTGATTGTGGCGGCGATGGAt                                    <  2:1497998/90‑1 (MQ=255)                                                                    ccTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGAt                           <  2:342124/100‑1 (MQ=255)                                                                    ccTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGAt                           <  1:374442/100‑1 (MQ=255)                                                                     cTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGAtt                          >  1:1509779/1‑100 (MQ=255)                                                                      ttGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTa                         <  1:1985339/100‑1 (MQ=255)                                                                      ttGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTa                         <  1:15225/100‑1 (MQ=255)                                                                        gAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTAGTAGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATc                       >  1:608879/1‑100 (MQ=255)                                                                         aGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATcc                      >  1:2771436/1‑100 (MQ=255)                                                                         aGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATcc                      <  2:1910996/100‑1 (MQ=255)                                                                           gCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCAt                    <  1:877956/100‑1 (MQ=255)                                                                            cTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCAt                    >  2:2078843/1‑99 (MQ=255)                                                                             tAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACATGATTGTGGCGGCGATGGATACCGCTGATTATCCATGt                  >  2:2023745/1‑100 (MQ=255)                                                                              aaaaCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTa                 >  1:973434/1‑100 (MQ=255)                                                                              aaaaCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTa                 >  1:318380/1‑100 (MQ=255)                                                                               aaaCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTa                 >  2:1545601/1‑99 (MQ=255)                                                                               aaaCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTAc                <  1:1198702/100‑1 (MQ=255)                                                                                 aCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCg              >  1:2404514/1‑100 (MQ=255)                                                                                   cATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTc            >  2:101358/1‑100 (MQ=255)                                                                                    aTCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCTGCGATGGATACCGTTGATTATCCATGTACCGTCg           <  2:283089/100‑1 (MQ=255)                                                                                    aTCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCg           <  2:1591933/100‑1 (MQ=255)                                                                                    aTCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCg           <  2:1679960/100‑1 (MQ=255)                                                                                     tCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGTTTGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCgg          <  1:125239/100‑1 (MQ=255)                                                                                     tCCCCGGGTTGATAGCCAATATGCAACATTGCACTAAACCTGGTGGTTACAACCTGATTGTGACGGCGATGGAAACCGCTGATTAACCATGTACCGTc            >  2:790861/1‑98 (MQ=255)                                                                                      ccccGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGc         >  1:2409121/1‑100 (MQ=255)                                                                                         caggTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTc      <  2:236962/98‑1 (MQ=255)                                                                                         cGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTc      <  1:527859/100‑1 (MQ=255)                                                                                           ggTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGc         <  1:559446/95‑1 (MQ=255)                                                                                            gTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTccc    <  2:2694464/99‑1 (MQ=255)                                                                                             ttGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGt  <  1:765880/100‑1 (MQ=255)                                                                                                  |                                                                                               TTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGT  >  NC_000913/1501815‑1502006

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