BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,938,374	C→T	D145N (GAC→AAC)	yiePs ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,938,374	0	C	T	97.9%	823.9 / 0.3	239	D145N (GAC→AAC)	yiePs	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (3/2); new base T (131/101); total (135/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-01

Only 100 of 241 total aligned reads displayed.

GCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATAG  >  NC_000913/3938284‑3938462

gcgcCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGc                                                                                        >  2:4375457/1‑93 (MQ=255)

  gcCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTg                                                                                      <  1:4645482/93‑1 (MQ=255)

      gACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTaa                                                                                  >  1:2202259/1‑93 (MQ=255)

      gACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTaa                                                                                  >  1:253726/1‑93 (MQ=255)

       aCAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAAt                                                                                 <  1:1363340/93‑1 (MQ=255)

       aCAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAAt                                                                                 <  1:4408920/93‑1 (MQ=255)

       aCAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAAt                                                                                 <  2:4252318/93‑1 (MQ=255)

       aCAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAAt                                                                                 >  1:895587/1‑93 (MQ=255)

        cAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAAtt                                                                                <  2:1080596/93‑1 (MQ=255)

          aGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGa                                                                              >  2:3739548/1‑93 (MQ=255)

            cTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACg                                                                            <  2:4435378/93‑1 (MQ=255)

            cTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACg                                                                            >  2:2948283/1‑93 (MQ=255)

                aaaTGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaa                                                                        >  1:4666431/1‑93 (MQ=255)

                 aaTGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaag                                                                       >  2:713827/1‑93 (MQ=255)

                  aTGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaagt                                                                      <  2:3397763/93‑1 (MQ=255)

                  aTGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaagt                                                                      >  2:3018154/1‑93 (MQ=255)

                    tccGCGTCTCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaagtaa                                                                    <  1:3957529/92‑1 (MQ=255)

                    gCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaagtaa                                                                    <  1:4589324/93‑1 (MQ=255)

                    gCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaagtaa                                                                    <  2:378695/93‑1 (MQ=255)

                    gCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaagtaa                                                                    <  2:1193637/93‑1 (MQ=255)

                    gCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaagtaa                                                                    >  2:205245/1‑93 (MQ=255)

                    gCCGCATCGTCTTCCCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGTTTTATCACTGTGTTGCTGGTAATTGACGTGaagtaa                                                                    <  1:1653640/93‑1 (MQ=255)

                     ccGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaagtaag                                                                   >  2:4632384/1‑93 (MQ=255)

                      cgcgTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaagtaagt                                                                  >  2:276224/1‑93 (MQ=255)

                       gcgTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAgtg                                                                 <  1:749488/93‑1 (MQ=255)

                       gcgTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAgtg                                                                 <  1:3988442/93‑1 (MQ=255)

                         gTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGtgtg                                                               <  1:3197234/93‑1 (MQ=255)

                          tCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGa                                                              <  2:1862282/93‑1 (MQ=255)

                          tCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGa                                                              <  1:553343/93‑1 (MQ=255)

                          tCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGa                                                              >  1:3880150/1‑93 (MQ=255)

                          tCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCCGGTAATTGACGTGAAGTAAGTGTGa                                                              <  2:144711/93‑1 (MQ=255)

                           cGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGAt                                                             <  2:2559306/93‑1 (MQ=255)

                             ccATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATaa                                                           >  2:1328457/1‑93 (MQ=255)

                              cATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATaaa                                                          >  1:827612/1‑93 (MQ=255)

                              cATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATaaa                                                          <  1:582656/93‑1 (MQ=255)

                              cATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATaaa                                                          <  1:3683044/93‑1 (MQ=255)

                              cATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATaaa                                                          <  1:3612353/93‑1 (MQ=255)

                               aTCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAAc                                                         >  2:152906/1‑93 (MQ=255)

                               aTCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAAc                                                         >  2:3502227/1‑93 (MQ=255)

                               aTCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAAc                                                         >  1:48675/1‑93 (MQ=255)

                                tCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAAcc                                                        <  2:1690660/93‑1 (MQ=255)

                                tCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAAcc                                                        >  2:3071889/1‑93 (MQ=255)

                                tCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAAcc                                                        <  1:1484343/93‑1 (MQ=255)

                                 cGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCg                                                       >  2:3661777/1‑93 (MQ=255)

                                  gCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGa                                                      >  1:1086257/1‑93 (MQ=255)

                                  gCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGa                                                      >  1:4034739/1‑93 (MQ=255)

                                   cTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGaa                                                     <  1:3452813/93‑1 (MQ=255)

                                    tttGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAAt                                                    >  1:4552406/1‑93 (MQ=255)

                                    tttGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAAt                                                    <  2:3619681/93‑1 (MQ=255)

                                       ggATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGa                                                 >  2:86134/1‑93 (MQ=255)

                                       ggATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGa                                                 >  1:3659919/1‑93 (MQ=255)

                                        gATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGaa                                                <  2:1815121/93‑1 (MQ=255)

                                           aaTCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAg                                             >  1:2790335/1‑93 (MQ=255)

                                             tCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCg                                           <  2:1756764/93‑1 (MQ=255)

                                              cGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGa                                          <  2:2860540/93‑1 (MQ=255)

                                               gCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAg                                         <  1:316636/93‑1 (MQ=255)

                                               gCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAg                                         >  1:2568451/1‑93 (MQ=255)

                                                 aTCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGc                                       <  1:4148581/93‑1 (MQ=255)

                                                 aTCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTCCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGc                                       <  2:3630391/93‑1 (MQ=255)

                                                  tCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCg                                      <  2:2248099/93‑1 (MQ=255)

                                                    gACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGaa                                    >  1:3035685/1‑93 (MQ=255)

                                                    gACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGaa                                    >  1:4318819/1‑93 (MQ=255)

                                                     aCAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAAt                                   >  1:3485945/1‑93 (MQ=255)

                                                     aCAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAAt                                   >  2:2115661/1‑93 (MQ=255)

                                                     aCAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAAt                                   <  1:4445529/93‑1 (MQ=255)

                                                       aaTCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGa                                 <  1:3775665/93‑1 (MQ=255)

                                                        aTCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAg                                >  2:2981951/1‑93 (MQ=255)

                                                         tCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAgg                               <  2:219498/93‑1 (MQ=255)

                                                           gCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTc                             <  1:771884/93‑1 (MQ=255)

                                                            ccTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCa                            <  1:1437336/93‑1 (MQ=255)

                                                             cagATGCAGATCCAGCTTTATCACTGTGTGGCTGGTAATTGACGTGATGTAAGTGTGATTAACCGAATTTAATAGCGAGGCGAATGAGGTCaa                           >  1:4216875/3‑93 (MQ=255)

                                                             cTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCaa                           >  2:221238/1‑93 (MQ=255)

                                                             cTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCaa                           >  2:1280185/1‑93 (MQ=255)

                                                                aTGCAGATCCAGCTTTATCACTGTGTTGCTGGTAAGTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCaaaaa                        <  1:1068424/93‑1 (MQ=255)

                                                                 tGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAAc                       <  2:586922/93‑1 (MQ=255)

                                                                 tGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAAc                       >  2:501236/1‑93 (MQ=255)

                                                                  gCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACg                      <  2:3811817/93‑1 (MQ=255)

                                                                   cAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACgg                     <  2:472991/93‑1 (MQ=255)

                                                                    aGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGa                    >  1:882928/1‑93 (MQ=255)

                                                                     gATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGAt                   >  1:1911549/1‑93 (MQ=255)

                                                                     gATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGGAAGTGTGATAAACCGAATGGCATAGCGAGGCGAATTGGATCAAAAACGGAt                   >  1:977969/1‑93 (MQ=255)

                                                                      aTCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGAtt                  <  2:2631384/93‑1 (MQ=255)

                                                                      aTCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGAtt                  <  2:4670009/93‑1 (MQ=255)

                                                                        ccAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTgg                >  1:2294925/1‑93 (MQ=255)

                                                                        ccAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTgg                >  2:4138933/1‑93 (MQ=255)

                                                                         cAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGgc               <  1:290345/93‑1 (MQ=255)

                                                                           gCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAACAACGGATTGgcgc             >  1:2624812/1‑93 (MQ=255)

                                                                            cTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCt            >  1:2038900/1‑93 (MQ=255)

                                                                             tttATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTc           <  1:2391253/93‑1 (MQ=255)

                                                                             tttATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTc           >  2:1936405/1‑93 (MQ=255)

                                                                              ttATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTct          >  1:1597604/1‑92 (MQ=255)

                                                                              ttATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCa          >  1:279997/1‑93 (MQ=255)

                                                                               tATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCAt         <  1:2119813/93‑1 (MQ=255)

                                                                                aTCGCTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCAtt        <  2:3431757/93‑1 (MQ=255)

                                                                                aTCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCAtt        >  1:2257992/1‑93 (MQ=255)

                                                                                 tCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCAttt       >  2:4068781/1‑93 (MQ=255)

                                                                                 tCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCAttt       <  2:2732718/93‑1 (MQ=255)

                                                                                   aCTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCa     >  1:3265976/1‑93 (MQ=255)

                                                                                   aCTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCa     <  1:1116060/93‑1 (MQ=255)

                                                                                    cTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCAt    >  2:2360278/1‑93 (MQ=255)

GCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATAG  >  NC_000913/3938284‑3938462

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