BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,787,263	C→G	G149A (GGC→GCC)	yibD ←	putative regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,787,263	0	C	G	99.5%	1342.7 / ‑3.4	432	G149A (GGC→GCC)	yibD	putative regulator
Reads supporting (aligned to +/- strand): ref base C (2/0); new base G (228/200); total (232/200)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 450 total aligned reads displayed.

TTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCAC  >  NC_000913/3787173‑3787352

tttttAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCg                                                                                         >  2:2393403/1‑93 (MQ=255)

tttttAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCg                                                                                         >  2:1421955/1‑93 (MQ=255)

 ttttAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCgg                                                                                        <  2:602905/93‑1 (MQ=255)

 ttttAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCgg                                                                                        >  2:3063356/1‑93 (MQ=255)

   ttAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGtt                                                                                      >  2:7307044/1‑93 (MQ=255)

       cAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATa                                                                                  >  2:2088752/1‑93 (MQ=255)

       cAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATa                                                                                  >  2:4161720/1‑93 (MQ=255)

        aaTAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATAc                                                                                 >  1:466402/1‑93 (MQ=255)

         aTAACATCATGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACg                                                                                >  2:1335559/1‑93 (MQ=255)

         aTAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACg                                                                                >  1:5586026/1‑93 (MQ=255)

         aTAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACg                                                                                <  1:3543468/93‑1 (MQ=255)

         aTAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACg                                                                                >  2:701603/1‑93 (MQ=255)

          tAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGc                                                                               <  1:190761/93‑1 (MQ=255)

           aaCATCACGGCGATAAACCGCCATCCAGACAACGGGAGTCCAACGGCGCGAAGAAATCCCCATCCGCAGCCAGTCCCGGGTGGTTAATACGcc                                                                              >  2:4722161/1‑93 (MQ=255)

           aaCATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGcc                                                                              >  1:2597481/1‑93 (MQ=255)

           aaCATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGcc                                                                              >  1:5066348/1‑93 (MQ=255)

           aaCATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGcc                                                                              >  1:5162363/1‑93 (MQ=255)

             cATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCgg                                                                            >  2:7830136/1‑93 (MQ=255)

              aTCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGt                                                                           <  2:2703648/93‑1 (MQ=255)

              aTCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGt                                                                           <  1:3008314/93‑1 (MQ=255)

              aTCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGt                                                                           <  1:6508475/93‑1 (MQ=255)

              aTCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGt                                                                           >  2:5405543/1‑93 (MQ=255)

               tCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGtt                                                                          <  1:4328119/93‑1 (MQ=255)

                 aCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCAGGGCGGTTAATACGCCGGTTGa                                                                        <  2:4982200/93‑1 (MQ=255)

                  cGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAg                                                                       >  1:7385939/1‑93 (MQ=255)

                  cGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAg                                                                       >  1:5461192/1‑93 (MQ=255)

                   ggCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGc                                                                      >  1:6229719/1‑93 (MQ=255)

                    gCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCg                                                                     >  2:7682138/1‑93 (MQ=255)

                     ggaaaaaCCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGa                                                                    <  1:4399488/89‑1 (MQ=255)

                       aTAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAg                                                                  >  2:3550656/1‑93 (MQ=255)

                        tAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAgg                                                                 >  2:1240672/1‑93 (MQ=255)

                          aaCCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCg                                                               >  1:1163363/1‑93 (MQ=255)

                           aCCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGa                                                              <  1:4531450/93‑1 (MQ=255)

                            cccccATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGac                                                             >  1:7293232/1‑92 (MQ=255)

                            cccccATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGAt                                                             <  2:6785332/93‑1 (MQ=255)

                            cccccATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGAt                                                             <  1:4883847/93‑1 (MQ=255)

                             ccccATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATc                                                            >  2:1140180/1‑93 (MQ=255)

                              cccATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCg                                                           <  1:3172245/93‑1 (MQ=255)

                              cccATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCg                                                           >  2:2929259/1‑93 (MQ=255)

                              cccATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCg                                                           >  1:3840707/1‑93 (MQ=255)

                              cccATCCAGACAACGTGAGTCCAAAGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCg                                                           >  2:7586340/1‑93 (MQ=255)

                               ccATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCgg                                                          <  1:1261654/93‑1 (MQ=255)

                                cATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGt                                                         >  1:5789364/1‑93 (MQ=255)

                                cATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGt                                                         <  2:897307/93‑1 (MQ=255)

                                cATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGt                                                         <  1:1667884/93‑1 (MQ=255)

                                cATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGt                                                         <  1:7171541/93‑1 (MQ=255)

                                 aTCCAGACAGCGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTg                                                        <  2:4752008/93‑1 (MQ=255)

                                 aTCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTg                                                        <  2:4427246/93‑1 (MQ=255)

                                  tCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTgg                                                       >  2:1010826/1‑93 (MQ=255)

                                  tCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTgg                                                       >  2:5138368/1‑93 (MQ=255)

                                     aGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGgga                                                    <  2:2619782/93‑1 (MQ=255)

                                      gACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGggat                                                   <  2:5667095/93‑1 (MQ=255)

                                      gACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGggat                                                   >  2:6883181/1‑93 (MQ=255)

                                         aaCGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGggatgga                                                <  1:2185319/93‑1 (MQ=255)

                                           cGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGAtt                                              >  1:6837176/1‑93 (MQ=255)

                                             tGAGTCCAACGGCGCGAACAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGCGCGGTGGGGATGGATTGc                                            >  2:65961/1‑93 (MQ=255)

                                              gAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGcc                                           <  2:4033700/93‑1 (MQ=255)

                                               aGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCa                                          >  1:406730/1‑93 (MQ=255)

                                                gTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAg                                         >  1:2857619/1‑93 (MQ=255)

                                                 tCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAgg                                        >  2:6486624/1‑93 (MQ=255)

                                                  ccAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGt                                       >  2:5780643/1‑93 (MQ=255)

                                                  ccAACGGCGAGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGt                                       <  2:1918127/93‑1 (MQ=255)

                                                   cAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGtt                                      >  2:3367605/1‑93 (MQ=255)

                                                   cAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGtt                                      <  2:3499146/93‑1 (MQ=255)

                                                    aaCGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGttt                                     >  2:1955543/1‑93 (MQ=255)

                                                     aCGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTtc                                    <  1:1302738/93‑1 (MQ=255)

                                                       ggCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTtctc                                  <  1:7668782/93‑1 (MQ=255)

                                                           cgAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGt                              <  2:7459980/93‑1 (MQ=255)

                                                            gaagaaAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGtt                             >  1:6571955/1‑93 (MQ=255)

                                                             aagaaAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGttt                            <  1:5992487/93‑1 (MQ=255)

                                                              agaaAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTc                           >  1:3320134/1‑93 (MQ=255)

                                                               gaaAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCCGGTTTCTCCCGTTTCa                          >  2:2216978/1‑93 (MQ=255)

                                                               gaaAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCa                          >  1:1042729/1‑93 (MQ=255)

                                                               aaaaGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCa                          <  1:6198424/92‑1 (MQ=255)

                                                                aaaGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCAc                         >  2:2547538/1‑93 (MQ=255)

                                                                  aGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGa                       >  1:7747670/1‑93 (MQ=255)

                                                                   gCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGaa                      >  1:2681873/1‑93 (MQ=255)

                                                                    ccccATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGaaa                     >  1:1836619/1‑93 (MQ=255)

                                                                    ccccATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGaaa                     >  2:3667629/1‑93 (MQ=255)

                                                                       cATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGCAAaca                  >  2:5709575/1‑93 (MQ=255)

                                                                        aTCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAacac                 >  1:153967/1‑93 (MQ=255)

                                                                         tCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACAcc                >  1:6141297/1‑93 (MQ=255)

                                                                         tCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACAcc                >  2:1665274/1‑93 (MQ=255)

                                                                         tCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACAcc                >  2:321517/1‑93 (MQ=255)

                                                                          ccGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCa               <  1:7843735/93‑1 (MQ=255)

                                                                            gCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGt             <  2:6520974/93‑1 (MQ=255)

                                                                             cAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTc            <  1:69579/93‑1 (MQ=255)

                                                                             cAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTc            <  2:3878006/93‑1 (MQ=255)

                                                                               gCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAg          >  1:3618932/1‑93 (MQ=255)

                                                                               gCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAg          <  1:3389171/93‑1 (MQ=255)

                                                                               gCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAg          >  2:6147781/1‑93 (MQ=255)

                                                                               gCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAg          <  1:2340210/93‑1 (MQ=255)

                                                                                 cAGTCCGGGGCGTTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCg        <  2:4547578/93‑1 (MQ=255)

                                                                                  aGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGt       >  2:5260984/1‑93 (MQ=255)

                                                                                   gTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGtt      <  1:5307427/93‑1 (MQ=255)

                                                                                   gTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGtt      >  1:4161138/1‑93 (MQ=255)

                                                                                    tCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTg     <  2:3925723/93‑1 (MQ=255)

                                                                                      cGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCa   <  2:3829217/93‑1 (MQ=255)

                                                                                      cGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCa   <  2:4899563/93‑1 (MQ=255)

                                                                                       ggggCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCAc  <  1:5382672/93‑1 (MQ=255)

TTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCAC  >  NC_000913/3787173‑3787352

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: `‑`