BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	178,128	A→G	I112T (ATT→ACT)	yadT ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	178,128	0	A	G	99.4%	536.7 / ‑5.7	163	I112T (ATT→ACT)	yadT	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (72/90); total (72/91)
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.98e-01

Only 100 of 167 total aligned reads displayed.

TGATCCAGCAGGGATTGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAAT  >  NC_000913/178038‑178216

tGATCCAGCAGGGATTGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTg                                                                                        >  1:949776/1‑93 (MQ=255)

 gATCCAGCAGGGATTGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGc                                                                                       >  1:2196833/1‑93 (MQ=255)

  aTCCAGCAGGGATTGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCt                                                                                      <  2:3234690/93‑1 (MQ=255)

   tCCAGCAGGGATTGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCtt                                                                                     <  2:4025544/93‑1 (MQ=255)

    ccAGCAGGGATTGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTg                                                                                    <  2:597386/93‑1 (MQ=255)

     cagcagGGATTGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGt                                                                                   >  2:739837/1‑93 (MQ=255)

     cagcagGGATTGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGt                                                                                   >  1:1732201/1‑93 (MQ=255)

      agcagGGATTGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTc                                                                                  <  1:3022154/93‑1 (MQ=255)

      agcagGGATTGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTc                                                                                  <  1:2947608/93‑1 (MQ=255)

        cagGGATTGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGc                                                                                <  2:135416/93‑1 (MQ=255)

         agGGATTGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCa                                                                               >  1:2298565/1‑93 (MQ=255)

         agGGATTGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCa                                                                               >  1:4033207/1‑93 (MQ=255)

              ttGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGAc                                                                          <  1:3747181/93‑1 (MQ=255)

              ttGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGAc                                                                          >  1:1182877/1‑93 (MQ=255)

                gcgcGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGAccc                                                                        <  2:1254220/93‑1 (MQ=255)

                gcgcGGCTTGTTCGGACTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAAAGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGAccc                                                                        >  2:3734811/1‑93 (MQ=255)

                   cgGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCcaca                                                                     >  2:1101965/1‑93 (MQ=255)

                   cgGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCcaca                                                                     >  1:4183699/1‑93 (MQ=255)

                   cgGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCcaca                                                                     <  2:4472422/93‑1 (MQ=255)

                   cgGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCcaca                                                                     <  1:1787766/93‑1 (MQ=255)

                     gCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATc                                                                   >  2:3938039/1‑93 (MQ=255)

                     gCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATc                                                                   >  2:104130/1‑93 (MQ=255)

                        tGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACt                                                                >  1:2037673/1‑93 (MQ=255)

                         gTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACtt                                                               >  1:2466171/1‑93 (MQ=255)

                         gTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACtt                                                               <  2:1615891/93‑1 (MQ=255)

                          ttCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACttt                                                              >  1:994595/1‑93 (MQ=255)

                             ggCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTAt                                                           >  2:4174091/1‑93 (MQ=255)

                              gCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTAtt                                                          >  1:812864/1‑93 (MQ=255)

                              gCCTTGTCTGGTGGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTAtt                                                          <  2:2244290/93‑1 (MQ=255)

                                  tGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCa                                                      <  1:4623357/93‑1 (MQ=255)

                                   ttCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAg                                                     <  1:1687052/92‑1 (MQ=255)

                                   gTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTTGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAg                                                     <  2:2183243/93‑1 (MQ=255)

                                      tGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGa                                                  <  1:552515/93‑1 (MQ=255)

                                      tGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGa                                                  <  2:3802086/93‑1 (MQ=255)

                                      tGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGa                                                  >  1:3153851/1‑93 (MQ=255)

                                       ggTTGCGGGCTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGaa                                                 <  1:1502894/93‑1 (MQ=255)

                                       ggTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGaa                                                 <  1:4136474/93‑1 (MQ=255)

                                        gTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAg                                                >  2:2036846/1‑93 (MQ=255)

                                             ggACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGGAGCCAGc                                           <  2:476258/93‑1 (MQ=255)

                                             ggACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGc                                           >  2:1520982/1‑93 (MQ=255)

                                             ggACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGc                                           <  1:2611713/93‑1 (MQ=255)

                                             cgACTCCCAGTGGCCCGGTGACGTAACCGCTTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGc                                           <  1:1273899/92‑1 (MQ=255)

                                              gACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGGCCCACATCACTTTTATTCCCAGCGAAGCCAGCt                                          <  1:629990/93‑1 (MQ=255)

                                              gACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCt                                          <  2:3553317/93‑1 (MQ=255)

                                              gACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCt                                          >  2:59618/1‑93 (MQ=255)

                                               aCTCCTGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTg                                         <  2:3947899/93‑1 (MQ=255)

                                               aCTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTg                                         >  2:3033574/1‑93 (MQ=255)

                                                cTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTgg                                        >  2:4344613/1‑93 (MQ=255)

                                                   cAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCa                                     >  2:3072450/1‑93 (MQ=255)

                                                    aGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCaa                                    >  1:226024/1‑93 (MQ=255)

                                                     gggggCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAAc                                   <  2:2574640/93‑1 (MQ=255)

                                                     gggggCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAAc                                   >  1:951196/1‑93 (MQ=255)

                                                     gggggCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAAc                                   <  2:195545/93‑1 (MQ=255)

                                                      ggggCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAAcc                                  <  2:2395513/93‑1 (MQ=255)

                                                      ggggCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAAcc                                  <  2:3522854/93‑1 (MQ=255)

                                                      ggggCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAAcc                                  <  2:1817235/93‑1 (MQ=255)

                                                         gCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGc                               <  1:3841204/93‑1 (MQ=255)

                                                         gCCAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGc                               <  1:2542159/93‑1 (MQ=255)

                                                           cAGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCg                             >  2:969558/1‑93 (MQ=255)

                                                            aGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGc                            <  1:143352/93‑1 (MQ=255)

                                                            aGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGc                            >  2:2850426/1‑93 (MQ=255)

                                                            aGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGc                            >  1:520610/1‑93 (MQ=255)

                                                            aGTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGc                            <  2:2609887/93‑1 (MQ=255)

                                                             gTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGcccat                           >  1:3692138/1‑90 (MQ=255)

                                                             gTTGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCt                           >  1:1663663/1‑93 (MQ=255)

                                                              ttGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTc                          <  1:2162022/93‑1 (MQ=255)

                                                              ttGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTc                          >  1:4157027/1‑93 (MQ=255)

                                                              ttGACGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTc                          >  2:4354757/1‑93 (MQ=255)

                                                                  cGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCa                      >  1:3253848/1‑93 (MQ=255)

                                                                  cGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCa                      <  1:1917992/93‑1 (MQ=255)

                                                                  cGTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCa                      >  2:872570/1‑93 (MQ=255)

                                                                   gTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCAt                     <  1:772207/93‑1 (MQ=255)

                                                                   gTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCAt                     >  2:2948541/1‑93 (MQ=255)

                                                                   gTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCAt                     >  2:765146/1‑93 (MQ=255)

                                                                   gTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCAt                     <  2:1078923/93‑1 (MQ=255)

                                                                   gTAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCAt                     >  1:3324535/1‑93 (MQ=255)

                                                                    tAACGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCAtt                    <  1:3384221/93‑1 (MQ=255)

                                                                    tAACCCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCAtt                    <  2:3386142/93‑1 (MQ=255)

                                                                     aaCGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTa                   <  2:2973976/93‑1 (MQ=255)

                                                                      aCGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTAc                  >  2:841623/1‑93 (MQ=255)

                                                                      aCGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTAc                  <  2:3835337/93‑1 (MQ=255)

                                                                      aCGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTAc                  >  1:123456/1‑93 (MQ=255)

                                                                       cGCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTAcc                 <  2:4405356/93‑1 (MQ=255)

                                                                        gCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCt                <  2:1884802/93‑1 (MQ=255)

                                                                        gCATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCt                <  1:965057/93‑1 (MQ=255)

                                                                         cATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTc               <  2:1765822/93‑1 (MQ=255)

                                                                         cATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTc               <  1:2670178/93‑1 (MQ=255)

                                                                         cATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTc               >  2:2000498/1‑93 (MQ=255)

                                                                         cATTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTc               <  1:4415544/93‑1 (MQ=255)

                                                                          aTTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTcc              <  2:4525106/93‑1 (MQ=255)

                                                                          aTTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTcc              <  2:1853724/93‑1 (MQ=255)

                                                                          aTTGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTcc              >  2:3210962/1‑93 (MQ=255)

                                                                           ttGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTCCa             <  1:1765942/93‑1 (MQ=255)

                                                                           ttGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTCCa             >  1:1952431/1‑93 (MQ=255)

                                                                            tGGCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTCCAc            <  1:1437310/93‑1 (MQ=255)

                                                                             ggCAATTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTCCACg           >  1:4481745/1‑93 (MQ=255)

                                                                                aaTTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTCCACGCCa        <  1:402436/93‑1 (MQ=255)

                                                                                aaTTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTCCACGCCa        <  1:4232717/93‑1 (MQ=255)

                                                                                 aTTTGTTCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAg       <  1:3132016/93‑1 (MQ=255)

                                                                                      ttCAGTGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAAt  <  1:2492314/93‑1 (MQ=255)

TGATCCAGCAGGGATTGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAAT  >  NC_000913/178038‑178216

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