BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,577,196	T→G	M133L (ATG→CTG)	gntK ←	gluconate kinase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,577,196	0	T	G	99.3%	512.5 / ‑4.1	143	M133L (ATG→CTG)	gntK	gluconate kinase 2
Reads supporting (aligned to +/- strand): ref base T (1/0); new base G (68/74); total (69/74)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

Only 100 of 145 total aligned reads displayed.

TGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACATTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTAC  >  NC_000913/3577098‑3577291

tGCCACAACACCTTCCAGCTGTTGATCGATATCTACCACCAGTACATCTGTTTCGTCCGGACTCGGCTCTTGTAGCTTTTATAACTGCGTCACCATCAtt                                                                                                >  2:2326797/1‑100 (MQ=255)

  ccACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTg                                                                                              <  1:528078/100‑1 (MQ=255)

  ccACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTg                                                                                              >  1:2658547/1‑100 (MQ=255)

    acaacaCCTTCGAGCGGTTGATCGATATCCACCACCAGTAAATCGGTTTCGTCCGCACCCGGCTCCTGCAGCATTTCAAACTGCGTCACCAACAGTTggg                                                                                            <  1:3595456/100‑1 (MQ=255)

          ccTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCGCCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTaa                                                                                      <  2:3841557/100‑1 (MQ=255)

          ccTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTaa                                                                                      <  2:2791819/100‑1 (MQ=255)

          ccTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTaa                                                                                      <  1:1489688/100‑1 (MQ=255)

            ttCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAg                                                                                    <  1:1899881/100‑1 (MQ=255)

             tCCAGCGGTTGCTCGATCTCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGa                                                                                   <  2:181029/100‑1 (MQ=255)

             tCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGa                                                                                   <  2:2808092/100‑1 (MQ=255)

             tCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGa                                                                                   <  1:1783518/100‑1 (MQ=255)

             tCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGa                                                                                   >  2:2733365/1‑100 (MQ=255)

             tCCAGCGGTTGATCGATATCCACCACCAGTACAGCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGa                                                                                   <  1:4037808/100‑1 (MQ=255)

               cAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGaaa                                                                                 >  1:1318814/1‑100 (MQ=255)

                aGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAAt                                                                                >  2:2428843/1‑100 (MQ=255)

                aGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAAt                                                                                >  2:1658626/1‑100 (MQ=255)

                   ggTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCCACAGTTGGGTTTTAAAGAAATGGc                                                                             >  2:3337190/1‑100 (MQ=255)

                   ggTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGc                                                                             >  1:1363265/1‑100 (MQ=255)

                   ggTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGc                                                                             >  2:670348/1‑100 (MQ=255)

                    gTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAGACTGCGTCACCATCAGTTGGGTTTTAAAGAAATGGcc                                                                            <  1:1060623/100‑1 (MQ=255)

                      tGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCtt                                                                          >  1:2236657/1‑100 (MQ=255)

                       gATCGATATCCACCACCAGTACATTGGTTTCGTCCGCACCCGGCTCTTGCAGGGTTTCAAACTGCGTCACAAACAGTTGGGTTTTAAAGCAATGGTCttt                                                                         >  2:360984/1‑100 (MQ=255)

                       gATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCttt                                                                         >  2:2056306/1‑100 (MQ=255)

                       gATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCttt                                                                         >  2:3431032/1‑100 (MQ=255)

                        aTCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTg                                                                        >  1:3373808/1‑100 (MQ=255)

                          cGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTgcg                                                                      <  1:2052249/100‑1 (MQ=255)

                          cGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTgcg                                                                      <  1:440724/100‑1 (MQ=255)

                          cGATATCCACCACCAGTACATCAGTTTCCTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTgcg                                                                      <  1:3710805/100‑1 (MQ=255)

                           gATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTgcgc                                                                     >  2:2527911/1‑100 (MQ=255)

                                 caccaccaGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTc                                                               >  1:1517015/1‑100 (MQ=255)

                                 cacaaccaGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTc                                                               >  2:3676461/1‑100 (MQ=255)

                                    caccaGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAgg                                                            >  1:2142669/1‑100 (MQ=255)

                                     accaGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAggc                                                           >  2:1188231/1‑100 (MQ=255)

                                      ccaGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAggcg                                                          >  2:3112996/1‑100 (MQ=255)

                                      ccaGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAggcg                                                          >  1:2163365/1‑100 (MQ=255)

                                      ccaGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAggcg                                                          >  1:2200335/1‑100 (MQ=255)

                                       caGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAggcgg                                                         <  2:2513927/100‑1 (MQ=255)

                                       caGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAggcgg                                                         <  1:562723/100‑1 (MQ=255)

                                        aGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAggcggc                                                        <  1:1189583/100‑1 (MQ=255)

                                         gTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCt                                                       >  2:544119/1‑100 (MQ=255)

                                         gTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCt                                                       <  2:2964144/100‑1 (MQ=255)

                                            cATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTc                                                    >  2:2061984/1‑100 (MQ=255)

                                            cATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTc                                                    <  1:3717027/100‑1 (MQ=255)

                                             aTCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCa                                                   >  2:1046785/1‑100 (MQ=255)

                                              tCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCaa                                                  >  1:3194715/1‑100 (MQ=255)

                                               cGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCTCTTTCAGGCGGCTTTCAAt                                                 <  2:3612013/100‑1 (MQ=255)

                                                   ttCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATcgca                                             <  1:3891749/100‑1 (MQ=255)

                                                   ttCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATcaca                                             <  1:2360420/100‑1 (MQ=255)

                                                    tCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACAt                                            <  2:3853731/100‑1 (MQ=255)

                                                    tCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACAt                                            <  1:1161233/100‑1 (MQ=255)

                                                      gTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCa                                          <  2:391328/100‑1 (MQ=255)

                                                      gTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCa                                          >  2:3185994/1‑100 (MQ=255)

                                                      gTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCa                                          >  1:3913571/1‑100 (MQ=255)

                                                      gTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCa                                          <  2:1490898/100‑1 (MQ=255)

                                                        ccGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCaaa                                        >  1:1953111/1‑100 (MQ=255)

                                                         cGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCaaaa                                       <  1:815276/100‑1 (MQ=255)

                                                           cACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATc                                     <  1:486621/100‑1 (MQ=255)

                                                             cccGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGc                                   >  2:1647000/1‑100 (MQ=255)

                                                              ccGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGcc                                  <  1:1632937/100‑1 (MQ=255)

                                                               cGGCTCCTACAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTGAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAACATCGCCt                                 <  2:3756892/100‑1 (MQ=255)

                                                                ggCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCtt                                >  1:3540544/1‑100 (MQ=255)

                                                                 gCTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCttt                               <  2:3714884/100‑1 (MQ=255)

                                                                  cTCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTc                              <  1:2357399/100‑1 (MQ=255)

                                                                   tCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCa                             <  1:1091689/100‑1 (MQ=255)

                                                                   tCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCa                             <  2:1529603/100‑1 (MQ=255)

                                                                   tCCTGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCa                             <  1:2837421/100‑1 (MQ=255)

                                                                      tGCAGCGTTTCAAACTGCGTCACCAACAGTTGTGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAAt                          <  2:1699923/100‑1 (MQ=255)

                                                                      tGCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAAt                          <  1:549207/100‑1 (MQ=255)

                                                                       gCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATa                         <  1:1686555/100‑1 (MQ=255)

                                                                       gCAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATa                         <  2:1126136/100‑1 (MQ=255)

                                                                       gCAGCGTTTCAAACTGCGGCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATa                         <  2:1163667/100‑1 (MQ=255)

                                                                        cAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAg                        >  2:128445/1‑100 (MQ=255)

                                                                        cAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAg                        <  1:1304296/100‑1 (MQ=255)

                                                                        cAGCGTTTCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAg                        >  1:813777/1‑100 (MQ=255)

                                                                             tttCAAACTGCGTCACCAACATTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGaa                   >  1:1915656/1‑100 (MQ=255)

                                                                             tttCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGaa                   <  1:569152/100‑1 (MQ=255)

                                                                              ttCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGaaa                  <  1:956690/100‑1 (MQ=255)

                                                                              ttCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGaaa                  <  2:928457/100‑1 (MQ=255)

                                                                              ttCAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCCTTCAAATAGGTGaaa                  <  2:1916154/100‑1 (MQ=255)

                                                                                cAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAaga                >  2:454922/1‑100 (MQ=255)

                                                                                cAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAaga                >  1:982179/1‑100 (MQ=255)

                                                                                cAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAaga                >  1:3473114/1‑100 (MQ=255)

                                                                                cAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAaga                <  1:3112208/100‑1 (MQ=255)

                                                                                cAAACTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAaga                >  1:2831531/1‑100 (MQ=255)

                                                                                 aaaCTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAagag               >  2:1728501/1‑100 (MQ=255)

                                                                                 aaaCTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAAACACATCAAAATCGCCTTTCAAATAGATGAAagag               >  2:2625071/1‑100 (MQ=255)

                                                                                  aaCTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAagaga              >  1:1186717/1‑100 (MQ=255)

                                                                                  aaCTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAagaga              <  1:2186342/100‑1 (MQ=255)

                                                                                  aaCTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAagaga              >  1:2744294/1‑100 (MQ=255)

                                                                                  aaCTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAagaga              <  2:1977551/100‑1 (MQ=255)

                                                                                  aaCTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAagaga              >  2:3919036/1‑100 (MQ=255)

                                                                                   aCTGCGTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAagagag             <  1:1874687/100‑1 (MQ=255)

                                                                                        gTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCg        >  2:2929905/1‑100 (MQ=255)

                                                                                        gTCACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCg        >  2:3967392/1‑100 (MQ=255)

                                                                                          cACCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGa      <  2:3160606/100‑1 (MQ=255)

                                                                                           aCCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGAt     >  1:3052758/1‑100 (MQ=255)

                                                                                           aCCAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGAt     <  2:2149926/100‑1 (MQ=255)

                                                                                            ccAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGAtt    <  2:956726/100‑1 (MQ=255)

                                                                                             cAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTa   >  2:2152719/1‑100 (MQ=255)

                                                                                             cAACAGTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTa   >  2:2662826/1‑100 (MQ=255)

TGCCACAACACCTTCCAGCGGTTGATCGATATCCACCACCAGTACATCGGTTTCGTCCGCACCCGGCTCCTGCAGCGTTTCAAACTGCGTCACCAACATTTGGGTTTTAAAGAAATGGCCTTTGCGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTAC  >  NC_000913/3577098‑3577291

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: `‑`