BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,608,927	T→C	D147D (GAT→GAC)	hyfG →	hydrogenase 4, subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,608,927	0	T	C	98.9%	622.7 / NA	176	D147D (GAT→GAC)	hyfG	hydrogenase 4, subunit
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (77/97); total (78/98)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 177 total aligned reads displayed.

TCTGGGGCGAGCGAGAAATTCGCGATATGTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGC  >  NC_000913/2608830‑2609023

tCTGGGGCGAGCGAGAAATTCGCGATATGTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGacat                                                                                                <  2:1796628/100‑1 (MQ=255)

 cTGGGGCGAGCGAGAAATTCGCGATATGTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATg                                                                                               >  2:2391470/1‑100 (MQ=255)

  tGGGGCGAGCGAGAAATTCGCGATATGTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGTCTGGCCGGAAGACATGc                                                                                              >  2:505122/1‑100 (MQ=255)

   ggggCGAGCGAGAAATTCGCGATATGTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCa                                                                                             <  2:2525764/100‑1 (MQ=255)

     ggCGAGCGAGAAATTCGCGATATGTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATc                                                                                           <  1:4764714/100‑1 (MQ=255)

     ggCGAGCGAGAAATTCGCGATATGTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATc                                                                                           >  1:354762/1‑100 (MQ=255)

      gcgagcgaGAAAGTCGCGATAGGTACGGGCTTATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACCGGCCGGAAGACATGCATcc                                                                                          <  2:1400138/100‑1 (MQ=255)

        tctcaagaAATTCGCGATATGTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGc                                                                                        <  1:2205080/95‑1 (MQ=255)

        gagcgagAAATTCGCGATATGTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGc                                                                                        <  1:4837956/100‑1 (MQ=255)

        gagcgagAAATTCGCGATATGTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGc                                                                                        <  2:3481619/100‑1 (MQ=255)

        gagcgagAAATTCGCGATATGTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGc                                                                                        <  2:3614258/100‑1 (MQ=255)

         agcgagAAATTCGCGATATGTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCt                                                                                       <  1:2951078/100‑1 (MQ=255)

          gcgagAAATTCGCGATATGTACGCGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGTAGACATGCATCCGCTg                                                                                      >  2:4251109/1‑100 (MQ=255)

           cgagAAATTCGCGATATGTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTgc                                                                                     >  2:4519797/1‑100 (MQ=255)

           cgagAAATTCGCGATATGTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCAgc                                                                                     >  2:2178552/1‑100 (MQ=255)

             agaAATTCGCGATATGTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTgcgc                                                                                   >  2:791243/1‑100 (MQ=255)

                 aTTCGCGATATGTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAg                                                                               <  1:3939102/100‑1 (MQ=255)

                   tCGCGATATGTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGAt                                                                             <  2:3464205/100‑1 (MQ=255)

                    cgcgATATGTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATg                                                                            <  1:3185100/100‑1 (MQ=255)

                        atatGTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGAt                                                                        <  1:1390094/100‑1 (MQ=255)

                        atatGTACGGGCTGATTCCGGCTGGCCTGCCGGATCCGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGAt                                                                        <  2:4806048/100‑1 (MQ=255)

                         tatGTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATg                                                                       <  2:2111692/100‑1 (MQ=255)

                            gTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGAt                                                                    >  1:4473305/1‑100 (MQ=255)

                                gggCTGATTGCGGTTGGCCTGCCAGATCAGCGTCGCCTGGTGTTGCCCGATGACTGACCGGTAGACATGCATGCGCTGCGCACAGATGCGATGGATTATc                                                                >  2:3133582/1‑100 (MQ=255)

                                 ggCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCg                                                               <  2:3645184/100‑1 (MQ=255)

                                  gCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGa                                                              >  2:1289489/1‑100 (MQ=255)

                                    tGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACt                                                            >  2:2811881/1‑100 (MQ=255)

                                       ttCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTgcg                                                         >  1:3228417/1‑100 (MQ=255)

                                       ttCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTgcg                                                         >  1:1244386/1‑100 (MQ=255)

                                       ttCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTgcg                                                         >  2:160186/1‑100 (MQ=255)

                                        tCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTgcgc                                                        >  1:980177/1‑100 (MQ=255)

                                          cGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGccc                                                      <  2:2130742/100‑1 (MQ=255)

                                           ggTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCt                                                     >  1:3710274/1‑100 (MQ=255)

                                              tGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGaa                                                  <  2:2991226/100‑1 (MQ=255)

                                              tGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGaa                                                  >  2:4537082/1‑100 (MQ=255)

                                               ggCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAAc                                                 <  2:3257037/100‑1 (MQ=255)

                                               ggCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAAc                                                 <  1:2041291/100‑1 (MQ=255)

                                               ggCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAAc                                                 <  2:5024542/100‑1 (MQ=255)

                                               ggCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAAc                                                 <  1:1298318/100‑1 (MQ=255)

                                                gCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAAcc                                                <  2:2584703/100‑1 (MQ=255)

                                                gCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAAcc                                                <  1:3138425/100‑1 (MQ=255)

                                                gCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAAcc                                                <  2:876678/100‑1 (MQ=255)

                                                  cTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACcga                                              >  1:1548813/1‑100 (MQ=255)

                                                   tGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACcgac                                             <  2:1880251/100‑1 (MQ=255)

                                                   tGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACcgac                                             <  1:1204529/100‑1 (MQ=255)

                                                    gCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACcgacg                                            >  1:1075606/1‑100 (MQ=255)

                                                      cGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCgacgac                                          <  1:1865763/100‑1 (MQ=255)

                                                      cGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCgacgac                                          <  2:2259384/100‑1 (MQ=255)

                                                      cGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACCGCGCCCTGAACCgacgac                                          <  2:1939512/100‑1 (MQ=255)

                                                       ggATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTCACGACTGCGCCCTGAACCGACGACt                                         >  1:4180132/1‑100 (MQ=255)

                                                       ggATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACt                                         >  2:4008337/1‑100 (MQ=255)

                                                        gATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTg                                        <  2:4415399/100‑1 (MQ=255)

                                                         aTCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGa                                       <  1:4053641/100‑1 (MQ=255)

                                                         aTCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGa                                       <  2:2089841/100‑1 (MQ=255)

                                                          tCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGAt                                      <  1:803877/100‑1 (MQ=255)

                                                            aGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTc                                    <  2:4399661/100‑1 (MQ=255)

                                                             gCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTcc                                   <  1:1708772/100‑1 (MQ=255)

                                                              cGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCg                                  <  1:700927/100‑1 (MQ=255)

                                                               gTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGa                                 <  2:1162590/100‑1 (MQ=255)

                                                               gTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGa                                 >  1:3884908/1‑100 (MQ=255)

                                                                tCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGaa                                <  1:4301014/100‑1 (MQ=255)

                                                                tCGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGaa                                <  1:3435753/100‑1 (MQ=255)

                                                                 cGCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGaaa                               <  1:2486858/100‑1 (MQ=255)

                                                                  gCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAAc                              >  2:4955360/1‑100 (MQ=255)

                                                                  gCCTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAAc                              <  1:3362242/100‑1 (MQ=255)

                                                                   ccTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCGAAGATGCGATGGATTATCCACTGCGCCCTGAACCGACGACTGATTCCGAAACg                             <  2:4690399/100‑1 (MQ=255)

                                                                   ccTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACg                             >  2:3892187/1‑100 (MQ=255)

                                                                   ccTGGTGTTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACg                             >  1:1742620/1‑100 (MQ=255)

                                                                         gtTGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCg                       >  2:1652317/1‑100 (MQ=255)

                                                                          ttGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGt                      <  1:4333289/100‑1 (MQ=255)

                                                                           tGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGtt                     >  1:1894538/1‑100 (MQ=255)

                                                                           tGCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGtt                     >  2:2645779/1‑100 (MQ=255)

                                                                            gCCCGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGttt                    >  1:1023005/1‑100 (MQ=255)

                                                                             cccGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTa                   >  1:871513/1‑100 (MQ=255)

                                                                             cccGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTa                   <  1:3117475/100‑1 (MQ=255)

                                                                             cccGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTa                   <  2:1666050/100‑1 (MQ=255)

                                                                              ccGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTCCGCCCTGAACCGACGACTGGTTCCGAAACGTATCCGTTTAt                  <  2:4834397/100‑1 (MQ=255)

                                                                               cGATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATc                 >  1:4268671/1‑100 (MQ=255)

                                                                                gATGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCa                >  2:4073031/1‑100 (MQ=255)

                                                                                 aTGACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCaa               >  1:2029351/1‑100 (MQ=255)

                                                                                   gACTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATg             <  1:4396396/100‑1 (MQ=255)

                                                                                    aCTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGa            >  1:2845692/1‑100 (MQ=255)

                                                                                    aCTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGa            >  1:2846375/1‑100 (MQ=255)

                                                                                     cTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATGCGTTTATCAATGAg           >  2:4637771/1‑100 (MQ=255)

                                                                                     cTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAg           <  2:951922/100‑1 (MQ=255)

                                                                                     cTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAg           <  2:3675751/100‑1 (MQ=255)

                                                                                     cTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAg           <  1:1235476/100‑1 (MQ=255)

                                                                                     cTGGCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCCCCCTGAATCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAg           <  1:4583456/100‑1 (MQ=255)

                                                                                     cTGGCCGGAAGACATGCAGCCGCTGCGGAAAGTGTCGATGGATTATCGGCGGAGCCCTTAACCGACGACTGATTCCGAAACGTATCCGTTTATCagcgag           <  2:4305979/100‑6 (MQ=255)

                                                                                       ggCCGGGAGACATGCATACGCTGCGCCAAAATGCGGAGGGATATCGACTGCGCCCTGAACCGACGACCGATTTCCAAACGGATCCCCTTATCCATGAggg         >  2:4511663/1‑100 (MQ=255)

                                                                                       ggCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGAATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAggg         <  1:3269828/100‑1 (MQ=255)

                                                                                        gCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGc        <  2:4498878/100‑1 (MQ=255)

                                                                                        gCCGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGc        <  2:4557588/100‑1 (MQ=255)

                                                                                        gCCGGAAGACATGCATCCGCTGCGCAAACATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGc        <  1:4684279/100‑1 (MQ=255)

                                                                                         ccGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCa       <  2:4605104/100‑1 (MQ=255)

                                                                                          cGGAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCaa      >  1:4008488/1‑100 (MQ=255)

                                                                                            gAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACa    >  1:2925896/1‑100 (MQ=255)

                                                                                            gAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACa    <  1:1681231/100‑1 (MQ=255)

                                                                                            gAAGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACa    >  2:3766974/1‑100 (MQ=255)

                                                                                              aGACATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGc  <  1:4411905/100‑1 (MQ=255)

TCTGGGGCGAGCGAGAAATTCGCGATATGTACGGGCTGATTCCGGTTGGCCTGCCGGATCAGCGTCGCCTGGTGTTGCCCGATGACTGGCCGGAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGC  >  NC_000913/2608830‑2609023

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