BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,554,118	+C	coding (367/396 nt)	yhgKJ ←	RNA phosphage cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,554,117	1	.	C	98.6%	438.6 / ‑1.1	143	coding (368/396 nt)	yhgKJ	RNA phosphage cyclase
Reads supporting (aligned to +/- strand): ref base . (2/0); new base C (57/84); total (59/84)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 163 total aligned reads displayed.

GGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGAC‑GGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGGC  >  NC_000913/3554027‑3554199

ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGt‑gccg‑cc                                                                                  >  1:693952/1‑92 (MQ=255)

ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGt‑gccg‑cc                                                                                  >  2:3004909/1‑92 (MQ=255)

 gTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCCGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGGTATCGGt‑gccg‑ccg                                                                                 >  1:4096567/1‑91 (MQ=255)

 gTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCCGCGGATAAAATCCGCAAGCGGCGCCGAC‑GGGCTATCGGt‑gccg‑ccg                                                                                 >  2:4302840/1‑91 (MQ=255)

  tttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGAAGGCGGCGCCGAC‑GGGTTATCGGt‑gccg‑ccgc                                                                                <  1:4251317/93‑4 (MQ=255)

   ttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGt‑gccg‑ccgct                                                                               <  1:3119883/93‑5 (MQ=255)

      cTGATGAATTCCTTTTTCCGCCAGCAGCGGCTCCAGCGCCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCAc                                                                             <  1:4654931/93‑1 (MQ=255)

      cTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCAc                                                                             <  2:3195842/93‑1 (MQ=255)

      cTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCAc                                                                             <  1:221190/93‑1 (MQ=255)

       tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACt                                                                            >  2:3313969/1‑93 (MQ=255)

       tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACt                                                                            >  1:1687361/1‑93 (MQ=255)

       tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACt                                                                            >  2:4180367/1‑93 (MQ=255)

        gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACtt                                                                           >  1:1256687/1‑93 (MQ=255)

        gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACtt                                                                           >  1:598811/1‑93 (MQ=255)

        gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGGTATCGGT‑GCCGCCGCTCACtt                                                                           >  1:4193217/1‑93 (MQ=255)

        gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGAGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACtt                                                                           >  2:4344219/1‑93 (MQ=255)

        gatgaATCCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACtt                                                                           <  1:890229/93‑1 (MQ=255)

         atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTc                                                                          >  2:3220426/1‑93 (MQ=255)

           gaATTCCTATTTTCGCCAGCAGCGGCTCTAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCaa                                                                        >  1:1507279/1‑93 (MQ=255)

           gaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCaa                                                                        <  2:1684647/93‑1 (MQ=255)

           gaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCCCCGCC‑GGGTTATCGGT‑GCCGCCGCTCACTTCaa                                                                        >  1:663363/1‑93 (MQ=255)

            aaTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAac                                                                       >  2:3833969/1‑93 (MQ=255)

             aTTCCTATTTTCGCGAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAaca                                                                      <  2:361236/93‑1 (MQ=255)

             aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGTTCCGC‑CGCTCACTTCAaca                                                                      >  1:3831162/1‑93 (MQ=255)

             aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAaca                                                                      >  2:714295/1‑93 (MQ=255)

             aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAaca                                                                      <  2:1890220/93‑1 (MQ=255)

              ttCCTATTTTCGCCAGCAGCGGCTCCAGCACCCCGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGGTATCGG‑CGCCGCCGCTCACTTCAacac                                                                     >  1:3076946/1‑93 (MQ=255)

                 cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGGTATCGGT‑GCCGCCGCTGACTTCAACAcgcg                                                                  >  1:3791446/1‑93 (MQ=255)

                  tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGa                                                                 <  1:2958011/93‑1 (MQ=255)

                  tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGa                                                                 >  2:3935684/1‑93 (MQ=255)

                   aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGaa                                                                >  1:1889856/1‑93 (MQ=255)

                    ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAg                                                               >  2:2285909/1‑93 (MQ=255)

                    ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAg                                                               <  2:4380258/93‑1 (MQ=255)

                    tgttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAg                                                               <  2:3368299/91‑1 (MQ=255)

                     tttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAgg                                                              <  1:2760043/93‑1 (MQ=255)

                     tttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAgg                                                              <  1:3353429/93‑1 (MQ=255)

                        caccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTcc                                                           <  2:4598739/91‑1 (MQ=255)

                         cccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCa                                                          <  2:2316461/92‑1 (MQ=255)

                          gcagcagCGGCGCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                                         <  1:786299/92‑1 (MQ=255)

                          ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                                         <  1:2457941/93‑1 (MQ=255)

                          ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                                         <  2:2732523/93‑1 (MQ=255)

                           cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGTGTCGC‑CGCTGACTTCAACACGCGAAGAGCCATc                                                        >  2:4547664/1‑93 (MQ=255)

                           cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATc                                                        <  2:2953682/93‑1 (MQ=255)

                           cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATc                                                        <  1:3583968/93‑1 (MQ=255)

                            agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                                                       <  1:3723252/93‑1 (MQ=255)

                             gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                                                      >  2:2525845/1‑93 (MQ=255)

                             gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                                                      >  1:3654856/1‑93 (MQ=255)

                               agcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCa                                                    <  2:1774752/93‑1 (MQ=255)

                                gcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaa                                                   >  2:482215/1‑93 (MQ=255)

                                   gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                                                >  2:1730319/1‑93 (MQ=255)

                                   gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                                                <  2:119497/93‑1 (MQ=255)

                                   gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                                                <  1:3937920/93‑1 (MQ=255)

                                    cTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCACACca                                               >  1:1112556/1‑93 (MQ=255)

                                     tCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcac                                              <  2:2133212/93‑1 (MQ=255)

                                      ccAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAACGTCCATCGGCAAACcaca                                             >  1:2376299/1‑93 (MQ=255)

                                       cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                                            <  1:1473037/93‑1 (MQ=255)

                                        aGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                                           >  2:3403407/1‑93 (MQ=255)

                                         gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcg                                          <  2:603072/93‑1 (MQ=255)

                                         gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcg                                          <  1:3302346/93‑1 (MQ=255)

                                            cccGGCGGGAAAAATCCGCAGGCGGCGCCGGCGGGGTTAT‑GGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCgg                                       <  1:4493988/93‑1 (MQ=255)

                                             ccGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCggg                                      <  1:707844/93‑1 (MQ=255)

                                             ccGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCggg                                      >  1:3454541/1‑93 (MQ=255)

                                             ccGGCGGATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCggg                                      >  1:4540056/1‑93 (MQ=255)

                                              cggcggATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc                                     >  2:3299450/1‑93 (MQ=255)

                                               ggcggATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgca                                    >  2:435680/1‑93 (MQ=255)

                                                gcggATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcag                                   >  2:821224/1‑93 (MQ=255)

                                                 cggATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagc                                  >  2:1498242/1‑93 (MQ=255)

                                                  ggATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGTCCCGC‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca                                 >  1:2954037/1‑93 (MQ=255)

                                                  ggATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca                                 >  2:1577049/1‑93 (MQ=255)

                                                   gATAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc                                >  2:308630/1‑93 (MQ=255)

                                                     tAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCg                              <  1:4048654/93‑1 (MQ=255)

                                                     tAAAATCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCg                              <  2:2241888/93‑1 (MQ=255)

                                                       aaaTCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGGAAACCACAGCGCGGGCCGCACCGtt                            >  2:249021/1‑93 (MQ=255)

                                                       aaaTCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAANCCACAGCGCGGGCAGCACCNtt                            >  1:3311862/1‑93 (MQ=255)

                                                       aaaTCCGCAGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt                            >  1:99594/1‑93 (MQ=255)

                                                             gcaGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgcagc                      <  1:687445/91‑1 (MQ=255)

                                                              caGGCGGCGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgcagca                     <  1:2901959/92‑1 (MQ=255)

                                                               cggcggcGCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAc                    <  1:1358616/92‑1 (MQ=255)

                                                                  cggcgCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAg                 <  2:2428065/93‑1 (MQ=255)

                                                                  cggcgCCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAg                 <  1:2704207/93‑1 (MQ=255)

                                                                    gcgcCGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCg               <  1:4273601/93‑1 (MQ=255)

                                                                       ccGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTAc            >  2:2046265/1‑93 (MQ=255)

                                                                       ccGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTAc            <  1:680741/93‑1 (MQ=255)

                                                                       ccGAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTAc            <  1:2877747/93‑1 (MQ=255)

                                                                         gAC‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACaa          >  1:226068/1‑93 (MQ=255)

                                                                           c‑GGTTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACt        <  2:3556924/93‑1 (MQ=255)

                                                                           c‑GGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACt        <  1:1742458/93‑1 (MQ=255)

                                                                             gggTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACtt       <  1:874692/93‑1 (MQ=255)

                                                                              ggTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTc      >  2:3574666/1‑93 (MQ=255)

                                                                              ggTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTc      <  2:3776255/93‑1 (MQ=255)

                                                                              ggTGATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTc      <  1:3714123/93‑1 (MQ=255)

                                                                               gTTATCGGTTCCTC‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTcc     <  2:3484842/93‑1 (MQ=255)

                                                                               gTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTcc     <  2:3493708/93‑1 (MQ=255)

                                                                               gTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTcc     <  2:2099889/93‑1 (MQ=255)

                                                                               gTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTcc     >  2:2789970/1‑93 (MQ=255)

                                                                                ttATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCg    <  1:2494374/93‑1 (MQ=255)

                                                                                ttATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCg    <  2:2383889/93‑1 (MQ=255)

                                                                                 tATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCgg   <  2:1879595/93‑1 (MQ=255)

                                                                                 tATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCgg   <  2:2575326/93‑1 (MQ=255)

                                                                                  aTCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGgc  <  2:2532608/93‑1 (MQ=255)

GGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGAC‑GGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGGC  >  NC_000913/3554027‑3554199

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