BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,579,447	T→C	S598P (TCG→CCG)	tktB →	transketolase 2 isozyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,579,447	0	T	C	97.8%	427.4 / ‑0.5	139	S598P (TCG→CCG)	tktB	transketolase 2 isozyme
Reads supporting (aligned to +/- strand): ref base T (3/0); new base C (70/66); total (73/66)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01

Only 100 of 143 total aligned reads displayed.

GAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAA  >  NC_000913/2579357‑2579537

gagaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCg                                                                                          >  1:3651982/1‑93 (MQ=255)

 agaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCgg                                                                                         <  2:445758/93‑1 (MQ=255)

 agaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCgg                                                                                         >  2:612712/1‑93 (MQ=255)

 agaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCgg                                                                                         >  2:1096336/1‑93 (MQ=255)

 agaAAATAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCgg                                                                                         <  1:1729595/93‑1 (MQ=255)

  gaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgt                                                                                        <  2:2512481/93‑1 (MQ=255)

   aaaTTAG‑AGGAGAAGGTCGCAATGCACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgtgt                                                                                      <  2:2315443/93‑1 (MQ=255)

     aTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGtt                                                                                     >  1:1017927/1‑93 (MQ=255)

         gCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCt                                                                                 <  2:2718609/93‑1 (MQ=255)

          cAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCtt                                                                                <  2:106618/93‑1 (MQ=255)

           aGGAGAAGGTGGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCAttt                                                                               >  2:1083401/1‑92 (MQ=255)

           aGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTc                                                                               >  1:3039484/1‑93 (MQ=255)

           aGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTc                                                                               <  2:4329021/93‑1 (MQ=255)

           aGGAGAAGGTCGCAATGCACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGTGAGCCGGTGTTGCCTTc                                                                               <  1:3602059/93‑1 (MQ=255)

            ggAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCt                                                                              >  1:4531345/1‑93 (MQ=255)

            ggAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCt                                                                              <  1:4169905/93‑1 (MQ=255)

              agaAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTaa                                                                            >  1:2323814/1‑93 (MQ=255)

                aaGGTCGCAATGGACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACg                                                                          <  1:1063759/93‑1 (MQ=255)

                 aGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGt                                                                         <  2:2740147/93‑1 (MQ=255)

                  ggTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGtt                                                                        <  1:2032704/93‑1 (MQ=255)

                  ggTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGtt                                                                        <  2:541593/93‑1 (MQ=255)

                   gTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTg                                                                       >  2:4600414/1‑93 (MQ=255)

                   gTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTg                                                                       <  2:2550775/93‑1 (MQ=255)

                     cGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGCTGCg                                                                     >  1:2261767/1‑93 (MQ=255)

                      gCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCgg                                                                    >  2:4428381/1‑93 (MQ=255)

                       cAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGc                                                                   <  2:2105639/93‑1 (MQ=255)

                       cAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGc                                                                   >  1:327774/1‑93 (MQ=255)

                       cAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGc                                                                   <  1:349869/93‑1 (MQ=255)

                        aaTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCt                                                                  >  2:1402824/1‑93 (MQ=255)

                        aaTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCt                                                                  >  2:3866890/1‑93 (MQ=255)

                        aaTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCt                                                                  >  1:1975767/1‑93 (MQ=255)

                         aTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTc                                                                 >  1:2889348/1‑93 (MQ=255)

                          tGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTcg                                                                >  1:1687752/1‑93 (MQ=255)

                          tGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTcg                                                                <  1:2367957/93‑1 (MQ=255)

                              cgcgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTg                                                            <  2:4167080/93‑1 (MQ=255)

                               gcgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCCCgggg                                                           >  2:4153658/1‑93 (MQ=255)

                                 gTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCg                                                         >  2:4559071/1‑93 (MQ=255)

                                   aGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGt                                                       <  2:879564/93‑1 (MQ=255)

                                    gTTTCCCTGCCCTCGACCGATATTTTCGACGCCCATGATGAGGAATATCGGGAGCCGGGGTTGCCTTCTAACGTTGCGGCTCGCATTGCGGTg                                                      >  2:4236408/1‑93 (MQ=255)

                                     tttCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTgg                                                     <  2:2209218/93‑1 (MQ=255)

                                     tttCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTgg                                                     <  1:2839486/93‑1 (MQ=255)

                                     tttCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGCGgggg                                                     >  1:755184/1‑93 (MQ=255)

                                     tttCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGCGGTgg                                                     >  2:1376789/1‑93 (MQ=255)

                                       tCCCTGCCCTCGACCGATATTTTCGACGCCCAGGAGGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGaa                                                   <  2:1910422/93‑1 (MQ=255)

                                        cccTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAg                                                  <  1:710962/93‑1 (MQ=255)

                                         ccTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCCCGTGGCGGGGGGAgg                                                 >  2:3874174/1‑92 (MQ=255)

                                          cTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCa                                                <  1:1085654/93‑1 (MQ=255)

                                          cTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCa                                                <  1:2126729/93‑1 (MQ=255)

                                              ccTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGCGGTGGGAGGCGGTa                                            >  2:1387571/1‑93 (MQ=255)

                                              ccTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCCCGTGGCGGGGGGaggaggga                                            >  2:245021/1‑91 (MQ=255)

                                               cTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTAt                                           <  1:2729778/93‑1 (MQ=255)

                                               cTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGAAGCCGGGAt                                           >  2:3461890/1‑93 (MQ=255)

                                                tCGACCGATATTTTTGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGGAtt                                          >  1:1964540/1‑93 (MQ=255)

                                                tCGACCGATATTTTCGACGCCGAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTAtt                                          <  2:732433/93‑1 (MQ=255)

                                                tCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTAtt                                          <  1:1191750/93‑1 (MQ=255)

                                                 cGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCCCGGGGGGGTGGGAGGAGGTATTg                                         >  2:4577447/1‑93 (MQ=255)

                                                  gACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGGAGGGAtttc                                        >  2:2010399/1‑91 (MQ=255)

                                                  gACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGc                                        <  2:3336799/93‑1 (MQ=255)

                                                  gACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGAGGCAGGGATTGc                                        >  1:4526221/1‑93 (MQ=255)

                                                   aCCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGcc                                       >  2:187280/1‑93 (MQ=255)

                                                    ccGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCg                                      <  1:3557619/93‑1 (MQ=255)

                                                    ccGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGCGGTGGAAGCAGGTATTGCCg                                      >  2:509444/1‑93 (MQ=255)

                                                     cGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGCGGGGGGAGCAGGTATTGCCGa                                     >  1:4390991/1‑93 (MQ=255)

                                                     cGATATTTGCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGa                                     <  1:1504253/93‑1 (MQ=255)

                                                      gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGGATTGCCGAt                                    >  2:1370278/1‑93 (MQ=255)

                                                       atatTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGTATTGCCGAtt                                   >  1:2327051/1‑93 (MQ=255)

                                                       atatTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGAtt                                   <  1:4260824/93‑1 (MQ=255)

                                                         atTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTAc                                 <  1:575662/93‑1 (MQ=255)

                                                          ttttCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACt                                <  2:222537/93‑1 (MQ=255)

                                                           tttCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGGATTGCCGATTACTg                               >  1:1304419/1‑93 (MQ=255)

                                                             tCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCCATTACTGGt                             >  1:924659/1‑93 (MQ=255)

                                                              cGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTa                            <  2:1510309/93‑1 (MQ=255)

                                                              cGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGCGGGGGAAGCAGGGATTGCCGATTACTggga                            >  1:856712/1‑91 (MQ=255)

                                                                   cccAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGt                       <  1:2269831/93‑1 (MQ=255)

                                                                    ccAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGGCGACTACTGGTACAAGTa                      >  1:377161/1‑93 (MQ=255)

                                                                     cAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTAt                     <  1:4358124/93‑1 (MQ=255)

                                                                       ggATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGt                   <  2:1416809/93‑1 (MQ=255)

                                                                        gATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGGATTGCCGATTACTTGTACAAGTATGtt                  >  2:934821/1‑93 (MQ=255)

                                                                        gATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGtt                  <  1:2695418/93‑1 (MQ=255)

                                                                         aTGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGTATTGCCGATTACTGGTACAAGTATGTTg                 >  1:3199969/1‑93 (MQ=255)

                                                                         aTGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGGATTGCCGATTACTGGTACAAGTATGTTg                 >  1:4100889/1‑93 (MQ=39)

                                                                         aTGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTg                 <  2:3638406/93‑1 (MQ=255)

                                                                         aTGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTg                 <  2:4431854/93‑1 (MQ=255)

                                                                         aTGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTg                 <  2:168824/93‑1 (MQ=255)

                                                                         aTGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGAAGCAGGAATTGCCGATTACTGGTACAAGTATGTTg                 >  1:419131/1‑93 (MQ=255)

                                                                          tGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGGATTGCCGATTACTGGTACAAGTATGTTgg                >  1:3592909/1‑93 (MQ=39)

                                                                          tGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGAAGCAGGTATTGCCGATTTCTGGTACAAGTATGTTgg                >  2:2223336/1‑93 (MQ=255)

                                                                            aGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTc              <  2:3834016/93‑1 (MQ=255)

                                                                            aGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTc              <  2:4422956/93‑1 (MQ=255)

                                                                               aaTATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGAAGCAGGTATTTCCGATTACTGGTACCAGTATGTTGGTCTGa           >  2:4049951/1‑93 (MQ=255)

                                                                                 tatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGGGGGGGGAGCAGGGATTGCCGATTACTGGTACAAGTATGTTGGTCTGtaa         >  2:4233813/1‑93 (MQ=38)

                                                                                 tatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaaa         >  1:1778732/1‑93 (MQ=255)

                                                                                 tatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTCGCGGTGGGAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaaa         >  1:576825/1‑93 (MQ=255)

                                                                                   tCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGGGGGGGAAGCAGGGATTGCCGATTACTGGTACAAGTATGTTTGTCTGGAAgg       >  1:2359473/1‑93 (MQ=38)

                                                                                   tCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGGTTACTGGTACAAGTATGTTGGTCTGaaaag       >  1:930552/1‑91 (MQ=255)

                                                                                   tCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAgg       <  2:1013807/93‑1 (MQ=255)

                                                                                   tCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAgg       <  1:1474311/93‑1 (MQ=255)

                                                                                    cGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAggg      <  1:557933/93‑1 (MQ=255)

                                                                                      ggAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGc    <  2:236156/93‑1 (MQ=255)

                                                                                      ggAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGc    <  1:1717809/93‑1 (MQ=255)

GAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAA  >  NC_000913/2579357‑2579537

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