BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,069,904	T→C	T247A (ACC→GCC)	pgk ←	phosphoglycerate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,069,904	0	T	C	95.2%	962.0 / 23.3	333	T247A (ACC→GCC)	pgk	phosphoglycerate kinase
Reads supporting (aligned to +/- strand): ref base T (9/7); new base C (140/175); total (151/182)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 350 total aligned reads displayed.

CATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATAAAGGTGTTAGCGA  >  NC_000913/3069814‑3069993

caTCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCca                                                                                         <  1:260307/93‑1 (MQ=255)

caTCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCca                                                                                         >  2:1001567/1‑93 (MQ=255)

 aTCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCcag                                                                                        >  1:2854229/1‑93 (MQ=255)

  tCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCcagc                                                                                       >  2:4831052/1‑93 (MQ=255)

  tCGTGAACAGATTTGAGGGTAGCCGGTGCAGTTTCATAGACCTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCcagc                                                                                       <  1:9268498/93‑1 (MQ=255)

     ttAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGa                                                                                    <  2:6727893/93‑1 (MQ=255)

     ttAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGa                                                                                    <  1:7558631/93‑1 (MQ=255)

     ttAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGa                                                                                    <  2:3790295/93‑1 (MQ=255)

      tAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGAc                                                                                   >  2:7711143/1‑93 (MQ=255)

        aCAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGt                                                                                 <  1:3513251/93‑1 (MQ=255)

           gATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTa                                                                              >  2:6036119/1‑93 (MQ=255)

           gATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTa                                                                              <  2:6643981/93‑1 (MQ=255)

           gATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTa                                                                              <  1:9072287/93‑1 (MQ=255)

           gATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTa                                                                              >  1:8685892/1‑93 (MQ=255)

           gATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTa                                                                              <  1:3768145/93‑1 (MQ=255)

           gATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTa                                                                              <  1:4728947/93‑1 (MQ=255)

            aTTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAg                                                                             <  1:1241641/93‑1 (MQ=255)

             tttCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGc                                                                            >  1:5265386/1‑93 (MQ=255)

              ttCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCt                                                                           >  1:174246/1‑93 (MQ=255)

              ttCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCt                                                                           <  2:2535989/93‑1 (MQ=255)

               tCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCtt                                                                          <  1:2428656/93‑1 (MQ=255)

               tCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCtt                                                                          <  1:6964730/93‑1 (MQ=255)

                cAGGGTAGCCGGTGCAGTTTCAGAGACCTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTc                                                                         <  2:316280/93‑1 (MQ=255)

                cAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTc                                                                         <  2:9821581/93‑1 (MQ=255)

                  gggTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGt                                                                       >  1:7249180/1‑93 (MQ=255)

                   ggTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTc                                                                      >  1:444194/1‑93 (MQ=255)

                   ggTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTc                                                                      <  2:5283860/93‑1 (MQ=255)

                     tAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCaa                                                                    >  2:6341758/1‑93 (MQ=255)

                       gCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAAcc                                                                  <  1:4135688/93‑1 (MQ=255)

                        ccGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGACCCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCa                                                                 <  2:3907156/93‑1 (MQ=255)

                        ccGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGTCAACCa                                                                 <  2:4013691/93‑1 (MQ=255)

                        ccGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCa                                                                 <  2:2072398/93‑1 (MQ=255)

                         cGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAg                                                                <  1:3199714/93‑1 (MQ=255)

                         cGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAg                                                                >  1:4000374/1‑93 (MQ=255)

                          ggTGAAGTCTCAGAGAACTCGGTTGCTACGCTAACATCAGACTGAACCGGGATGTTGCATGTGGCCAGCAGACGTTTAGCTTCGTCAACCAgg                                                               <  2:2945599/93‑1 (MQ=255)

                           gTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGt                                                              <  2:4442418/93‑1 (MQ=255)

                           gTGCAGTTTCAGAGAACTCGGTCCCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGt                                                              <  2:5590419/93‑1 (MQ=255)

                            tGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTc                                                             <  1:5113189/93‑1 (MQ=255)

                            tGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTc                                                             <  2:9236497/93‑1 (MQ=255)

                             gCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCa                                                            >  1:9797964/1‑93 (MQ=255)

                                 tttCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCtt                                                        <  1:4345447/93‑1 (MQ=255)

                                 tttCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCtt                                                        >  1:717799/1‑93 (MQ=255)

                                 tttCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTNANCTTCGTCAACcagntcagntn                                                        >  1:7999453/1‑92 (MQ=255)

                                  ttCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTc                                                       <  2:3409308/93‑1 (MQ=255)

                                  ttCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTc                                                       <  2:7213983/93‑1 (MQ=255)

                                    cAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGt                                                     >  2:3639247/1‑93 (MQ=255)

                                    cAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGt                                                     <  1:8141077/93‑1 (MQ=255)

                                      gagaACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGGCAACCAGGTCAGCTTCATAc                                                   >  1:2147312/1‑93 (MQ=255)

                                       agaACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACa                                                  <  1:6419122/93‑1 (MQ=255)

                                        gaACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAg                                                 >  2:9386732/1‑93 (MQ=255)

                                         aaCTCGGTTGCTACGCGAACATCAAACGGAACCGGGATGTTGCAGGTGGACAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACacg                                                >  2:7035469/1‑91 (MQ=255)

                                          aCTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAggg                                               <  1:5471327/93‑1 (MQ=255)

                                           cTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGa                                              >  1:2560663/1‑93 (MQ=255)

                                               gTTGCTACGCGAACATCAGACGGCACCGGGATGGTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCACCTTCGGACAGGGAtttc                                          >  1:9487685/1‑92 (MQ=255)

                                               gTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTa                                          >  1:7899788/1‑93 (MQ=255)

                                                ttGCTACGCGAACATCAGACGGAACCGGTATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTAc                                         <  1:1074447/93‑1 (MQ=255)

                                                 tGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTAcc                                        <  2:3728177/93‑1 (MQ=255)

                                                  gCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTAccc                                       >  1:8482257/1‑93 (MQ=255)

                                                   cTACGCGAACATCAGACGGAACCGGGATGTTGCAGGGGGCCAGCAGACGGTTAGCTTCGGCAACCAGGTCAGCTTCGTACAGGGATTTACCca                                      >  1:710569/1‑93 (MQ=255)

                                                    tACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCcac                                     <  1:3118986/93‑1 (MQ=255)

                                                     aCGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCcaca                                    <  1:6830149/93‑1 (MQ=255)

                                                       gcgAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATc                                  <  2:7861778/93‑1 (MQ=255)

                                                       gcgAACATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATc                                  >  2:1149808/1‑93 (MQ=255)

                                                          aaCATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTg                               <  1:1911650/93‑1 (MQ=255)

                                                          aaCATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTg                               >  2:6118564/1‑93 (MQ=255)

                                                          aaCATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTg                               >  1:8368847/1‑93 (MQ=255)

                                                            cATCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGc                             <  2:7411934/93‑1 (MQ=255)

                                                             aTCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGcc                            <  2:8042504/93‑1 (MQ=255)

                                                              tCAGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCt                           >  2:2757765/1‑93 (MQ=255)

                                                                aGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTtg                         <  1:1636840/93‑1 (MQ=255)

                                                                aGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTtg                         >  2:9623114/1‑93 (MQ=255)

                                                                aGACGGAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTtg                         <  1:8888725/93‑1 (MQ=255)

                                                                    ggAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGcc                     >  1:6570001/1‑93 (MQ=255)

                                                                    ggAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGcc                     <  2:8740015/93‑1 (MQ=255)

                                                                     gAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCcg                    >  1:4672878/1‑93 (MQ=255)

                                                                     gAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCcg                    >  1:5890016/1‑93 (MQ=255)

                                                                     gAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCcg                    <  2:5075704/93‑1 (MQ=255)

                                                                     gAACCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCcg                    >  2:1862270/1‑93 (MQ=255)

                                                                      aaCCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCcgc                   <  1:4897629/93‑1 (MQ=255)

                                                                       aCCGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCcgcg                  >  2:1470327/1‑93 (MQ=255)

                                                                        ccGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGa                 >  2:521518/1‑93 (MQ=255)

                                                                        ccGGGATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGa                 >  2:6498041/1‑93 (MQ=255)

                                                                           ggATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATaa              >  2:8295113/1‑93 (MQ=255)

                                                                           ggATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATaa              >  2:3267512/1‑93 (MQ=255)

                                                                            gATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATaaa             >  1:5952894/1‑93 (MQ=255)

                                                                            gATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATaaa             <  1:3316123/93‑1 (MQ=255)

                                                                            gATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATaaa             <  1:1545974/93‑1 (MQ=255)

                                                                            gATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATaaa             <  2:6901417/93‑1 (MQ=255)

                                                                            gATGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATaaa             >  2:2745679/1‑93 (MQ=255)

                                                                             aTGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATAAAg            <  2:5341359/93‑1 (MQ=255)

                                                                              tGTTGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATAAAgg           <  1:5551607/93‑1 (MQ=255)

                                                                                ttGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATAAAGgtg         >  2:4619363/1‑93 (MQ=255)

                                                                                ttGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATAAAGgtg         >  2:1644222/1‑93 (MQ=255)

                                                                                 tGCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATAAAGgtgt        <  2:868443/93‑1 (MQ=255)

                                                                                  gCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATAAAGGTGtt       <  1:6795863/93‑1 (MQ=255)

                                                                                  gCAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATAAAGGTGtt       >  2:4176658/1‑93 (MQ=255)

                                                                                   cAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATAAAGGTGTTa      <  1:6494362/93‑1 (MQ=255)

                                                                                   cAGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATAAAGGTGTTa      <  2:242855/93‑1 (MQ=255)

                                                                                    aGGTGGCCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATAAAGGTGTTAg     <  2:4208027/93‑1 (MQ=255)

                                                                                       tggtCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATAAAGGTGTTAGCGa  >  2:5840193/1‑93 (MQ=255)

CATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATAAAGGTGTTAGCGA  >  NC_000913/3069814‑3069993

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