BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	788,889	A→G	L105P (CTG→CCG)	galK ←	galactokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	788,889	0	A	G	78.1%	580.7 / 196.4	320	L105P (CTG→CCG)	galK	galactokinase
Reads supporting (aligned to +/- strand): ref base A (37/33); new base G (123/127); total (160/160)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

Only 100 of 321 total aligned reads displayed.

GACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGA  >  NC_000913/788798‑788979

gACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAg                                                                                           >  2:3465312/1‑93 (MQ=255)

 aCTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGt                                                                                          <  2:5915514/93‑1 (MQ=255)

 aCTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGt                                                                                          <  2:6892522/93‑1 (MQ=255)

 aCTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGt                                                                                          >  1:4941993/1‑93 (MQ=255)

  cTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGtt                                                                                         <  2:1117817/93‑1 (MQ=255)

  cTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGtt                                                                                         <  1:5754839/93‑1 (MQ=255)

     ccAGTGAAGCTGAAGAACATAGAGCGGCCGCCTGCGGCACATTGGCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCa                                                                                      <  1:4794149/93‑1 (MQ=255)

      cAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAg                                                                                     <  2:2810107/93‑1 (MQ=255)

       aGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGa                                                                                    <  1:3142345/93‑1 (MQ=255)

         tGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATg                                                                                  >  1:1118048/1‑93 (MQ=255)

           aaGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGtt                                                                                <  2:4319732/93‑1 (MQ=255)

            aGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGttt                                                                               <  1:3497580/93‑1 (MQ=255)

               ggAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTcac                                                                            >  2:5528653/1‑93 (MQ=255)

               ggAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTcac                                                                            <  2:5340929/93‑1 (MQ=255)

                gaagaaCTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTcacc                                                                           >  2:2451769/1‑93 (MQ=255)

                gaagaaCTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTcacc                                                                           >  2:5665235/1‑93 (MQ=255)

                gaagaaCTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTcacc                                                                           >  1:3698777/1‑93 (MQ=255)

                 aagaaCTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTcacca                                                                          <  1:7479583/93‑1 (MQ=255)

                 aagaaCTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTcacca                                                                          <  2:4962973/93‑1 (MQ=255)

                  agaaCTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTcaccac                                                                         >  2:6013602/1‑93 (MQ=255)

                   gaaCTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacg                                                                        >  2:2120373/1‑93 (MQ=255)

                    aaCTTAACCCGGCTCCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCAccacgc                                                                       <  1:3268138/93‑1 (MQ=255)

                     aCTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCAccacgcc                                                                      <  2:2861031/93‑1 (MQ=255)

                      cTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacgcca                                                                     >  2:1047501/1‑93 (MQ=255)

                       ttAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCAccacgccac                                                                    <  2:7524197/93‑1 (MQ=255)

                       ttAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCAccacgccac                                                                    <  1:71877/93‑1 (MQ=255)

                         aaCCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGa                                                                  <  1:6259030/93‑1 (MQ=255)

                         aaCCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGa                                                                  <  1:1592264/93‑1 (MQ=255)

                          aCCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGaa                                                                 <  2:5213916/93‑1 (MQ=255)

                            ccGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACg                                                               <  1:235482/93‑1 (MQ=255)

                            ccGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACg                                                               <  1:6123665/93‑1 (MQ=255)

                               gCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAg                                                            >  2:3308532/1‑93 (MQ=255)

                               gCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAg                                                            <  2:6402804/93‑1 (MQ=255)

                               gCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAg                                                            <  2:3996577/93‑1 (MQ=255)

                                cACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGt                                                           <  2:6523054/93‑1 (MQ=255)

                                 aCCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGtt                                                          <  1:7368650/93‑1 (MQ=255)

                                     tGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGcc                                                      >  1:6073423/1‑93 (MQ=255)

                                       cGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCa                                                    <  1:7336925/93‑1 (MQ=255)

                                        ggCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCAt                                                   >  2:134289/1‑93 (MQ=255)

                                        ggCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCAt                                                   <  2:611625/93‑1 (MQ=255)

                                         gCACATTGCCGCTGTTCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCAtt                                                  <  1:6834031/93‑1 (MQ=255)

                                         gCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCAtt                                                  >  2:2274366/1‑93 (MQ=255)

                                           acaTTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGa                                                >  2:2521630/1‑93 (MQ=255)

                                           acaTTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGa                                                >  2:1754356/1‑93 (MQ=255)

                                           acaTTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGa                                                >  1:3426068/1‑93 (MQ=255)

                                           acaTTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCAttnn                                                >  2:382629/1‑91 (MQ=255)

                                            caTTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGAt                                               >  1:5612515/1‑93 (MQ=255)

                                            caTTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGAt                                               >  1:3570174/1‑93 (MQ=255)

                                             aTTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATa                                              <  1:1719514/93‑1 (MQ=255)

                                              ttGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAg                                             >  1:6999274/1‑93 (MQ=255)

                                              ttGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAg                                             <  1:53004/93‑1 (MQ=255)

                                              ttGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAg                                             >  1:4337992/1‑93 (MQ=255)

                                              ttGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAg                                             <  1:4016293/93‑1 (MQ=255)

                                              ttGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAg                                             >  2:720158/1‑93 (MQ=255)

                                               tGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGt                                            <  2:393722/93‑1 (MQ=255)

                                               tGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGt                                            <  1:2272775/93‑1 (MQ=255)

                                                 ccGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGttt                                          <  1:6232937/93‑1 (MQ=255)

                                                 ccGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGttt                                          <  2:1641574/93‑1 (MQ=255)

                                                 ccGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGttt                                          >  1:1479585/1‑93 (MQ=255)

                                                    cTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCa                                       <  2:880813/93‑1 (MQ=255)

                                                    cTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCa                                       <  1:986251/93‑1 (MQ=255)

                                                    cTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCa                                       <  1:750693/93‑1 (MQ=255)

                                                      gATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCAtg                                     <  2:7414869/93‑1 (MQ=255)

                                                      gATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCAtg                                     <  1:3950902/93‑1 (MQ=255)

                                                        tCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCAtgtg                                   >  2:1567294/1‑93 (MQ=255)

                                                          aCCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCg                                 >  2:7092825/1‑93 (MQ=255)

                                                           ccATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGa                                >  1:2499986/1‑93 (MQ=255)

                                                            cATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGAc                               <  2:840935/93‑1 (MQ=255)

                                                            cATGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGAc                               <  1:4049126/93‑1 (MQ=255)

                                                            cATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGAc                               >  1:1889688/1‑93 (MQ=255)

                                                            cATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTAAGCCCATTGATAGTTTTCATGTGCGAc                               <  2:4832132/93‑1 (MQ=255)

                                                             aTGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACa                              >  1:466789/1‑93 (MQ=255)

                                                              tGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACaa                             <  1:5708208/93‑1 (MQ=255)

                                                              tGTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACaa                             >  1:2094749/1‑93 (MQ=255)

                                                               gTCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAAt                            >  1:2423391/1‑93 (MQ=255)

                                                                tCCACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATg                           <  1:279406/93‑1 (MQ=255)

                                                                  cACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATggg                         >  2:4204770/1‑93 (MQ=255)

                                                                  cACGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATggg                         >  2:2069561/1‑93 (MQ=255)

                                                                   aCGCCGCCGAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGgc                        <  1:290288/93‑1 (MQ=255)

                                                                   aCGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGgc                        >  2:1993340/1‑93 (MQ=255)

                                                                      ccgccgAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCa                     <  1:1867058/93‑1 (MQ=255)

                                                                      ccgccgAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCa                     <  2:3380658/93‑1 (MQ=255)

                                                                      ccgccgAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCa                     <  2:4636110/93‑1 (MQ=255)

                                                                        gccgAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATc                   <  1:1365427/93‑1 (MQ=255)

                                                                        gccgAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATc                   >  1:126760/1‑93 (MQ=255)

                                                                         ccgAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAGTGGGCGCATCg                  >  1:2215623/1‑93 (MQ=255)

                                                                          cgAAGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGa                 >  1:6327614/1‑93 (MQ=255)

                                                                            aaGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAgg               >  2:3141174/1‑93 (MQ=255)

                                                                            aaGCTGTTGTTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAgg               >  1:4549940/1‑93 (MQ=255)

                                                                             aGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAggg              >  2:2922409/1‑93 (MQ=255)

                                                                             aGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAggg              >  2:4488762/1‑93 (MQ=255)

                                                                                tgttgtTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGaaa           >  2:6102841/1‑93 (MQ=255)

                                                                                tgttgtTACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGaaa           >  2:4326967/1‑93 (MQ=255)

                                                                                 gttgttACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGaaaa          >  1:4131881/1‑93 (MQ=255)

                                                                                  ttgttACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAAc         >  1:1009004/1‑93 (MQ=255)

                                                                                      tACGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACtcgt     <  1:6741473/93‑1 (MQ=255)

                                                                                      tACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACtcgt     <  2:6227498/93‑1 (MQ=255)

                                                                                        cGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACtcgtcg   <  2:5064445/93‑1 (MQ=255)

                                                                                         gCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGa  <  1:1630109/93‑1 (MQ=255)

                                                                                         gCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGa  >  1:4614988/1‑93 (MQ=255)

GACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGA  >  NC_000913/788798‑788979

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