BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,923,951	T→C	H194H (CAT→CAC)	yqcD →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,923,951	0	T	C	84.9%	623.4 / 104.5	271	H194H (CAT→CAC)	yqcD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (21/20); new base C (114/116); total (135/136)
Fisher's exact test for biased strand distribution p-value = 8.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

Only 100 of 281 total aligned reads displayed.

ACTGACTATCTGGAGAATGCCACCTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCA  >  NC_000913/2923859‑2924043

actGACTATCTGGAGAATGCCACCTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCAt                                                                                              >  2:4738731/1‑93 (MQ=255)

  tGACTATCTGGAGAATGCCACCTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACccacca                                                                                            >  1:1086724/1‑93 (MQ=255)

  tGACTATCTGGAGAATGCCACCTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACccacca                                                                                            >  2:1850428/1‑93 (MQ=255)

      tATCTGGAGAATGCCACCTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCaccaacca                                                                                        <  1:5850437/93‑1 (MQ=255)

      tATCTGGAGAATGCCACCTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCa                                                                                        >  2:5603924/1‑93 (MQ=255)

        tCTGGAGAATGCCACCTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGa                                                                                      >  1:668373/1‑93 (MQ=255)

         cTGGAGAATGCCACCTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGAt                                                                                     >  1:4285695/1‑93 (MQ=255)

         cTGGAGAATGCCACCTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGAt                                                                                     <  1:3416341/93‑1 (MQ=255)

         cTGGAGAATGCCACCTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGAt                                                                                     <  1:260698/93‑1 (MQ=255)

          tGGAGAATGCCACCTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGAtt                                                                                    <  1:2901546/93‑1 (MQ=255)

          tGGAGAATGCCACCTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGAtt                                                                                    >  2:3171673/1‑93 (MQ=255)

           ggAGAATGCCACCTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTg                                                                                   >  1:2120795/1‑93 (MQ=255)

            gagaATGCCACCTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTgg                                                                                  <  1:3095069/93‑1 (MQ=255)

             agaATGCCACCTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTggg                                                                                 <  1:7242759/93‑1 (MQ=255)

                aTGCCACCTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGtt                                                                              <  2:5503908/93‑1 (MQ=255)

                 tGCCACCCGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTc                                                                             <  1:666415/93‑1 (MQ=255)

                    cACCTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCt                                                                          >  1:2055935/1‑93 (MQ=255)

                     aCCTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTc                                                                         <  2:4707843/93‑1 (MQ=255)

                     aCCTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTc                                                                         >  1:4814285/1‑93 (MQ=255)

                      ccTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTcc                                                                        <  1:2535001/93‑1 (MQ=255)

                      ccTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTcc                                                                        <  2:7109761/93‑1 (MQ=255)

                      ccTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTcc                                                                        >  1:4417688/1‑93 (MQ=255)

                       cTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCa                                                                       >  2:2567598/1‑93 (MQ=255)

                         gtggtgAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCaaa                                                                     >  2:4992763/1‑93 (MQ=255)

                            gtgAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTc                                                                  >  2:6676186/1‑93 (MQ=255)

                                aaaaGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTa                                                              >  1:5928893/1‑93 (MQ=255)

                                aaaaGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTa                                                              <  2:2891976/93‑1 (MQ=255)

                                aaaaGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTa                                                              <  1:1051087/93‑1 (MQ=255)

                                 aaaGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTAt                                                             <  1:4939143/93‑1 (MQ=255)

                                 aaaGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTAt                                                             >  1:90761/1‑93 (MQ=255)

                                 aaaGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTAt                                                             <  1:4369969/93‑1 (MQ=255)

                                 aaaGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTAt                                                             <  1:1652552/93‑1 (MQ=255)

                                 aaaGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTAt                                                             >  2:5259439/1‑93 (MQ=255)

                                 aaaGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTAt                                                             <  2:5793564/93‑1 (MQ=255)

                                  aaGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATc                                                            <  2:366803/93‑1 (MQ=255)

                                  aaGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATc                                                            <  1:4230433/93‑1 (MQ=255)

                                  aaGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATc                                                            >  1:3825595/1‑93 (MQ=255)

                                  aaGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATc                                                            <  2:4120318/93‑1 (MQ=255)

                                    gtagtGGAAGAGACGCTTGTCAGCCCCCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGt                                                          >  1:349378/1‑93 (MQ=255)

                                    gtagtGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGt                                                          >  1:1657777/1‑93 (MQ=255)

                                     tagtGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTg                                                         <  1:4147041/93‑1 (MQ=255)

                                     tagtGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTg                                                         <  1:402931/93‑1 (MQ=255)

                                       gtGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGa                                                       >  1:586592/1‑93 (MQ=255)

                                       gtGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGa                                                       <  1:3636261/93‑1 (MQ=255)

                                        tGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGAc                                                      >  2:3286918/1‑93 (MQ=255)

                                          gAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGc                                                    <  2:3073115/93‑1 (MQ=255)

                                          gAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGc                                                    >  2:3512137/1‑93 (MQ=255)

                                          gAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGc                                                    >  1:1265796/1‑93 (MQ=255)

                                           aaGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGcc                                                   <  2:1147980/93‑1 (MQ=255)

                                           aaGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGcc                                                   <  2:1707824/93‑1 (MQ=255)

                                            agagACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCa                                                  >  1:3550479/1‑93 (MQ=255)

                                             gagaCGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCaa                                                 >  2:6253802/1‑93 (MQ=255)

                                               gaCGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAAt                                               >  2:1775868/1‑93 (MQ=255)

                                                aCGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAAtt                                              >  1:2732892/1‑93 (MQ=255)

                                                aCGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAAtt                                              <  2:5085388/93‑1 (MQ=255)

                                                 cGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTg                                             >  1:835169/1‑93 (MQ=255)

                                                  gCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGa                                            <  1:5189850/93‑1 (MQ=255)

                                                  gCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGa                                            <  2:5346991/93‑1 (MQ=255)

                                                   cTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGAc                                           >  1:3910631/1‑93 (MQ=255)

                                                   cTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGAc                                           <  1:6104185/93‑1 (MQ=255)

                                                    ttGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACa                                          >  1:654697/1‑93 (MQ=255)

                                                    ttGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACa                                          >  2:5907337/1‑93 (MQ=255)

                                                    ttGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACa                                          >  2:7443532/1‑93 (MQ=255)

                                                      gTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACaga                                        >  1:534232/1‑93 (MQ=255)

                                                      gTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACaga                                        <  2:3825539/93‑1 (MQ=255)

                                                      gTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACaga                                        <  2:726241/93‑1 (MQ=255)

                                                            cACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGaaaaa                                  >  1:1750436/1‑93 (MQ=255)

                                                            cACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGaaaaa                                  <  1:4778634/93‑1 (MQ=255)

                                                            cACCTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGaaaaa                                  <  2:4507731/93‑1 (MQ=255)

                                                             aCCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAAc                                 >  2:4389217/1‑93 (MQ=255)

                                                              ccTGCTGAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAAct                                <  2:1380995/93‑1 (MQ=255)

                                                               ctgctgAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAActg                               <  2:1189646/93‑1 (MQ=255)

                                                               ctgctgAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAActg                               >  2:1025250/1‑93 (MQ=255)

                                                                tgctgAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTTAGTATCGTGGACGCCAAATTGACAGAGAAAAActgc                              >  2:6863822/1‑93 (MQ=255)

                                                                tgctgAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAActgc                              <  2:2389174/93‑1 (MQ=255)

                                                                tgctgAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAActgc                              <  1:2281638/93‑1 (MQ=255)

                                                                tgctgAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAActgc                              <  1:6473007/93‑1 (MQ=255)

                                                                tgctgAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAActgc                              >  1:756681/1‑93 (MQ=255)

                                                                 gctgAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAActgct                             >  1:223886/1‑93 (MQ=255)

                                                                 gctgAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAActgct                             >  1:559402/1‑93 (MQ=255)

                                                                 gctgAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAActgct                             >  1:5794087/1‑93 (MQ=255)

                                                                 gctgAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAActgct                             <  1:1463141/93‑1 (MQ=255)

                                                                  ctgAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAActgctg                            <  1:741192/93‑1 (MQ=255)

                                                                   tgAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACtgctgc                           <  1:5416300/93‑1 (MQ=255)

                                                                   tgAAATCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACtgctgc                           <  1:6970519/93‑1 (MQ=255)

                                                                     aaaTCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGt                         <  2:5425227/93‑1 (MQ=255)

                                                                      aaTCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGtt                        <  2:2341398/93‑1 (MQ=255)

                                                                      aaTCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGtt                        >  1:1222747/1‑93 (MQ=255)

                                                                      aaTCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGtt                        <  1:777692/93‑1 (MQ=255)

                                                                       aTCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTa                       >  1:2963940/1‑93 (MQ=255)

                                                                        tCAAACTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTAc                      >  2:4250734/1‑93 (MQ=255)

                                                                          aaaCTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCt                    >  1:6622418/1‑93 (MQ=255)

                                                                          aaaCTGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCt                    >  1:3688841/1‑93 (MQ=255)

                                                                              tGCCTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGtc                <  2:899871/93‑1 (MQ=255)

                                                                                 cTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCa             <  2:6036332/93‑1 (MQ=255)

                                                                                 cTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCa             <  2:5724371/93‑1 (MQ=255)

                                                                                 cTGATCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCa             <  2:5040622/93‑1 (MQ=255)

                                                                                  tGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCAt            >  1:6037538/1‑93 (MQ=255)

                                                                                    aTCACCCACCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTc          >  2:6386389/1‑93 (MQ=255)

                                                                                         ccATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGtca     >  2:3764672/1‑93 (MQ=255)

ACTGACTATCTGGAGAATGCCACCTGTGGTGAAAAAGTAGTGGAAGAGACGCTTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCA  >  NC_000913/2923859‑2924043

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