BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,284,457	C→T	G165D (GGC→GAC)	acs ←	acetyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,284,457	0	C	T	99.0%	1023.3 / NA	306	G165D (GGC→GAC)	acs	acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand): ref base C (1/0); major base T (151/151); minor base G (3/0); total (155/151)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 345 total aligned reads displayed.

TACTGCGCCCGGCACGCACACCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGCCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCC  >  NC_000913/4284365‑4284546

tACTGCGCCCGGCACGCACACCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAgc                                                                                           >  1:1077726/1‑93 (MQ=255)

   tGCGCCCGGCACGCACACCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAgtcgc                                                                                        <  2:7189021/93‑1 (MQ=255)

    gcgcCCGGCACGCACACCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAgtcgcc                                                                                       >  2:6407121/1‑93 (MQ=255)

    gcgcCCGGCACGCACACCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAgtcgcc                                                                                       >  2:6964156/1‑93 (MQ=255)

       cccGGCACGCACACCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGtcgccgga                                                                                    >  2:1738824/1‑91 (MQ=255)

       cacGGCACGCACACCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGTAGtcgccgga                                                                                    >  2:2984768/3‑91 (MQ=255)

         cGGCACGCACACCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGaaaa                                                                                  <  1:4941947/93‑1 (MQ=255)

         cGGCACGCACACCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGaaaa                                                                                  <  1:2528422/93‑1 (MQ=255)

          ggCACGCACACCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAAt                                                                                 >  1:260581/1‑93 (MQ=255)

          ggCACGCACACCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAAt                                                                                 <  1:522155/93‑1 (MQ=255)

            cacgcacACCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCa                                                                               <  1:2075094/93‑1 (MQ=255)

             acgcacACCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTTCGGCGAGGAGTCGCCGAAAAacac                                                                              >  2:4310580/1‑93 (MQ=255)

              cgcacACCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCAcc                                                                             <  1:1441817/93‑1 (MQ=255)

              cgcacACCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCAcc                                                                             >  1:1868287/1‑93 (MQ=255)

              cgcacACCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCAcc                                                                             <  2:1336816/93‑1 (MQ=255)

                 acacCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGaa                                                                          <  2:3121591/93‑1 (MQ=255)

                  cacCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAAt                                                                         >  2:6520198/1‑93 (MQ=255)

                   acCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATg                                                                        >  2:2292123/1‑93 (MQ=255)

                   acCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATg                                                                        >  2:4670275/1‑93 (MQ=255)

                   acCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATg                                                                        >  1:3239426/1‑93 (MQ=255)

                   acCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATg                                                                        >  1:1781822/1‑93 (MQ=255)

                     cTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCa                                                                      <  1:5978143/93‑1 (MQ=255)

                       tcgtcgGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCAcc                                                                    >  2:6296656/1‑93 (MQ=255)

                       tcgtcgGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCAcc                                                                    <  1:2777002/93‑1 (MQ=255)

                       tcgtcgGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCAcc                                                                    <  1:5848640/93‑1 (MQ=255)

                        cgtcgGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACcg                                                                   >  1:6970104/1‑93 (MQ=255)

                             gAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATAACCGAATGCCCCGCGCCa                                                              >  1:6830910/1‑93 (MQ=255)

                              aaGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCaa                                                             <  2:52616/93‑1 (MQ=255)

                               aGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAAt                                                            <  1:1025272/93‑1 (MQ=255)

                                gTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATg                                                           <  1:3141330/93‑1 (MQ=255)

                                gTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATg                                                           <  1:5688719/93‑1 (MQ=255)

                                   aTCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCgg                                                        <  1:7184198/93‑1 (MQ=255)

                                   aTCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCgg                                                        <  1:3896355/93‑1 (MQ=255)

                                   aTCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTAGCCGAAAATCACCGAATGCACCGCGCCAATGCgg                                                        <  1:5537272/93‑1 (MQ=255)

                                      aCCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGgcg                                                     >  2:3282960/1‑93 (MQ=255)

                                      aCCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGgcg                                                     >  2:7706959/1‑93 (MQ=255)

                                        cAGGCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCa                                                   <  1:7153911/93‑1 (MQ=255)

                                         aGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAg                                                  >  2:3668450/1‑93 (MQ=255)

                                          gTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAgg                                                 >  1:6139805/1‑93 (MQ=255)

                                          gTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAgg                                                 >  2:6182442/1‑93 (MQ=255)

                                           tCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGc                                                >  2:5070114/1‑93 (MQ=255)

                                           tCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAAGc                                                >  2:6067378/1‑93 (MQ=255)

                                            cGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGcc                                               <  2:2244637/93‑1 (MQ=255)

                                            cGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGcc                                               <  2:5528362/93‑1 (MQ=255)

                                            cGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGcc                                               <  1:6769161/93‑1 (MQ=255)

                                             gggAGTTGGAATCAATAAGGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCa                                              <  2:7595395/91‑1 (MQ=255)

                                             gTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCa                                              >  2:5885452/1‑93 (MQ=255)

                                              tGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGGCGCCGAAAAACACCGAATGCACCGCGCCAAAGCGGGCGCAGGCCAg                                             >  1:7046463/1‑93 (MQ=255)

                                               gAGTTGGAATCAATAATGCGCCCGGCAACGGCTTTCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGGCAATGCGGGCGCAGGCCAGc                                            >  1:1570726/1‑93 (MQ=255)

                                               gAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGc                                            <  2:2747502/93‑1 (MQ=255)

                                                 gTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCAt                                          <  2:2648079/93‑1 (MQ=255)

                                                     gAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCa                                      <  2:4849317/93‑1 (MQ=255)

                                                     gAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCa                                      <  1:2277939/93‑1 (MQ=255)

                                                       aTCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAAc                                    <  1:753614/93‑1 (MQ=255)

                                                       aTCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAAc                                    <  1:1257811/93‑1 (MQ=255)

                                                        tCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAAcc                                   <  1:4585923/93‑1 (MQ=255)

                                                        tCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAAcc                                   >  1:3615727/1‑93 (MQ=255)

                                                         cAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACcg                                  <  1:6204402/93‑1 (MQ=255)

                                                           ataatGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACcgcg                                >  2:117919/1‑93 (MQ=255)

                                                            taatGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCgg                               >  1:3419889/1‑93 (MQ=255)

                                                             aatGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGc                              <  2:1979990/93‑1 (MQ=255)

                                                              atGCGCCCGGCAACGGCTTCGGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCt                             <  2:4091738/93‑1 (MQ=255)

                                                              atGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCt                             <  1:4783533/93‑1 (MQ=255)

                                                               tGCGCCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCtt                            >  2:2677457/1‑93 (MQ=255)

                                                                gcgcCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTc                           <  2:1937467/93‑1 (MQ=255)

                                                                gcgcCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTc                           <  2:4352725/93‑1 (MQ=255)

                                                                 ggcCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTcc                          <  1:4674971/92‑1 (MQ=255)

                                                                 ggcCCGGCAACGGCTCCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTcc                          <  1:4242535/92‑1 (MQ=255)

                                                                 cgcCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCGCCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTcc                          <  1:2950151/93‑1 (MQ=255)

                                                                 cgcCCGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTcc                          <  1:2501803/93‑1 (MQ=255)

                                                                    ccGGCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGc                       >  2:5665235/1‑93 (MQ=255)

                                                                       gCAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCAcc                    >  2:6603256/1‑93 (MQ=255)

                                                                        cAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCa                   <  1:7355319/93‑1 (MQ=255)

                                                                        cAACGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCa                   <  2:163670/93‑1 (MQ=255)

                                                                          aCGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCTCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATc                 <  1:6304864/93‑1 (MQ=255)

                                                                          aCGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATc                 >  2:4575064/1‑93 (MQ=255)

                                                                          aCGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATc                 >  2:3510400/1‑93 (MQ=255)

                                                                          aCGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCCTCCTGCACCATc                 >  2:7128400/1‑93 (MQ=255)

                                                                           cGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCg                <  1:6402754/93‑1 (MQ=255)

                                                                           cGGCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCg                <  2:1492733/93‑1 (MQ=255)

                                                                            ggCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCgg               >  1:4507186/1‑93 (MQ=255)

                                                                            ggCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCgg               >  1:68712/1‑93 (MQ=255)

                                                                            ggCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCgg               <  2:7922334/93‑1 (MQ=255)

                                                                             tcTTCCGCCGAGAGGTCGCCGAAAGTCTCCGAATCCCCCGCGCCAATGAGAGCGCAGGCCAGCTGCGCAACCACGGCTTCCGGCACCCTCggt              <  2:2478182/92‑2 (MQ=255)

                                                                             gCTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGc              >  2:899411/1‑93 (MQ=255)

                                                                              cTTCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCa             <  1:4060952/93‑1 (MQ=255)

                                                                               ttCCGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCat            >  2:6733018/1‑93 (MQ=255)

                                                                               tcccGGCGAGAAGTCGCCAAAAATCACCGATTGCCCCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCat            <  2:5737921/91‑1 (MQ=255)

                                                                                gccGGCGAGAAGTCGCCTAAAATCCCCGAATGCACCGCGCGAAGGCGGGCGCAGGCCAGCATCGCACCCGCGGCTTCCGGCACCATCGGCata           <  2:7398522/92‑1 (MQ=255)

                                                                                  cGGCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCAtata         <  1:3514555/93‑1 (MQ=255)

                                                                                   ggCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATaa        <  1:3252388/93‑1 (MQ=255)

                                                                                    gCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATaaa       >  1:4327481/1‑93 (MQ=255)

                                                                                    gCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATaaa       >  1:3954047/1‑93 (MQ=255)

                                                                                    gCGAGAAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATaaa       <  2:445482/93‑1 (MQ=255)

                                                                                      gagaAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATc     >  2:4762273/1‑93 (MQ=255)

                                                                                      gagaAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATc     <  2:3827781/93‑1 (MQ=255)

                                                                                      gagaAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATc     >  1:892303/1‑93 (MQ=255)

                                                                                      gagaAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATc     >  2:1394758/1‑93 (MQ=255)

                                                                                        gaAGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGc   >  1:5167350/1‑93 (MQ=255)

                                                                                         aaGTCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCTGCACCATCGGCATATAAATCGcc  >  2:7374688/1‑93 (MQ=255)

TACTGCGCCCGGCACGCACACCTTCGTCGGAAGTGATCACCAGTCGTGAGTTGGAATCAATAATGCGCCCGGCAACGGCTTCCGGCGAGAAGCCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCC  >  NC_000913/4284365‑4284546

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