BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,754,634	T→C	N274S (AAC→AGC)	yiaY ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,754,634	0	T	C	87.4%	718.7 / 75.9	286	N274S (AAC→AGC)	yiaY	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (22/14); new base C (119/130); total (141/145)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 289 total aligned reads displayed.

GGCACAGTCACGCAGACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGTTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCGAGGAACTGGGCATAA  >  NC_000913/3754543‑3754723

ggCACAGTCACGCAGACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGtt                                                                                          >  1:1065893/1‑93 (MQ=255)

 gCACAGTCACGCAGACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTg                                                                                         >  1:4809434/1‑93 (MQ=255)

   acaGTCACGCAGACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTa                                                                                       >  1:6639146/1‑93 (MQ=255)

   acaGTCACGCAGACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTa                                                                                       <  2:4840615/93‑1 (MQ=255)

     aGTCACGCAGACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGTTGTAGa                                                                                     >  2:3213317/1‑93 (MQ=255)

     aGTCACGCAGACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGa                                                                                     >  1:6589839/1‑93 (MQ=255)

      gTCACGCAGACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGTTGTAGaa                                                                                    >  2:4392206/1‑93 (MQ=255)

      gTCACGCAGACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGaa                                                                                    >  2:1478173/1‑93 (MQ=255)

       tCACGCAGACGTGCAGCGGCGACTTTGCTGTTGAATGCCTGAACGTGCGGCAGCAAAACGGCGTTACATTCACCGTGTGGCAGGCTGTAGaaa                                                                                   <  2:2192201/93‑1 (MQ=255)

       tCACGCAGACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGTTGTAGaaa                                                                                   <  1:3523960/93‑1 (MQ=255)

       tCACGCAGACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGaaa                                                                                   <  2:5743790/93‑1 (MQ=255)

       tCACGCAGACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGaaa                                                                                   <  2:7211505/93‑1 (MQ=255)

        cACGCAGACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGTTGTAGAAAc                                                                                  <  1:5191370/93‑1 (MQ=255)

        cACGCAGACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAAc                                                                                  <  1:1110628/93‑1 (MQ=255)

          cGCAGACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCg                                                                                <  1:3636809/93‑1 (MQ=255)

            cAGACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGcc                                                                              >  1:3282073/1‑93 (MQ=255)

             aGACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGccc                                                                             >  1:1463751/1‑93 (MQ=255)

             aGACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGccc                                                                             <  2:6619305/93‑1 (MQ=255)

             aGACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGccc                                                                             <  1:549742/93‑1 (MQ=255)

              gACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCa                                                                            <  1:4340437/93‑1 (MQ=255)

              gACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCa                                                                            <  2:4137435/93‑1 (MQ=255)

               aCGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAg                                                                           >  1:6720817/1‑93 (MQ=255)

                cGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGc                                                                          <  1:2312682/93‑1 (MQ=255)

                  tGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTg                                                                        >  1:5095325/1‑93 (MQ=255)

                   gCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTgg                                                                       <  1:4554030/93‑1 (MQ=255)

                   gCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTgg                                                                       <  1:2495604/93‑1 (MQ=255)

                    cAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGt                                                                      >  2:469093/1‑93 (MQ=255)

                      gcggcgACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTgc                                                                    <  2:264609/93‑1 (MQ=255)

                      gcggcgACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTgc                                                                    >  1:1811362/1‑93 (MQ=255)

                       cggcgACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTgcg                                                                   >  2:4323651/1‑93 (MQ=255)

                       cggcgACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTgcg                                                                   >  2:3156736/1‑93 (MQ=255)

                        ggcgACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTgcgc                                                                  <  1:4241791/93‑1 (MQ=255)

                           gACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCAt                                                               >  1:4988642/1‑93 (MQ=255)

                           gACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCAt                                                               <  2:3492366/93‑1 (MQ=255)

                            aCTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATc                                                              >  2:2968200/1‑93 (MQ=255)

                            aCTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATc                                                              >  2:2074490/1‑93 (MQ=255)

                            aCTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATc                                                              <  1:1706029/93‑1 (MQ=255)

                             cTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCg                                                             >  2:1694803/1‑93 (MQ=255)

                              tttGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGc                                                            >  1:6877977/1‑93 (MQ=255)

                               ttGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGTTGTAGAAACCGCCCAGCTGGTGCGCCATCGCa                                                           >  2:727247/1‑93 (MQ=255)

                                tGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCAt                                                          <  1:6194079/93‑1 (MQ=255)

                                   tGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGaa                                                       <  2:4501889/93‑1 (MQ=255)

                                   tGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGaa                                                       <  2:1815080/93‑1 (MQ=255)

                                   tGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGaa                                                       >  1:4725154/1‑93 (MQ=255)

                                    gTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAAc                                                      <  2:332777/93‑1 (MQ=255)

                                    gTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAAc                                                      <  2:2310944/93‑1 (MQ=255)

                                     ttGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACa                                                     <  1:7451422/93‑1 (MQ=255)

                                     ttGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACa                                                     <  1:6523953/93‑1 (MQ=255)

                                       gAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGTTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATa                                                   >  2:6059191/1‑93 (MQ=255)

                                          tACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGTTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAAcc                                                >  1:1955497/1‑93 (MQ=255)

                                          tACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGTTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAAcc                                                >  2:4466902/1‑93 (MQ=255)

                                           aCCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAAccc                                               >  2:36852/1‑93 (MQ=255)

                                            ccTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACTTAACCCa                                              >  1:4824985/1‑93 (MQ=255)

                                            ccTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCa                                              >  2:6693741/1‑93 (MQ=255)

                                             cTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGTTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCag                                             >  2:3891922/1‑93 (MQ=255)

                                             cTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGTTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCag                                             >  2:2477992/1‑93 (MQ=255)

                                              tGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGTTGTAGAAACCGCCCAGCTGGTGTGCCATCGCATGAACATAACCCaga                                            >  2:4788113/1‑93 (MQ=255)

                                              tGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCaga                                            >  1:1002349/1‑93 (MQ=255)

                                              nnaaCGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCaga                                            <  1:2221570/91‑1 (MQ=255)

                                                 aCGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGTTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGaa                                         >  2:6168990/1‑93 (MQ=255)

                                                   gTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGTTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGc                                       >  1:990890/1‑93 (MQ=255)

                                                    tGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCa                                      <  1:5915315/93‑1 (MQ=255)

                                                    tGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCa                                      <  2:6316649/93‑1 (MQ=255)

                                                     ccGGCAGCAATACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCCGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCat                                     <  1:6221590/92‑1 (MQ=255)

                                                      cGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCatt                                    >  2:1358023/1‑93 (MQ=255)

                                                        gcagcaAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCattat                                  <  1:5246131/93‑1 (MQ=255)

                                                        gcagcaAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCattat                                  <  2:5516265/93‑1 (MQ=255)

                                                        gcagcaAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCattat                                  >  1:4923167/1‑93 (MQ=255)

                                                           gcaAAACGGCGTTACATACACCGTGTGGCAGGTTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGa                               <  1:6391405/93‑1 (MQ=255)

                                                           gcaAAACGGCGTTACATACACCGTGTGGCAGGTTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGa                               <  1:1366874/93‑1 (MQ=255)

                                                           gcaAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGa                               <  2:2821400/93‑1 (MQ=255)

                                                            caAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGaa                              <  1:7195853/93‑1 (MQ=255)

                                                            caAAACGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGaa                              <  1:475096/93‑1 (MQ=255)

                                                                aCGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGcc                          <  2:1971732/93‑1 (MQ=255)

                                                                 cGGCGTTACATACACCGTGTGGCAGGTTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCa                         <  2:5085030/93‑1 (MQ=255)

                                                                 cGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCa                         >  2:792450/1‑93 (MQ=255)

                                                                 cGGCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCa                         >  1:7226211/1‑93 (MQ=255)

                                                                  ggCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCAt                        >  1:5667713/1‑93 (MQ=255)

                                                                   gCGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATa                       >  1:5963299/1‑93 (MQ=255)

                                                                    cGTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATAc                      <  2:1162960/93‑1 (MQ=255)

                                                                     gTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATAcc                     >  1:3385426/1‑93 (MQ=255)

                                                                     gTTACATACACCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATAcc                     <  2:5337522/93‑1 (MQ=255)

                                                                         catacaCCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGca                 >  2:5915587/1‑92 (MQ=255)

                                                                          atacaCCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGcaa                <  2:7151564/93‑3 (MQ=255)

                                                                           tacaCCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAg               <  2:6866458/93‑1 (MQ=255)

                                                                           tacaCCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAg               <  1:4149941/93‑1 (MQ=255)

                                                                           tacaCCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAg               >  2:1914075/1‑93 (MQ=255)

                                                                             cacCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGa             >  1:5066121/1‑93 (MQ=255)

                                                                              acCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGaa            <  2:406222/93‑1 (MQ=255)

                                                                               ccGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAAc           >  2:7428852/1‑93 (MQ=255)

                                                                                cGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACt          <  1:5371816/93‑1 (MQ=255)

                                                                                cGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACt          >  1:878875/1‑93 (MQ=255)

                                                                                cGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACt          <  1:2904104/93‑1 (MQ=255)

                                                                                 gtgtGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTg         >  2:5141190/1‑93 (MQ=255)

                                                                                 gtgtGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTg         <  2:2760009/93‑1 (MQ=255)

                                                                                 gtgtGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTg         >  2:3028500/1‑93 (MQ=255)

                                                                                  tgtgGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTgg        >  2:3931502/1‑93 (MQ=255)

                                                                                      gCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCAt    >  1:3938249/1‑93 (MQ=255)

                                                                                      gCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCAt    <  2:4218337/93‑1 (MQ=255)

                                                                                        aGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATaa  >  1:7067225/1‑93 (MQ=255)

GGCACAGTCACGCAGACGTGCAGCGGCGACTTTGCTGTTGAATACCTGAACGTGCGGCAGCAAAACGGCGTTACATACACCGTGTGGCAGGTTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCGAGGAACTGGGCATAA  >  NC_000913/3754543‑3754723

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