BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,965,001	T→C	M486V (ATG→GTG)	ptsP ←	PtsP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,965,001	0	T	C	99.5%	662.8 / ‑4.8	213	M486V (ATG→GTG)	ptsP	PtsP
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (95/116); total (97/116)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

Only 100 of 216 total aligned reads displayed.

CTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCAT  >  NC_000913/2964911‑2965090

  tCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACcccc                                                                                       <  1:911040/93‑1 (MQ=255)

    gCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTg                                                                                     >  1:1460281/1‑93 (MQ=255)

     cTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGa                                                                                    <  1:1494323/93‑1 (MQ=255)

     cTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGa                                                                                    <  1:510888/93‑1 (MQ=255)

       gATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATa                                                                                  >  1:2340464/1‑93 (MQ=255)

       gATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATa                                                                                  >  1:1652864/1‑93 (MQ=255)

       gATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATa                                                                                  >  1:1390991/1‑93 (MQ=255)

       gATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATa                                                                                  <  2:4216968/93‑1 (MQ=255)

        gtCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCCCCCCCTGATAc                                                                                 <  2:3178841/92‑1 (MQ=255)

        aTCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATAc                                                                                 <  1:3719708/93‑1 (MQ=255)

          cGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTg                                                                               >  2:2814238/1‑93 (MQ=255)

               tGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCa                                                                          >  1:1146637/1‑93 (MQ=255)

               tGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCa                                                                          >  1:241419/1‑93 (MQ=255)

                gTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCAc                                                                         >  1:3218635/1‑93 (MQ=255)

                 tAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCAcc                                                                        <  1:826451/93‑1 (MQ=255)

                  aaGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCt                                                                       >  2:2571025/1‑93 (MQ=255)

                   aGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCCAGGGGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTg                                                                      >  2:4191715/1‑93 (MQ=255)

                   aGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTg                                                                      >  1:3771257/1‑93 (MQ=255)

                    ggCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGt                                                                     >  1:1122758/1‑93 (MQ=255)

                     gCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGtt                                                                    >  2:3654773/1‑93 (MQ=255)

                     gCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGtt                                                                    <  1:204294/93‑1 (MQ=255)

                         ctgctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTGTCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttc                                                                <  1:2588328/93‑1 (MQ=255)

                          tgctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttct                                                               >  2:3026245/1‑93 (MQ=255)

                          tgctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttct                                                               <  2:3681156/93‑1 (MQ=255)

                          tgctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttct                                                               >  1:4667718/1‑93 (MQ=255)

                          tgctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttct                                                               <  2:385418/93‑1 (MQ=255)

                          tgctTATCTGCTCCGACATCCAGCGTACGCAAGGGGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttct                                                               >  1:3908241/1‑93 (MQ=255)

                           gctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttctt                                                              >  1:4484299/1‑93 (MQ=255)

                           gctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttctt                                                              >  1:3810932/1‑93 (MQ=255)

                             ttATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTcc                                                            <  1:1702512/93‑1 (MQ=255)

                             ttATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTcc                                                            >  1:981096/1‑93 (MQ=255)

                             ttATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTcc                                                            <  1:56757/93‑1 (MQ=255)

                              tATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCg                                                           >  2:454817/1‑93 (MQ=255)

                               aTCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGa                                                          >  1:3245917/1‑93 (MQ=255)

                                tCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCATCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGAc                                                         <  2:2235876/93‑1 (MQ=255)

                                tCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGAc                                                         <  1:1001826/93‑1 (MQ=255)

                                tCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGAc                                                         >  2:1875057/1‑93 (MQ=255)

                                  tGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACgg                                                       <  2:1054078/93‑1 (MQ=255)

                                      ccGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGaaa                                                   <  2:2177864/93‑1 (MQ=255)

                                       cGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAAc                                                  <  1:4437300/93‑1 (MQ=255)

                                        gACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAAcc                                                 <  1:1998456/93‑1 (MQ=255)

                                        gACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAAcc                                                 >  1:792688/1‑93 (MQ=255)

                                           aTCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACt                                              <  2:1426533/93‑1 (MQ=255)

                                           aTCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACt                                              >  2:1027541/1‑93 (MQ=255)

                                           aTCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACt                                              <  2:2739988/93‑1 (MQ=255)

                                           aTCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACt                                              >  2:38250/1‑93 (MQ=255)

                                           aTCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACt                                              <  2:4245399/93‑1 (MQ=255)

                                           aTCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACt                                              >  1:102746/1‑93 (MQ=255)

                                            tCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACtt                                             >  1:4240761/1‑93 (MQ=255)

                                             ccAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACttt                                            <  1:1069138/93‑1 (MQ=255)

                                             ccAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACttt                                            <  2:2723159/93‑1 (MQ=255)

                                                gCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTgca                                         >  1:3732125/1‑93 (MQ=255)

                                                gCGTACGCAAGGTGACAGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTgca                                         <  2:2490702/93‑1 (MQ=255)

                                                  gTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTgcagc                                       >  2:2978200/1‑93 (MQ=255)

                                                  gTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTgcagc                                       <  2:2002542/93‑1 (MQ=255)

                                                   tACGCAAGGTGACGGGTTTATCATTGAACATTCGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTgcagca                                      <  1:4401242/93‑1 (MQ=255)

                                                    aCGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCAt                                     <  2:881097/93‑1 (MQ=255)

                                                    aCGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCAt                                     >  2:555433/1‑93 (MQ=255)

                                                    aCGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCAt                                     <  1:3406821/93‑1 (MQ=255)

                                                       cAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGaa                                  <  1:268616/93‑1 (MQ=255)

                                                       cAAGGTGACGGGCTTATCATTGAACATTTGCAGCAACCCCAGATACTGCGCCAACTGTTCTTCTTCGGACGGGAAACCACTATGAATCACGaa                                  >  2:1572820/1‑93 (MQ=255)

                                                        aaGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAAt                                 >  2:2215386/1‑93 (MQ=255)

                                                         aGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATg                                >  2:3748414/1‑93 (MQ=255)

                                                         aGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATg                                <  2:3239288/93‑1 (MQ=255)

                                                           gTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATggg                              <  1:3594792/93‑1 (MQ=255)

                                                            tGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGa                             <  2:1722860/93‑1 (MQ=255)

                                                             gACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGAt                            <  2:332597/93‑1 (MQ=255)

                                                              aCGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGAtt                           <  2:781018/93‑1 (MQ=255)

                                                                gggTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTc                         <  2:616006/93‑1 (MQ=255)

                                                                 gttttATCATTGAACTTTTGCAGCTCCCCCTGATACGTCCCCACCTGTTCTTCTTCCGACGGGAAACCACTTGGCAGCATGAATGGGATTTCa                        <  2:1153366/91‑1 (MQ=255)

                                                                  gTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCCCTTTGCAGCATGAATGGGATTTCAg                       <  2:4052307/93‑1 (MQ=255)

                                                                  gTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAg                       <  2:4326455/93‑1 (MQ=255)

                                                                  gTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAg                       <  1:1970207/93‑1 (MQ=255)

                                                                   tttATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGt                      >  2:2141627/1‑93 (MQ=255)

                                                                    ttATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTg                     >  1:685716/1‑93 (MQ=255)

                                                                       tCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGa                  >  1:3446464/1‑93 (MQ=255)

                                                                        cATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGAt                 <  1:2175647/93‑1 (MQ=255)

                                                                        cATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGAt                 >  2:351070/1‑93 (MQ=255)

                                                                           tGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATaaa              <  1:2892495/93‑1 (MQ=255)

                                                                           tGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATaaa              <  1:1291013/93‑1 (MQ=255)

                                                                            gAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAg             <  1:3234229/93‑1 (MQ=255)

                                                                            gAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAg             >  2:1836073/1‑93 (MQ=255)

                                                                            gAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAg             <  2:217114/93‑1 (MQ=255)

                                                                             aaCATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGa            >  2:2415852/1‑93 (MQ=255)

                                                                             aaCATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGa            >  1:1593484/1‑93 (MQ=255)

                                                                             aaCATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGa            >  2:3912731/1‑93 (MQ=255)

                                                                             aaCATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGa            <  1:1251/93‑1 (MQ=255)

                                                                              aCATTTGCAGCCCCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGAc           <  2:3365538/93‑1 (MQ=255)

                                                                              aCATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGAc           <  2:2037930/93‑1 (MQ=255)

                                                                              aCATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGAc           <  1:1958352/93‑1 (MQ=255)

                                                                              aCATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGAc           >  1:419804/1‑93 (MQ=255)

                                                                                aTTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCg         >  1:3678591/1‑93 (MQ=255)

                                                                                 tttGCAGCACCCCCTGATGCTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGTGATTTCAGTGCGATAAAGACCGa        <  2:2331618/93‑1 (MQ=255)

                                                                                 tttGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGa        >  1:4040804/1‑93 (MQ=255)

                                                                                 tttGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGa        <  1:337442/93‑1 (MQ=255)

                                                                                  ttGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGAt       >  2:3195984/1‑93 (MQ=255)

                                                                                  ttGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGAt       <  2:4665182/93‑1 (MQ=255)

                                                                                  ttGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGAt       <  2:2011042/93‑1 (MQ=255)

                                                                                   tGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATg      >  1:2487706/1‑93 (MQ=255)

                                                                                      agcaCCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCa   <  1:1341012/93‑1 (MQ=255)

CTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCAT  >  NC_000913/2964911‑2965090

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