BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,323,280	0	A	G	22.6%	205.3 / 58.1	102	intergenic (+29/‑167)	atoC/atoD	fused response regulator of ato operon, in two‑component system with AtoS: response regulator/sigma54 interaction protein/acetyl‑CoA:acetoacetyl‑CoA transferase, alpha subunit
Reads supporting (aligned to +/- strand): ref base A (39/40); new base G (10/13); total (49/53)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-01

Only 100 of 110 total aligned reads displayed.

TGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATC  >  NC_000913/2323181‑2323379

tGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACATCAAAACTTGCTATGCAGAAATTTGca                                                                                                     <  2:1823269/100‑1 (MQ=255)

tGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGca                                                                                                     >  2:912340/1‑100 (MQ=255)

 gCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGcac                                                                                                    >  1:1000838/1‑100 (MQ=255)

            gTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAtt                                                                                         >  2:2303283/1‑100 (MQ=255)

               gcggTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTc                                                                                      <  1:1277319/97‑1 (MQ=255)

               gCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTc                                                                                      <  2:1606692/100‑1 (MQ=37)

                  gTGCATTGATGTATAAGCTCCAGGAATACGGTATCGATCCGGCGAATGTATAACAGCAAAGCTTGTTATGCAGAATATTGCGCACTGTGCAAATTTCTGc                                                                                   <  1:586771/100‑1 (MQ=21)

                  gTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTg                                                                                    >  2:2924135/1‑99 (MQ=255)

                  gTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGc                                                                                   <  2:1630651/100‑1 (MQ=255)

                  gTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGc                                                                                   <  2:826102/100‑1 (MQ=25)

                  gTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGc                                                                                   <  1:2654761/100‑1 (MQ=25)

                    gCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCAt                                                                                 <  2:138246/100‑1 (MQ=255)

                    gCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCAt                                                                                 <  1:566046/100‑1 (MQ=25)

                    gCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCAt                                                                                 <  1:428786/100‑1 (MQ=25)

                      aTTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAg                                                                               <  2:1617389/100‑1 (MQ=25)

                      aTTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAg                                                                               <  2:1863744/100‑1 (MQ=25)

                      aTTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAg                                                                               <  1:950845/100‑1 (MQ=25)

                       ttGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGc                                                                              >  1:100442/1‑100 (MQ=255)

                       ttGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAgn                                                                              <  2:686290/100‑2 (MQ=25)

                       ttGATGTATAAACTCCAGGAATACGGTAGCGATCCGGCGGATGTATAACACCAAAACTTTCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGc                                                                              <  2:1310086/100‑1 (MQ=21)

                        tGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGca                                                                             <  1:1684215/100‑2 (MQ=25)

                         gATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCg                                                                            >  1:1075515/1‑100 (MQ=255)

                           tGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCt                                                                          <  2:2528816/100‑1 (MQ=255)

                           ngTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCCATCTTCTGCATacaagct                                                                          >  2:2049887/2‑94 (MQ=255)

                             nataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCtaa                                                                        >  2:2701369/2‑98 (MQ=255)

                              ataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAACCTTGCTATGCCGACATTTGCACAATGCGCAATTTTCTGCATAGCCGCTCAt                                                                       >  1:583995/1‑100 (MQ=255)

                               taAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCAtt                                                                      >  1:1323799/1‑100 (MQ=255)

                                aaaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtc                                                                     <  2:1937763/100‑1 (MQ=255)

                                 aaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGAGCGCAATTTTCTGCATAGCCGCTCATtct                                                                    >  1:1486210/1‑100 (MQ=255)

                                  aCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtctc                                                                   <  2:2865292/100‑1 (MQ=255)

                                   cTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTcc                                                                  <  2:1010175/100‑1 (MQ=255)

                                   cTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTcc                                                                  >  2:842327/1‑100 (MQ=255)

                                    tCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCt                                                                 <  1:1998759/100‑1 (MQ=255)

                                    tCCAGGAAGACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCt                                                                 <  1:2924786/100‑1 (MQ=255)

                                     ccAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCtt                                                                >  2:1463343/1‑100 (MQ=255)

                                     ccAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCtt                                                                >  1:190750/1‑100 (MQ=255)

                                       aGGAATACGGTATCGATCCGGCGGATGTATAACACCACAACTTGCTATGCAGAAATTTGCACAGTGCGCATTTTTCTGCATAGCCGCTCATTCTCCTtat                                                              <  2:190751/100‑1 (MQ=255)

                                        ggAATACGGTATCGATCCGGCGGATGTATAACACCAAGACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTtata                                                             >  1:1597301/1‑100 (MQ=255)

                                         gAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATaa                                                            <  2:2647399/100‑1 (MQ=255)

                                            tACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTACAAAt                                                          <  1:1838795/99‑1 (MQ=255)

                                              cGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATccc                                                       <  2:304937/100‑1 (MQ=255)

                                              cGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATccc                                                       >  2:1190382/1‑100 (MQ=255)

                                              cGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATccc                                                       <  2:1587992/100‑1 (MQ=255)

                                                ggaTCGATCCGTCGGATGTATAACACCAAAACTTGCTATGCAGAGATTTGCATAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATccc                                                       <  1:1737595/96‑1 (MQ=255)

                                                 tATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCAt                                                    >  1:1373100/1‑100 (MQ=255)

                                                  aTCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCaac                                                   >  1:1820021/1‑98 (MQ=255)

                                                  aTCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATc                                                   <  1:1662568/100‑1 (MQ=255)

                                                     gATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCaa                                                >  1:1202893/1‑100 (MQ=255)

                                                     gATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCaa                                                >  2:2185034/1‑100 (MQ=255)

                                                      aTCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAAt                                               <  2:1351277/100‑1 (MQ=255)

                                                        ccGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAAttt                                             >  2:638214/1‑100 (MQ=255)

                                                          ggcggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTAt                                           >  1:1172380/1‑100 (MQ=255)

                                                            cggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATcc                                         >  2:1985244/1‑100 (MQ=255)

                                                            cggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATcc                                         >  1:1714598/1‑100 (MQ=255)

                                                            cggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATcc                                         <  1:1248640/100‑1 (MQ=255)

                                                             ggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATccc                                        <  2:1000843/100‑1 (MQ=255)

                                                                tGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTc                                     >  2:2510736/1‑100 (MQ=255)

                                                                 gTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCa                                    <  2:1820029/100‑1 (MQ=255)

                                                                 gTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCa                                    <  2:388196/100‑1 (MQ=255)

                                                                  tataACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCat                                   <  2:1323805/100‑1 (MQ=255)

                                                                  tataACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCat                                   >  2:2188051/1‑100 (MQ=255)

                                                                   ataACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCata                                  >  2:2655244/1‑100 (MQ=255)

                                                                   ataACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCata                                  >  2:742539/1‑100 (MQ=255)

                                                                    taACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCatat                                 >  1:531054/1‑100 (MQ=255)

                                                                    taACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCatat                                 >  2:2530128/1‑100 (MQ=255)

                                                                     aaCACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATAtt                                >  2:639603/1‑100 (MQ=255)

                                                                       cacCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCa                              <  1:2598829/100‑1 (MQ=255)

                                                                         ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                            >  1:2362900/1‑100 (MQ=255)

                                                                         ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                            >  2:2715337/1‑100 (MQ=255)

                                                                         ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                            <  2:1373106/100‑1 (MQ=255)

                                                                         ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                            >  2:378278/1‑100 (MQ=255)

                                                                           aaaaCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAttn                          <  2:2751132/100‑2 (MQ=255)

                                                                           aaaaCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAtt                           <  1:912336/99‑1 (MQ=255)

                                                                           aaaaCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTa                          <  2:2084214/100‑1 (MQ=255)

                                                                            aaaCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAg                         <  2:173685/100‑1 (MQ=255)

                                                                             aaCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGt                        <  2:2486136/100‑1 (MQ=255)

                                                                              ttccgatctGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGtt                       <  1:2131412/92‑1 (MQ=255)

                                                                               gcggCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTa                      >  2:1113984/4‑100 (MQ=21)

                                                                               cTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTa                      >  2:2859494/1‑100 (MQ=255)

                                                                                 tGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTaaa                    <  1:2957560/100‑1 (MQ=255)

                                                                                  gCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAAt                   <  2:1532247/100‑1 (MQ=255)

                                                                                   cTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATa                  <  2:1771857/100‑1 (MQ=255)

                                                                                   cTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATa                  >  2:55839/1‑100 (MQ=25)

                                                                                    tATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATaa                 <  2:1595727/100‑1 (MQ=255)

                                                                                     aTGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAAc                <  1:2185025/100‑1 (MQ=255)

                                                                                     aTGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAAc                >  2:2971097/1‑100 (MQ=25)

                                                                                       gCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTagatcggaa              >  2:2131421/1‑92 (MQ=255)

                                                                                       gCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa              >  2:2430589/1‑100 (MQ=255)

                                                                                       gCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa              >  2:1879039/1‑100 (MQ=255)

                                                                                       gCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa              >  1:780477/1‑100 (MQ=25)

                                                                                        cAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa             >  2:1675382/1‑100 (MQ=37)

                                                                                          gAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAAt           >  1:1897691/1‑100 (MQ=37)

                                                                                          gAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAAt           >  1:1897568/1‑100 (MQ=37)

                                                                                            aaTTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATcc         <  2:1178687/100‑1 (MQ=255)

                                                                                            aaTTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATcc         <  1:2942528/100‑1 (MQ=255)

                                                                                             gatctcacaGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCa        <  1:226712/95‑1 (MQ=255)

                                                                                               ttGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaat      >  2:638375/1‑100 (MQ=255)

                                                                                                  cacaGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaagatc   >  2:226713/1‑96 (MQ=255)

                                                                                                   acaGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAAtc  <  2:2648076/100‑1 (MQ=255)

                                                                                                   acaGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAAtc  <  2:882148/100‑1 (MQ=255)

TGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATC  >  NC_000913/2323181‑2323379

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