BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,754,991	A→G	T424A (ACT→GCT)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,754,991	0	A	G	83.9%	520.6 / 88.2	211	T424A (ACT→GCT)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand): ref base A (21/13); new base G (77/100); total (98/113)
Fisher's exact test for biased strand distribution p-value = 6.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

Only 100 of 214 total aligned reads displayed.

GCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGG  >  NC_000913/1754901‑1755082

gCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCt                                                                                           <  2:3398926/93‑1 (MQ=255)

 ccaccaTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCtt                                                                                          >  1:2159025/1‑93 (MQ=255)

  caccaTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTc                                                                                         <  1:1711867/93‑1 (MQ=255)

   accaTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCt                                                                                        <  2:4477818/93‑1 (MQ=255)

   accaTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCt                                                                                        >  1:6534312/1‑93 (MQ=255)

     caTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTAc                                                                                      <  2:5034584/93‑1 (MQ=255)

     caTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTAc                                                                                      <  2:2654735/93‑1 (MQ=255)

      aTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACt                                                                                     <  2:1630715/93‑1 (MQ=255)

      aTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACt                                                                                     >  1:307672/1‑93 (MQ=255)

       tCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGTCGCAGCTTGTTAAAGAAATCGCTTCTACtg                                                                                    >  2:1927051/1‑93 (MQ=255)

        ccTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACtga                                                                                   <  1:5550678/93‑1 (MQ=255)

         cTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACtgat                                                                                  <  1:877058/93‑1 (MQ=255)

          tGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACtgatg                                                                                 >  2:2863299/1‑93 (MQ=255)

           ggCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACtgatga                                                                                <  2:5047703/93‑1 (MQ=255)

            gCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTgatgat                                                                               >  1:4577240/1‑93 (MQ=255)

             cACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGAtt                                                                              >  1:6347617/1‑93 (MQ=255)

             cACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGAtt                                                                              <  2:1832153/93‑1 (MQ=255)

              aCTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGAttt                                                                             >  2:2025945/1‑93 (MQ=255)

                tGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCt                                                                           <  1:791161/93‑1 (MQ=255)

                 gACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTa                                                                          >  1:6122529/1‑93 (MQ=255)

                   ccaccaACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTAcc                                                                        <  2:5171301/93‑1 (MQ=255)

                    caccaACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCg                                                                       <  2:6568856/93‑1 (MQ=255)

                      ccaACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTc                                                                     <  1:2927692/93‑1 (MQ=255)

                      ccaACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTc                                                                     <  1:6426734/93‑1 (MQ=255)

                      ccaACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTc                                                                     >  2:1857622/1‑93 (MQ=255)

                       caACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCt                                                                    <  2:4169431/93‑1 (MQ=255)

                         aCGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTgg                                                                  >  1:3977906/1‑93 (MQ=255)

                          cGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTggg                                                                 <  2:6007956/93‑1 (MQ=255)

                          cGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTggg                                                                 <  1:5798104/93‑1 (MQ=255)

                           gAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGt                                                                <  2:1716131/93‑1 (MQ=255)

                           gAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGt                                                                >  1:621991/1‑93 (MQ=255)

                           gAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGt                                                                >  1:1786190/1‑93 (MQ=255)

                            aaaaaaCGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTa                                                               <  2:492772/93‑1 (MQ=255)

                             aaaaaCGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTaa                                                              >  1:4598522/1‑93 (MQ=255)

                              aaaaCGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTaaa                                                             >  2:1064595/1‑93 (MQ=255)

                               aaaCGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAg                                                            <  2:2419284/93‑1 (MQ=255)

                                aaCGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGa                                                           <  2:2576984/93‑1 (MQ=255)

                                 aCGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGaa                                                          >  1:4048626/1‑93 (MQ=255)

                                  cGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAAc                                                         <  1:3807947/93‑1 (MQ=255)

                                   ggCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACt                                                        <  1:2017716/93‑1 (MQ=255)

                                   ggCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACt                                                        <  1:2629118/93‑1 (MQ=255)

                                   ggCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACt                                                        <  1:3992543/93‑1 (MQ=255)

                                   ggCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACt                                                        <  2:4855731/93‑1 (MQ=255)

                                     cTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTgg                                                      <  1:23161/93‑1 (MQ=255)

                                      tcatcaGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGc                                                     >  1:5895187/1‑93 (MQ=255)

                                        atcaGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGctc                                                   <  1:2431179/93‑1 (MQ=255)

                                        atcaGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGctc                                                   >  2:2098471/1‑93 (MQ=255)

                                           aGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGATATCGCTTCTTCTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGc                                                <  1:3299714/93‑1 (MQ=255)

                                            gTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCa                                               >  2:1150650/1‑93 (MQ=255)

                                            gTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCa                                               >  1:2961264/1‑93 (MQ=255)

                                            gTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCa                                               <  2:4034275/93‑1 (MQ=255)

                                            gTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCa                                               <  2:2494758/93‑1 (MQ=255)

                                              tGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCaga                                             <  2:3612837/93‑1 (MQ=255)

                                                gTACTGAGGAAAGTCGTTGTCCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGc                                           <  1:5365275/93‑1 (MQ=255)

                                                gTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGc                                           <  1:6031341/93‑1 (MQ=255)

                                                gTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGc                                           <  2:4321831/93‑1 (MQ=255)

                                                   cTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGt                                        >  2:1567433/1‑93 (MQ=255)

                                                   cTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGt                                        <  2:1905017/93‑1 (MQ=255)

                                                     gAGCGAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCt                                      <  2:6302427/93‑1 (MQ=255)

                                                     gAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCt                                      <  2:6246466/93‑1 (MQ=255)

                                                      aGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTg                                     <  1:5647368/93‑1 (MQ=255)

                                                         aaaGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGca                                  <  1:3870430/93‑1 (MQ=255)

                                                         aaaGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGca                                  <  1:4945596/93‑1 (MQ=255)

                                                          aaGGCGTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGcac                                 >  1:4463800/1‑93 (MQ=255)

                                                               gTTGTGCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACaaa                            >  1:2122153/1‑93 (MQ=255)

                                                                 tgtgCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAgg                          <  2:2685003/93‑1 (MQ=255)

                                                                 tgtgCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAgg                          <  2:124330/93‑1 (MQ=255)

                                                                  gtgCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGt                         <  2:5201140/93‑1 (MQ=255)

                                                                   tgCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTg                        <  1:3704832/93‑1 (MQ=255)

                                                                   tgCCGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTg                        >  1:6105859/1‑93 (MQ=255)

                                                                    gCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGa                       <  1:1472480/93‑1 (MQ=255)

                                                                      cGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACg                     >  2:3767877/1‑93 (MQ=255)

                                                                      cGCAGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACg                     <  2:1395764/93‑1 (MQ=255)

                                                                       gCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGt                    <  2:5828255/93‑1 (MQ=255)

                                                                         aGCTTGTTAAAGAGATCGCTTCTGCTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTg                  <  2:1040070/93‑1 (MQ=255)

                                                                         aGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTg                  >  2:4375642/1‑93 (MQ=255)

                                                                         aGCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTg                  <  1:1565621/93‑1 (MQ=255)

                                                                          gCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGt                 >  1:2383861/1‑93 (MQ=255)

                                                                          gCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGt                 >  2:6065728/1‑93 (MQ=255)

                                                                          gCTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGt                 >  1:3414338/1‑93 (MQ=255)

                                                                           cTTGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTa                <  2:508066/93‑1 (MQ=255)

                                                                            ttGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAg               >  2:3352247/1‑93 (MQ=255)

                                                                             tGTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGt              <  2:2156170/93‑1 (MQ=255)

                                                                             tGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGt              >  1:3555719/1‑93 (MQ=255)

                                                                              gTTAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGtt             <  1:1090219/93‑1 (MQ=255)

                                                                               ttAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTa            <  2:4634598/93‑1 (MQ=255)

                                                                                tAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTAt           >  1:1982643/1‑93 (MQ=255)

                                                                                tAAAGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGAGACGTTGTAGGTAt           >  2:5387444/1‑93 (MQ=255)

                                                                                tAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGCAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGCTAt           >  1:4165731/1‑93 (MQ=255)

                                                                                 aaaGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATg          <  2:3047050/93‑1 (MQ=255)

                                                                                  aaGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATgg         >  2:874577/1‑93 (MQ=255)

                                                                                  aaGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATgg         >  2:3217761/1‑93 (MQ=255)

                                                                                  aaGAGATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATgg         >  1:4725739/1‑93 (MQ=255)

                                                                                   agagATCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGt        >  2:5117530/1‑93 (MQ=255)

                                                                                    gagaTCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGtt       >  2:3171817/1‑93 (MQ=255)

                                                                                    gagaTCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGtt       >  1:5074205/1‑93 (MQ=255)

                                                                                      gaTCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTc     >  1:3365090/1‑93 (MQ=255)

                                                                                       aTCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCt    <  1:4415633/93‑1 (MQ=255)

                                                                                       aTCGCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCt    <  1:2475153/93‑1 (MQ=255)

                                                                                         cACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTgg  >  1:1142144/1‑93 (MQ=255)

GCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGG  >  NC_000913/1754901‑1755082

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