BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,379,796	C→T	76.1%	R283H (CGC→CAC)	menF ←	isochorismate synthase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,379,796	0	C	T	76.1%	266.6 / 94.4	160	R283H (CGC→CAC)	menF	isochorismate synthase 2
Reads supporting (aligned to +/- strand): ref base C (19/19); new base T (64/57); total (84/76)
Fisher's exact test for biased strand distribution p-value = 8.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 164 total aligned reads displayed.

CTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCTTCAGTACGC  >  NC_000913/2379698‑2379894

cTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTcgcg                                                                                                   >  2:1924985/1‑100 (MQ=255)

    acgcagacgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTgg                                                                                               >  1:3788997/1‑100 (MQ=255)

     cgcagacgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGATAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGt                                                                                              >  2:1176852/1‑100 (MQ=255)

     cgcagacgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGt                                                                                              >  1:1031638/1‑100 (MQ=255)

     cgcagacgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGt                                                                                              >  2:2630046/1‑100 (MQ=255)

       cagacgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCTACCACCAGCATGTTCTCGCGCTGGttt                                                                                            >  2:752650/1‑100 (MQ=255)

       cagacgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGttt                                                                                            >  1:161911/1‑100 (MQ=255)

       cagacgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGttt                                                                                            >  2:4095647/1‑100 (MQ=255)

        agacgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGtttt                                                                                           <  2:373105/100‑1 (MQ=255)

         gacgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGttttt                                                                                          >  1:2402853/1‑100 (MQ=255)

           cgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTat                                                                                        <  2:368633/100‑1 (MQ=255)

           cgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTat                                                                                        <  2:1272986/100‑1 (MQ=255)

           cgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTat                                                                                        >  1:3441772/1‑100 (MQ=255)

           cgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTat                                                                                        <  1:3243298/100‑1 (MQ=255)

           cgcagTACCTGCGGCGGTAAAACATCCAGCGTCTCGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTat                                                                                        <  2:2197741/100‑1 (MQ=255)

            gcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTatc                                                                                       <  2:1104961/100‑1 (MQ=255)

             cagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATGTCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTatca                                                                                      >  1:1743187/1‑100 (MQ=255)

               gTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTatcatc                                                                                    >  2:588672/1‑100 (MQ=255)

                 aCCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCg                                                                                  <  1:147901/100‑1 (MQ=255)

                  ccTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGc                                                                                 >  1:2131101/1‑100 (MQ=255)

                    tGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCa                                                                               <  1:105060/100‑1 (MQ=255)

                    tGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCa                                                                               >  2:1512707/1‑100 (MQ=255)

                       ggcggTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCa                                                                            <  1:3060052/100‑1 (MQ=255)

                       ggcggTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCa                                                                            <  1:1127798/100‑1 (MQ=255)

                        gcggTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAg                                                                           >  2:427058/1‑100 (MQ=255)

                        gcggTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAg                                                                           <  1:3451218/100‑1 (MQ=255)

                         cggTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGc                                                                          <  2:664963/100‑1 (MQ=255)

                          ggTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGGGCGGGTTTTTATCATCCGCCATCAGcc                                                                         >  1:189322/1‑100 (MQ=255)

                          ggTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGcc                                                                         <  1:1152304/100‑1 (MQ=255)

                          ggTAAAACATCCAGCGGCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGcc                                                                         <  1:5489/100‑1 (MQ=255)

                          ggTAAAACATACAGCGTCTGTGTATTCGCTTGTAATCGTTGACAGATATCTTCTACTTCTATCATGTTCTCTTTGTTGTTTTTATCATCCGCCATAAg                                                                           <  2:668579/98‑1 (MQ=255)

                          cgTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGcc                                                                         <  1:1466002/99‑1 (MQ=255)

                              aaaCATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCActc                                                                     <  1:321904/100‑1 (MQ=255)

                               aaCATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCActct                                                                    <  1:1922762/100‑1 (MQ=255)

                               aaCATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCActct                                                                    <  2:3896896/100‑1 (MQ=255)

                                aCATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACtctc                                                                   >  1:3154792/1‑100 (MQ=255)

                                 cATCCAGCGTCTTCGTATCGGCCGGTAATCGTTGCCACATATTTTCCACCATCAGCATGTTCTCGTGCTGGTTTGTATCATCCGCCATCAGCCACTCTcc                                                                  <  2:2342917/100‑1 (MQ=255)

                                 cATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTcc                                                                  <  2:85304/100‑1 (MQ=255)

                                  aTCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCt                                                                 >  1:409016/1‑100 (MQ=255)

                                  aTCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCt                                                                 <  2:2611966/100‑1 (MQ=255)

                                   tCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTa                                                                >  1:1307126/1‑100 (MQ=255)

                                    ccAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTaa                                                               >  2:3022327/1‑100 (MQ=255)

                                    ccAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTaa                                                               <  2:560642/100‑1 (MQ=255)

                                     cAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAAc                                                              >  1:3171570/1‑100 (MQ=255)

                                        cGTCTGGGTACCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTTTCATCCGCCATCAGCCACTCTCCTAActgc                                                           <  2:2735862/100‑1 (MQ=255)

                                             gggTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTgcgc                                                      >  2:1969817/1‑100 (MQ=255)

                                               gTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCt                                                    >  2:3292253/1‑100 (MQ=255)

                                                tATCGGCCTGTAATCGTTGACAGATATCTTCCACCCCCAGCATGTTCTAGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTg                                                   >  1:2426229/1‑100 (MQ=255)

                                                  tCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCt                                                 <  1:799277/100‑1 (MQ=255)

                                                  tCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCt                                                 <  1:990933/100‑1 (MQ=255)

                                                    ggCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTa                                               >  1:2577736/1‑100 (MQ=255)

                                                        tGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTatcat                                           <  1:999982/100‑1 (MQ=255)

                                                        tGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTatca                                            <  2:2748025/99‑1 (MQ=255)

                                                         gTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTatcatc                                          >  1:154479/1‑100 (MQ=255)

                                                         gTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTatcatc                                          >  1:2303055/1‑100 (MQ=255)

                                                           aaTCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAg                                        <  2:1138170/100‑1 (MQ=255)

                                                            aTCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAgg                                       >  2:1206488/1‑100 (MQ=255)

                                                            aTCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAgg                                       <  2:3788968/100‑1 (MQ=255)

                                                             tCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGa                                      <  2:2115944/100‑1 (MQ=255)

                                                               gTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGatt                                    >  2:1584561/1‑100 (MQ=255)

                                                               gTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGatt                                    >  1:1186754/1‑100 (MQ=255)

                                                                 tGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGattat                                  >  2:655867/1‑100 (MQ=255)

                                                                    cAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGattattgg                               >  1:2146330/1‑98 (MQ=255)

                                                                    cAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTg                               >  2:2774716/1‑100 (MQ=255)

                                                                    cAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTg                               >  2:4021469/1‑100 (MQ=255)

                                                                     aGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGc                              <  1:3704049/100‑1 (MQ=255)

                                                                      gATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGc                              <  2:3474160/99‑1 (MQ=255)

                                                                      gATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCt                             <  1:3769802/100‑1 (MQ=255)

                                                                      gATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCt                             >  1:2209995/1‑100 (MQ=255)

                                                                       atatCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTa                            >  1:3932158/1‑100 (MQ=255)

                                                                       atatCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTa                            >  1:1728340/1‑100 (MQ=255)

                                                                       atatCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTa                            >  2:601297/1‑100 (MQ=255)

                                                                       atatCTTCCACCACCAGCATGCTCTCGCGCGGGTTTCTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCt                             <  2:2818345/99‑1 (MQ=255)

                                                                          tCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTg                         >  1:2466738/1‑100 (MQ=255)

                                                                          tCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTg                         >  1:365763/1‑100 (MQ=255)

                                                                           cTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGt                        >  2:1636316/1‑100 (MQ=255)

                                                                            ttCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGtt                       >  1:1291216/1‑100 (MQ=255)

                                                                             tCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTc                      >  1:402827/1‑100 (MQ=255)

                                                                               caccaccaGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTccc                    >  1:801841/1‑100 (MQ=255)

                                                                               caccaccaGCATGTTCTCGTGCCGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTccc                    >  2:2411561/1‑100 (MQ=255)

                                                                                accaccaGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAGCTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCg                   >  2:804445/1‑100 (MQ=255)

                                                                                 ccaccaGCATGTTCTCGCGCTGGTTTTTATCATCCGCCCTCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGc                  >  1:89890/1‑100 (MQ=255)

                                                                                 ccaccaGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGc                  >  2:974302/1‑100 (MQ=255)

                                                                                  caccaGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCg                   >  2:4143290/1‑98 (MQ=255)

                                                                                  caccaGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGcc                 <  2:2348904/100‑1 (MQ=255)

                                                                                   accaGCATGTTCTCGTGCTAGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCa                >  2:1653525/1‑100 (MQ=255)

                                                                                   accaGCATGTTCTCGCGCTGGTTTTTCTCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTCCTTCTCCTCCGGCTTATTTACTACTGTTcccccc                 >  2:3684305/1‑99 (MQ=255)

                                                                                    ccaGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAg               <  2:1084394/100‑1 (MQ=255)

                                                                                    ccaGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCa                <  1:3532540/99‑1 (MQ=255)

                                                                                      aGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAgcg             <  2:2134911/100‑1 (MQ=255)

                                                                                       gCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCCCTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAgcgc            <  1:3269977/100‑1 (MQ=255)

                                                                                         aTGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCtt          >  2:1294079/1‑100 (MQ=255)

                                                                                          tGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCTTc         >  2:3205918/1‑100 (MQ=255)

                                                                                          tGTTCTCGTGCTAGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCTTc         <  1:1653534/100‑1 (MQ=255)

                                                                                            ttCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCTTCAg       <  1:2749286/100‑1 (MQ=255)

                                                                                             tctcGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCTTCAGt      >  1:111783/1‑100 (MQ=255)

                                                                                             tctcGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCGTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCTTCAGt      >  2:2667677/1‑100 (MQ=255)

                                                                                              ctcGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCTTCAGTa     <  2:3924419/100‑1 (MQ=255)

                                                                                                cgcgCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCTTCAGTACg   >  1:2869777/1‑100 (MQ=255)

                                                                                                 gcgcTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCTTCAGTACGc  <  2:1754142/100‑1 (MQ=255)

CTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCTTCAGTACGC  >  NC_000913/2379698‑2379894

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