BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,739,740	A→G	M141T (ATG→ACG)	yiaJ ←	putative regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,739,740	0	A	G	99.6%	687.9 / ‑4.8	220	M141T (ATG→ACG)	yiaJ	putative regulator
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (109/110); total (109/111)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

Only 100 of 235 total aligned reads displayed.

TAGTCCGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCATCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATGTTC  >  NC_000913/3739650‑3739829

tAGTCCGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTc                                                                                         >  1:4473072/1‑93 (MQ=255)

tAGTCCGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTc                                                                                         <  1:1927343/93‑1 (MQ=255)

 aGTCCGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTcc                                                                                        <  2:720291/93‑1 (MQ=255)

   tCCGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCg                                                                                      <  1:1108341/93‑1 (MQ=255)

    ccGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCgg                                                                                     <  1:2805842/93‑1 (MQ=255)

     cGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGt                                                                                    <  1:2922966/93‑1 (MQ=255)

       ggTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTATAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtg                                                                                  <  1:1944183/93‑1 (MQ=255)

       ggTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGGtg                                                                                  >  2:1828349/1‑93 (MQ=255)

        gTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtgt                                                                                 <  2:1866316/93‑1 (MQ=255)

         tGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtgtg                                                                                <  1:1828145/93‑1 (MQ=255)

         tGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtgtg                                                                                <  1:4182745/93‑1 (MQ=255)

         tGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtgtg                                                                                >  1:4427931/1‑93 (MQ=255)

           aCCAAACGCCATGTAGATCTTGCTCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTggg                                                                              >  1:967165/1‑93 (MQ=255)

           aCCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTggg                                                                              >  2:2745987/1‑93 (MQ=255)

           aCCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTggg                                                                              <  1:4332516/93‑1 (MQ=255)

            ccAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGt                                                                             <  1:782680/93‑1 (MQ=255)

             cAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGtt                                                                            >  1:517513/1‑93 (MQ=255)

              aaaCGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTc                                                                           <  1:3618438/93‑1 (MQ=255)

              aaaCGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTc                                                                           >  2:2661843/1‑93 (MQ=255)

                  gCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGc                                                                       >  2:4571210/1‑93 (MQ=255)

                     aTGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTta                                                                    >  2:809849/1‑93 (MQ=255)

                       gcaGATCTTGCCCATAGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTtata                                                                  <  2:47221/91‑1 (MQ=255)

                       gTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTtata                                                                  <  2:1762505/93‑1 (MQ=255)

                       gTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTtata                                                                  >  1:6200/1‑93 (MQ=255)

                         aGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATaaa                                                                <  1:2964589/93‑1 (MQ=255)

                          gATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAAt                                                               >  2:4498047/1‑93 (MQ=255)

                            tCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCa                                                             >  2:1161814/1‑93 (MQ=255)

                            tCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCa                                                             >  2:1991467/1‑93 (MQ=255)

                            tCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCa                                                             >  2:2344422/1‑93 (MQ=255)

                            tCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCa                                                             >  1:2683706/1‑93 (MQ=255)

                             cTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCaa                                                            >  1:3428897/1‑93 (MQ=255)

                             cTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCaa                                                            >  2:3499776/1‑93 (MQ=255)

                              ttGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTGCAGCTTATAAATCaaa                                                           >  2:4161348/1‑93 (MQ=255)

                              ttGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATaaaa                                                           >  2:3954482/1‑93 (MQ=255)

                              ttGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCaaa                                                           >  2:1832996/1‑93 (MQ=255)

                               tGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCaaaa                                                          <  1:1737010/93‑1 (MQ=255)

                                gCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAAt                                                         >  2:3220713/1‑93 (MQ=255)

                                 cccATTGCGGATCAGTAGAGCGGCATACGCTGGCCAATATATGCACGGGTTCGCAGCCTCCCGGTTGTGGGTTCCATCTTATAAATCAAAATa                                                        >  2:4339899/1‑93 (MQ=255)

                                 cccATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATa                                                        <  2:1455308/93‑1 (MQ=255)

                                  ccATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAg                                                       <  1:3215437/93‑1 (MQ=255)

                                  ccATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAg                                                       >  2:3423688/1‑93 (MQ=255)

                                  ccATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAg                                                       >  1:2837313/1‑93 (MQ=255)

                                  ccATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAg                                                       >  2:2371235/1‑93 (MQ=255)

                                   cATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGc                                                      <  1:2744134/93‑1 (MQ=255)

                                   cATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGc                                                      >  2:1792714/1‑93 (MQ=255)

                                    aTTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCg                                                     >  1:929876/1‑93 (MQ=255)

                                    aTTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCg                                                     >  1:1463069/1‑93 (MQ=255)

                                     ttGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGt                                                    >  2:2834358/1‑93 (MQ=255)

                                     ttGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGt                                                    <  2:4174096/93‑1 (MQ=255)

                                         ggAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATc                                                <  2:3251838/93‑1 (MQ=255)

                                         ggAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATc                                                >  2:210552/1‑93 (MQ=255)

                                          gAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCg                                               >  1:3621953/1‑93 (MQ=255)

                                          gAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCg                                               >  1:2340801/1‑93 (MQ=255)

                                             cAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCt                                            <  2:2119701/93‑1 (MQ=255)

                                             cAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCt                                            >  1:2091015/1‑93 (MQ=255)

                                              aGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCtt                                           >  1:384010/1‑93 (MQ=255)

                                              aGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCtt                                           <  1:840653/93‑1 (MQ=255)

                                               gTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTc                                          >  1:1709667/1‑93 (MQ=255)

                                               gTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTc                                          >  2:3513025/1‑93 (MQ=255)

                                                  gagCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTcgcg                                       >  2:938939/1‑93 (MQ=255)

                                                  gagCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTcgcg                                       >  1:4220424/1‑93 (MQ=255)

                                                      ggCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTg                                   >  2:1576288/1‑93 (MQ=255)

                                                          tatGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGgaga                               >  2:3882437/1‑93 (MQ=255)

                                                          tatGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGgaga                               >  2:2408253/1‑93 (MQ=255)

                                                           atGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGaa                              <  1:1040178/93‑1 (MQ=255)

                                                           atGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGaa                              >  2:230508/1‑93 (MQ=255)

                                                           atGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGaa                              <  1:2047555/93‑1 (MQ=255)

                                                            tGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAg                             >  2:1150691/1‑93 (MQ=255)

                                                               tGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCCGGGGAAGTTa                          >  2:2898524/1‑93 (MQ=255)

                                                                ggCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGGGGAGTTaa                         >  2:299903/1‑93 (MQ=255)

                                                                 gCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAAt                        <  1:2201655/93‑1 (MQ=255)

                                                                 gCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAAt                        <  2:2986049/93‑1 (MQ=255)

                                                                 gCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAAt                        >  2:2856292/1‑93 (MQ=255)

                                                                   cAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATgg                      >  2:2460063/1‑93 (MQ=255)

                                                                   cAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATgg                      >  1:2537972/1‑93 (MQ=255)

                                                                   cAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATgg                      <  1:1395232/93‑1 (MQ=255)

                                                                   cAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATgg                      >  1:1177125/1‑93 (MQ=255)

                                                                    aaTATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGt                     <  1:4625078/93‑1 (MQ=255)

                                                                     atatAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGtt                    <  1:2408869/93‑1 (MQ=255)

                                                                     atatAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGtt                    >  1:1583307/1‑93 (MQ=255)

                                                                       ataGGCACGGGTTCGCAGCGGCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTc                  <  1:2986203/93‑1 (MQ=255)

                                                                        taGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTAGAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCa                 >  1:3299541/1‑93 (MQ=255)

                                                                          ggCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCAcc               <  1:3560588/93‑1 (MQ=255)

                                                                           gCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCa              <  2:160099/93‑1 (MQ=255)

                                                                           gCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCa              <  2:2750738/93‑1 (MQ=255)

                                                                           gCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCa              <  1:1971663/93‑1 (MQ=255)

                                                                            cACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAg             <  2:1015574/93‑1 (MQ=255)

                                                                             aCGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGt            <  2:2549976/93‑1 (MQ=255)

                                                                             aCGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGt            >  1:2298129/1‑93 (MQ=255)

                                                                                gtttCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTGCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGc         <  1:4315535/91‑1 (MQ=255)

                                                                                ggTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGc         <  2:4034052/93‑1 (MQ=255)

                                                                                ggTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGc         <  2:516779/93‑1 (MQ=255)

                                                                                 gTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCg        <  2:1953185/93‑1 (MQ=255)

                                                                                  ttCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGa       <  2:2196842/93‑1 (MQ=255)

                                                                                  ttCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGGTTCACCAGTGGCGa       >  2:1786925/1‑93 (MQ=255)

                                                                                   tCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGAt      >  1:4070940/1‑93 (MQ=255)

                                                                                    cGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATg     <  2:1359169/93‑1 (MQ=255)

                                                                                      cagcgTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATGtt   >  2:4442922/1‑93 (MQ=255)

                                                                                      cagcgTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATGtt   >  2:3035105/1‑93 (MQ=255)

                                                                                       agcgTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATGTTc  >  1:3791122/1‑93 (MQ=255)

TAGTCCGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCATCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATGTTC  >  NC_000913/3739650‑3739829

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