BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,077,905	A→G	S584P (TCG→CCG)	tktA ←	transketolase 1 isozyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,077,905	0	A	G	99.4%	1108.2 / ‑5.0	341	S584P (TCG→CCG)	tktA	transketolase 1 isozyme
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (149/190); minor base T (1/1); total (150/191)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 355 total aligned reads displayed.

CCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGACATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGA  >  NC_000913/3077815‑3077994

ccGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCa                                                                                         <  1:8704579/93‑1 (MQ=255)

ccGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCa                                                                                         <  2:2330911/93‑1 (MQ=255)

ccGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCa                                                                                         <  1:3493127/93‑1 (MQ=255)

ccGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCa                                                                                         <  1:3078783/93‑1 (MQ=255)

 cGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCAt                                                                                        <  1:7356481/93‑1 (MQ=255)

 cGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCAt                                                                                        <  1:3961107/93‑1 (MQ=255)

 cGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCAt                                                                                        <  2:7115956/93‑1 (MQ=255)

 cGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCAt                                                                                        <  1:6349597/93‑1 (MQ=255)

  gCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATg                                                                                       <  2:5168854/93‑1 (MQ=255)

   cTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATgg                                                                                      <  2:5573631/93‑1 (MQ=255)

     tCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGac                                                                                    >  2:155947/1‑93 (MQ=255)

      cTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGaca                                                                                   <  1:6505457/93‑1 (MQ=255)

      cTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGaca                                                                                   <  1:2051522/93‑1 (MQ=255)

       tACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGacac                                                                                  >  2:3014373/1‑93 (MQ=255)

        aCAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGAcacc                                                                                 >  2:6697759/1‑93 (MQ=255)

         cAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGAcacca                                                                                <  2:4173766/93‑1 (MQ=255)

          cgCAACGCCTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGAtaccac                                                                               <  1:1625374/92‑1 (MQ=255)

          aGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGAcaccac                                                                               <  1:1541738/93‑1 (MQ=255)

           gCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCAcg                                                                              >  2:6949343/1‑93 (MQ=255)

            cAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCAcgc                                                                             <  2:3810558/93‑1 (MQ=255)

             aaCGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCAcgcg                                                                            <  1:9664249/93‑1 (MQ=255)

             aaCGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCAcgcg                                                                            >  2:9443558/1‑93 (MQ=255)

              aCGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCAcgcgc                                                                           <  1:2433088/93‑1 (MQ=255)

              aCGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCAcgcgc                                                                           >  2:4759852/1‑93 (MQ=255)

               cgcgTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCAcgcgcg                                                                          <  2:6614591/93‑1 (MQ=255)

                  gTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCtt                                                                       >  2:7328942/1‑93 (MQ=255)

                   tGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCttt                                                                      >  2:6573462/1‑93 (MQ=255)

                   tGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCttt                                                                      >  2:82707/1‑93 (MQ=255)

                    gCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACAGCACGCGCGCTTTc                                                                     >  1:9948919/1‑93 (MQ=255)

                    gCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTc                                                                     >  1:2949072/1‑93 (MQ=255)

                     cAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCa                                                                    <  2:8703848/93‑1 (MQ=255)

                      aGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCAc                                                                   <  1:5173574/93‑1 (MQ=255)

                        tAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGc                                                                 >  2:1817285/1‑93 (MQ=255)

                        tAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGc                                                                 >  1:4103730/1‑93 (MQ=255)

                          aCCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCt                                                               <  2:5281692/93‑1 (MQ=255)

                          aCCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCt                                                               >  2:4081626/1‑93 (MQ=255)

                          aCCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCt                                                               >  1:8810884/1‑93 (MQ=255)

                            cGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTc                                                             >  2:6112258/1‑93 (MQ=255)

                               tttCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGc                                                          >  1:5039389/1‑93 (MQ=255)

                                ttCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGca                                                         <  2:8875768/93‑1 (MQ=255)

                                      aGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGcagtcagt                                                   >  2:2728685/1‑93 (MQ=255)

                                      aGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGcagtcagt                                                   >  1:3423201/1‑93 (MQ=255)

                                       gTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGtt                                                  >  2:5913519/1‑93 (MQ=255)

                                        tACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGTCTTCGGCAGTCAGttt                                                 >  2:3537206/1‑93 (MQ=255)

                                        tACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGttt                                                 >  2:6831188/1‑93 (MQ=255)

                                        tACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGttt                                                 >  1:275262/1‑93 (MQ=255)

                                        tACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGttt                                                 >  1:9535491/1‑93 (MQ=255)

                                         aCGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGtttt                                                <  1:1219507/93‑1 (MQ=255)

                                         aCGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGtttt                                                <  1:7993454/93‑1 (MQ=255)

                                          cGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGttttt                                               <  1:7487629/93‑1 (MQ=255)

                                           ggATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTc                                              >  1:9367823/1‑93 (MQ=255)

                                           ggATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTc                                              <  1:1781648/93‑1 (MQ=255)

                                            gATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCg                                             <  2:631028/93‑1 (MQ=255)

                                            gATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCg                                             >  2:2281040/1‑93 (MQ=255)

                                              ttCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTa                                           >  1:7006548/1‑93 (MQ=255)

                                              ttCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTa                                           >  2:262847/1‑93 (MQ=255)

                                                 aCGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGgc                                        <  2:7808711/93‑1 (MQ=255)

                                                  cGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGgca                                       <  2:7980483/93‑1 (MQ=255)

                                                   ggTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGgcag                                      >  1:10001998/1‑93 (MQ=255)

                                                    gTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGgcagc                                     <  1:3765023/93‑1 (MQ=255)

                                                    gTAAGCAGCATCCTGCTTGTCAAATGCGTCGGGAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGgcagc                                     >  2:8419242/1‑93 (MQ=255)

                                                      aaGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCaa                                   <  1:5992591/93‑1 (MQ=255)

                                                       agcagcATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAAc                                  <  1:7676128/93‑1 (MQ=255)

                                                        gcagcaTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACa                                 <  2:1118668/93‑1 (MQ=255)

                                                         cagcaTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAg                                <  1:6733746/93‑1 (MQ=255)

                                                          agcaTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGc                               <  2:1383064/93‑1 (MQ=255)

                                                          agcaTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGc                               >  2:4860154/1‑93 (MQ=255)

                                                          agcaTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGc                               >  2:8521922/1‑93 (MQ=255)

                                                           gcaTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGcc                              <  1:4542393/93‑1 (MQ=255)

                                                            caTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCa                             <  1:3855382/93‑1 (MQ=255)

                                                             aTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAg                            >  2:3282816/1‑93 (MQ=255)

                                                             aTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTATCGTAGGCAGCAACAGCCAg                            <  2:1925187/93‑1 (MQ=255)

                                                              tCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGt                           <  1:5658125/93‑1 (MQ=255)

                                                              tCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGt                           >  1:8961046/1‑93 (MQ=255)

                                                                  gCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCaa                       <  1:9232557/93‑1 (MQ=255)

                                                                  gCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCaa                       <  2:7392138/93‑1 (MQ=255)

                                                                    ttGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACt                     <  1:585539/93‑1 (MQ=255)

                                                                        cAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTg                 >  1:2221603/1‑93 (MQ=255)

                                                                        cAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTg                 >  2:3184974/1‑93 (MQ=255)

                                                                         aaaTGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGa                <  1:5632508/93‑1 (MQ=255)

                                                                         aaaTGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGa                >  1:4747555/1‑93 (MQ=255)

                                                                          aaTGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGaa               <  2:5106751/93‑1 (MQ=255)

                                                                          aaTGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGaa               >  2:2474656/1‑93 (MQ=255)

                                                                          aaTGCGTCGATAGACGGCATGGACACCACGCGCGCTTTCACTCCTTCGGCAGTCAGTTTTTCTTAGGCAGCAACAGCCAGTTCAACTTCTGaa               <  1:5773139/93‑1 (MQ=255)

                                                                            ttcGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAAcc             <  2:9653275/91‑1 (MQ=255)

                                                                            tGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAAcc             >  2:7657597/1‑93 (MQ=255)

                                                                               ntNGGTANACGGCATGGACACNCCGCGGTCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGt          <  1:8355824/92‑1 (MQ=255)

                                                                               gTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGt          <  2:4277412/93‑1 (MQ=255)

                                                                               gTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGt          >  2:1262330/1‑93 (MQ=255)

                                                                                tCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTa         <  2:4686137/93‑1 (MQ=255)

                                                                                 cGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAg        <  1:6098227/93‑1 (MQ=255)

                                                                                 cGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAg        <  1:7099236/93‑1 (MQ=255)

                                                                                  ggTAGGCGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGc       <  1:5411665/93‑1 (MQ=255)

                                                                                  ggTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGc       <  1:4313455/93‑1 (MQ=255)

                                                                                  ggTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGc       <  1:1369031/93‑1 (MQ=255)

                                                                                  ggTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGc       <  1:9788708/93‑1 (MQ=255)

                                                                                   gTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCg      >  1:4911113/1‑93 (MQ=255)

                                                                                    tAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGa     <  2:8052943/93‑1 (MQ=255)

                                                                                      gacggcATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACATCCAGTTCAACTTCTGAACCGGTAGCGATg   >  2:9207749/1‑93 (MQ=255)

                                                                                      gacggcATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTACGGATg   >  1:3279505/1‑93 (MQ=255)

CCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGACATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGA  >  NC_000913/3077815‑3077994

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