BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,505,916	T→C	M105V (ATG→GTG)	yhfV ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,505,916	0	T	C	86.1%	809.8 / 116.5	339	M105V (ATG→GTG)	yhfV	putative hydrolase
Reads supporting (aligned to +/- strand): ref base T (25/22); new base C (154/138); total (179/160)
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.40e-01

Only 100 of 360 total aligned reads displayed.

GGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGA  >  NC_000913/3505825‑3506005

 gcgTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCAcct                                                                                         >  2:3723184/1‑93 (MQ=255)

  cgTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCAcctc                                                                                        <  1:6449005/93‑1 (MQ=255)

   gTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTcc                                                                                       >  1:5164729/1‑93 (MQ=255)

   gTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTcc                                                                                       >  2:132376/1‑93 (MQ=255)

   gTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTcc                                                                                       >  2:2189113/1‑93 (MQ=255)

    tAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCAcctcct                                                                                      >  1:1784393/1‑93 (MQ=255)

    tAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCAcctcct                                                                                      >  2:5151671/1‑93 (MQ=255)

    tAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCAcctcct                                                                                      >  2:2885635/1‑93 (MQ=255)

     aaTCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTg                                                                                     >  2:3192795/1‑93 (MQ=255)

          ttCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCa                                                                                >  1:809291/1‑93 (MQ=255)

              ttCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTc                                                                            <  1:1072638/93‑1 (MQ=255)

               tCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTcc                                                                           >  2:2480913/1‑93 (MQ=255)

               tCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTcc                                                                           >  1:4228070/1‑93 (MQ=255)

                 gCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTg                                                                         <  2:4136255/93‑1 (MQ=255)

                 gCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTg                                                                         <  2:6675231/93‑1 (MQ=255)

                  cTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGc                                                                        >  2:5979875/1‑93 (MQ=255)

                  cTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGc                                                                        >  1:3907723/1‑93 (MQ=255)

                   tGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCa                                                                       >  1:2805456/1‑93 (MQ=255)

                   tGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCa                                                                       >  2:7128891/1‑93 (MQ=255)

                   tGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCa                                                                       <  1:4942543/93‑1 (MQ=255)

                    ggTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCAc                                                                      <  2:5084015/93‑1 (MQ=255)

                    ggTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCAc                                                                      <  1:2188608/93‑1 (MQ=255)

                      tGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGc                                                                    >  1:2456585/1‑93 (MQ=255)

                       gCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCt                                                                   >  2:2076513/1‑93 (MQ=255)

                        ccGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTg                                                                  <  2:942383/93‑1 (MQ=255)

                         cGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGc                                                                 <  1:6965902/93‑1 (MQ=255)

                          gATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGTTACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCg                                                                <  2:4643375/93‑1 (MQ=255)

                          gATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCg                                                                >  1:2671638/1‑93 (MQ=255)

                           aTCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCgg                                                               <  2:241191/93‑1 (MQ=255)

                           aTCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCgg                                                               >  2:65701/1‑93 (MQ=255)

                            tctcCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCggg                                                              >  2:4977013/1‑93 (MQ=255)

                            tctcCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCggg                                                              >  2:4725687/1‑93 (MQ=255)

                            tctcCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCggg                                                              <  1:657905/93‑1 (MQ=255)

                             ctcCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGt                                                             <  2:6754100/93‑1 (MQ=255)

                             ctcCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGt                                                             <  1:7219805/93‑1 (MQ=255)

                                cgcgATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGc                                                          >  1:5640951/1‑93 (MQ=255)

                                 gcgATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGcc                                                         >  2:2207570/1‑93 (MQ=255)

                                 gcgATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGcc                                                         >  2:2676298/1‑93 (MQ=255)

                                   gatgatCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCcac                                                       <  2:3834567/93‑1 (MQ=255)

                                     tgatCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACAt                                                     >  1:6389832/1‑93 (MQ=255)

                                       atCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGt                                                   >  2:6164272/1‑93 (MQ=255)

                                       atCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGt                                                   >  2:2513147/1‑93 (MQ=255)

                                        tCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGtt                                                  <  1:101016/93‑1 (MQ=255)

                                         cccGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTc                                                 <  1:6766607/93‑1 (MQ=255)

                                         cccGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTc                                                 <  2:4856020/93‑1 (MQ=255)

                                          ccGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTcc                                                <  1:4281590/93‑1 (MQ=255)

                                            ggCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCgg                                              <  1:308629/93‑1 (MQ=255)

                                            ggCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCACCACATGTTCCgg                                              >  1:3182667/1‑93 (MQ=255)

                                             gCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCggg                                             >  2:5752332/1‑93 (MQ=255)

                                                ttCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGaaa                                          >  1:7183757/1‑93 (MQ=255)

                                                 tCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGaaaa                                         <  1:4848267/93‑1 (MQ=255)

                                                  cAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGaaaaa                                        <  2:5398684/93‑1 (MQ=255)

                                                  cAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGaaaaa                                        <  1:7358057/93‑1 (MQ=255)

                                                  cAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGaaaaa                                        <  2:2808878/93‑1 (MQ=255)

                                                   aGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAAc                                       >  2:799257/1‑93 (MQ=255)

                                                   aGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAAc                                       <  2:3442013/93‑1 (MQ=255)

                                                    gCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAAcg                                      <  2:4355429/93‑1 (MQ=255)

                                                     cTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAAcgc                                     <  1:5908290/93‑1 (MQ=255)

                                                      tCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAAcgcg                                    <  2:889812/93‑1 (MQ=255)

                                                      tCCGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAAcgcg                                    >  2:6064742/1‑93 (MQ=255)

                                                       ccGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGt                                   <  1:906653/93‑1 (MQ=255)

                                                        cGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTc                                  <  1:4794383/93‑1 (MQ=255)

                                                        cGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTc                                  <  1:1947900/93‑1 (MQ=255)

                                                        cGTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTc                                  <  1:418790/93‑1 (MQ=255)

                                                         gTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTcc                                 <  2:2581158/93‑1 (MQ=255)

                                                         gTGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTcc                                 <  1:1614529/93‑1 (MQ=255)

                                                          tGCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCt                                >  2:6593141/1‑93 (MQ=255)

                                                           gCCATCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTg                               <  2:6860368/93‑1 (MQ=255)

                                                             cATCGATACCCTGTTCAATTTCATCGACCCCCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGt                             <  1:4468298/93‑1 (MQ=255)

                                                               tCGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGtaa                           >  2:3565657/1‑93 (MQ=255)

                                                                cGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGtaat                          <  1:55224/93‑1 (MQ=255)

                                                                cGATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGtaat                          <  1:3511532/93‑1 (MQ=255)

                                                                 gATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGtaata                         <  2:631405/93‑1 (MQ=255)

                                                                 gATACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGtaata                         >  2:3987963/1‑93 (MQ=255)

                                                                  aTACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGtaataa                        <  2:1953381/93‑1 (MQ=255)

                                                                   tACCCTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAAc                       <  2:7323707/93‑1 (MQ=255)

                                                                     cccTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCg                     <  1:5313996/93‑1 (MQ=255)

                                                                     cccTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCg                     >  2:1767634/1‑93 (MQ=255)

                                                                     cccTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCg                     <  2:1508949/93‑1 (MQ=255)

                                                                      ccTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCgg                    >  1:3335574/1‑93 (MQ=255)

                                                                       cTGTTCAATTTCATCGACCCCCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGt                   <  2:1142267/93‑1 (MQ=255)

                                                                       cTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGt                   >  1:4411533/1‑93 (MQ=255)

                                                                       cTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGt                   >  1:6027839/1‑93 (MQ=255)

                                                                       cTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGt                   <  1:2982920/93‑1 (MQ=255)

                                                                       cTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGt                   >  1:3551546/1‑93 (MQ=255)

                                                                       cTGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGt                   >  1:6112760/1‑93 (MQ=255)

                                                                        tGTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTa                  <  1:2331139/93‑1 (MQ=255)

                                                                         gTTCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTAc                 >  1:1220459/1‑93 (MQ=255)

                                                                           tCAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAg               >  2:558010/1‑93 (MQ=255)

                                                                            cAATTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAgg              >  2:5449241/1‑93 (MQ=255)

                                                                             aaTTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGc             <  1:1420368/93‑1 (MQ=255)

                                                                             aaTTTCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGc             <  1:3360353/93‑1 (MQ=255)

                                                                                 tCATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGccacc         >  2:4511606/1‑93 (MQ=255)

                                                                                  cATCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGccacca        >  2:1074219/1‑93 (MQ=255)

                                                                                   aTCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCcaccac       >  1:869884/1‑93 (MQ=255)

                                                                                   aTCGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCcaccac       >  1:6248033/1‑93 (MQ=255)

                                                                                     cGACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGt     <  2:5911855/93‑1 (MQ=255)

                                                                                      gACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACATGCCACCACGtt    >  2:6125/1‑93 (MQ=255)

                                                                                      gACCACCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGtt    <  2:1577374/93‑1 (MQ=255)

                                                                                       accaccTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTg   >  2:5708665/1‑93 (MQ=255)

GGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGA  >  NC_000913/3505825‑3506005

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