BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,405,126	T→C	M252V (ATG→GTG)	b1341 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,405,126	0	T	C	98.3%	546.7 / ‑0.9	173	M252V (ATG→GTG)	b1341	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (3/0); new base C (87/83); total (90/83)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.47e-01

Only 100 of 185 total aligned reads displayed.

GTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAG  >  NC_000913/1405036‑1405215

gTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGa                                                                                         >  1:2504128/1‑93 (MQ=255)

gTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGa                                                                                         >  1:1988185/1‑93 (MQ=255)

gTAATCCGGTCATCGCGTCATGGCGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGa                                                                                         <  2:2579167/93‑1 (MQ=255)

   aTCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCt                                                                                      <  1:3358384/93‑1 (MQ=255)

       ggTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTc                                                                                  <  1:3164105/93‑1 (MQ=255)

         tCATCGCGTCATGGTGAGCAGCCTGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTACTATCATGCACAATACATCACACGAGCTTGTCGc                                                                                >  2:1869450/1‑93 (MQ=255)

          cATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGcc                                                                               <  1:1145819/93‑1 (MQ=255)

           agcgcgTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCa                                                                              <  1:4468492/91‑1 (MQ=255)

           aTCGCGTCTTGGTGAGCAGCATGTCCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCa                                                                              <  2:4526236/93‑1 (MQ=255)

              gcgTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATaa                                                                           >  1:3196920/1‑93 (MQ=255)

               cgTCATGGTGAGCAGCATGTCCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATaaa                                                                          <  1:1178191/93‑1 (MQ=255)

                gTCATGGTGAGCAGCATGTT‑CAGCTGCTCCTCCAGCCGGTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAAtt                                                                        >  2:3848631/1‑93 (MQ=255)

                 tCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAAtt                                                                        >  1:1556649/1‑93 (MQ=255)

                 tCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAAtt                                                                        >  1:4315976/1‑93 (MQ=255)

                 tCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAAtt                                                                        <  2:3687240/93‑1 (MQ=255)

                  cATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAAttt                                                                       <  1:2015078/93‑1 (MQ=255)

                     ggTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCaa                                                                    <  2:3324294/93‑1 (MQ=255)

                      gTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCaac                                                                   >  2:2427770/1‑92 (MQ=255)

                      gTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAAt                                                                   >  2:2427768/1‑93 (MQ=255)

                        gAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCg                                                                 <  1:4157978/93‑1 (MQ=255)

                          gcagcgTGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGt                                                               >  2:2420747/1‑93 (MQ=255)

                           cagcaTGTACCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTc                                                              <  2:2086065/93‑1 (MQ=255)

                            agcaTGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTcc                                                             <  1:4364269/93‑1 (MQ=255)

                            agcaTGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTcc                                                             >  1:383030/1‑93 (MQ=255)

                             tcctGTTCCAGCTGCTCCTCCAGCCGTTGTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCg                                                            <  2:4619169/90‑1 (MQ=255)

                             gcaTGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCg                                                            >  1:4530281/1‑93 (MQ=255)

                               aTGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGc                                                          >  2:3382951/1‑93 (MQ=255)

                               aTGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGc                                                          >  1:737123/1‑93 (MQ=255)

                                  ttCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATa                                                       >  1:483357/1‑93 (MQ=255)

                                  ttCCACCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATa                                                       <  2:2010451/93‑1 (MQ=255)

                                   cccAACTCCTCCCCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAg                                                      <  1:1342048/92‑1 (MQ=255)

                                    ccAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAgg                                                     >  2:4092880/1‑93 (MQ=255)

                                     cAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAgct                                                    >  1:4500508/1‑91 (MQ=255)

                                     cAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGt                                                    >  2:4581488/1‑93 (MQ=255)

                                      aGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTc                                                   >  1:4038467/1‑93 (MQ=255)

                                      aGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTc                                                   <  1:777669/93‑1 (MQ=255)

                                       gctgctCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCt                                                  >  1:2551832/1‑93 (MQ=255)

                                         tgctCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGc                                                >  1:2851105/1‑93 (MQ=255)

                                           ctcctcCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGcac                                              >  2:3275484/1‑93 (MQ=255)

                                            tcctccAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGcaca                                             >  1:4720985/1‑93 (MQ=255)

                                             cctccAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCCATCGGTCCGGCGTAGGTCTGCACAt                                            >  2:3033264/1‑93 (MQ=255)

                                             cctccAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACAt                                            >  2:4407787/1‑93 (MQ=255)

                                             cctccAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACAt                                            >  2:357335/1‑93 (MQ=255)

                                              ctccAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATg                                           <  2:3100437/93‑1 (MQ=255)

                                              ctccAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATg                                           >  1:1054838/1‑93 (MQ=255)

                                               tccAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGa                                          >  1:1232803/1‑93 (MQ=255)

                                               tccAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGa                                          >  2:3704978/1‑93 (MQ=255)

                                               cccAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGa                                          <  2:29402/92‑1 (MQ=255)

                                                 cAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACg                                        <  2:4438304/93‑1 (MQ=255)

                                                  atccGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGa                                       <  2:718564/91‑1 (MQ=255)

                                                  aGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGa                                       >  2:4159448/1‑93 (MQ=255)

                                                  aGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGa                                       >  1:4373333/1‑93 (MQ=255)

                                                   gCCGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAg                                      >  1:1713864/1‑93 (MQ=255)

                                                    ccGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGt                                     >  1:349051/1‑93 (MQ=255)

                                                     cGTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTc                                    <  2:2616920/93‑1 (MQ=255)

                                                      gTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCg                                   >  2:2302374/1‑93 (MQ=255)

                                                      gTTTTTGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCg                                   >  2:3152810/1‑93 (MQ=255)

                                                       tttttGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGa                                  >  1:4181977/1‑93 (MQ=255)

                                                       tttttGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGa                                  <  1:1975498/93‑1 (MQ=255)

                                                        ttttGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGACGaa                                 >  1:1559992/1‑93 (MQ=255)

                                                          ttGCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAAcc                               <  2:2039354/93‑1 (MQ=255)

                                                            gCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCAt                             >  1:2983805/1‑93 (MQ=255)

                                                            gCTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCAt                             <  1:4581488/93‑1 (MQ=255)

                                                             cTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATc                            <  2:939682/93‑1 (MQ=255)

                                                             cTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATc                            <  2:4712233/93‑1 (MQ=255)

                                                             cTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATc                            <  1:2508017/93‑1 (MQ=255)

                                                             cTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATc                            <  1:2742244/93‑1 (MQ=255)

                                                             cTCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATc                            <  2:1288321/93‑1 (MQ=255)

                                                              tCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATcc                           >  1:233115/1‑93 (MQ=255)

                                                              tCAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATcc                           >  2:1224659/1‑93 (MQ=255)

                                                               cAGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCg                          >  1:1731689/1‑93 (MQ=255)

                                                                aGTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGc                         >  1:1829065/1‑93 (MQ=255)

                                                                 gTAATATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGcc                        >  1:2729487/1‑93 (MQ=255)

                                                                   aaTATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAg                      >  1:393793/1‑93 (MQ=255)

                                                                   aaTATCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAg                      >  1:2110152/1‑93 (MQ=255)

                                                                     tatGATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGtt                    <  1:2372713/93‑1 (MQ=255)

                                                                     tatCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGtt                    <  1:4082036/93‑1 (MQ=255)

                                                                     tatCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGtt                    >  1:672330/1‑93 (MQ=255)

                                                                     tatCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGtt                    <  2:3658393/93‑1 (MQ=255)

                                                                     tatCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCCGCACATGACGAGTCGAACCATCCGCCAGtt                    >  2:1122749/1‑93 (MQ=255)

                                                                      atCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGttt                   >  2:1243945/1‑93 (MQ=255)

                                                                      atCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGttt                   <  1:2053421/93‑1 (MQ=255)

                                                                       tCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTa                  <  2:2051131/93‑1 (MQ=255)

                                                                       tCATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTa                  <  1:3733882/93‑1 (MQ=255)

                                                                        cATGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTAt                 <  2:4086873/93‑1 (MQ=255)

                                                                         aTGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATg                <  2:3175143/93‑1 (MQ=255)

                                                                          tGCACAATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGa               <  2:4416368/93‑1 (MQ=255)

                                                                            cacaATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAAc             <  2:2999877/93‑1 (MQ=255)

                                                                             acaATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACa            <  2:1650559/93‑1 (MQ=255)

                                                                              caATACATAACACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACaa           <  1:2809600/93‑1 (MQ=255)

                                                                                aTACATAACACGAGCTTGTAGCCATAAATTTCAATCGGTCCGGCATAGGGCTGCTCATGACGAGTGGACCAATCCGCTCGTTTATGAACaaaa         >  2:3248166/1‑93 (MQ=255)

                                                                                      aaCATGAGCATGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCaaa   <  2:3920523/93‑1 (MQ=255)

                                                                                      aaCACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCaaa   <  2:4747780/93‑1 (MQ=255)

                                                                                      aaCACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCaaa   <  1:865390/93‑1 (MQ=255)

                                                                                      aaCACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCaaa   <  1:4680362/93‑1 (MQ=255)

                                                                                      aaCACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCaaa   >  2:2411967/1‑93 (MQ=255)

                                                                                      aaCACGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCaaa   <  1:1951389/93‑1 (MQ=255)

                                                                                       acacGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAg  <  2:4545432/93‑1 (MQ=255)

                                                                                       acacGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAg  >  1:1261799/1‑93 (MQ=255)

                                                                                       acacGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAg  >  1:1973182/1‑93 (MQ=255)

GTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAG  >  NC_000913/1405036‑1405215

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