BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,093,283	G→A	G116D (GGC→GAC)	yqgE →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,093,283	0	G	A	99.1%	345.8 / NA	106	G116D (GGC→GAC)	yqgE	uncharacterized protein
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (50/55); total (50/56)
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 106 total aligned reads displayed.

TATTTTGCATACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAACGCGTGGC  >  NC_000913/3093188‑3093378

tATTTTGCATACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACAcc                                                                                             >  2:2240866/1‑100 (MQ=255)

 aTTTTGCATACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCg                                                                                            <  1:2509061/100‑1 (MQ=255)

   tttGCATCCTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGAt                                                                                          <  2:94878/100‑1 (MQ=255)

   tttGCATACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGa                                                                                           <  2:398978/99‑1 (MQ=255)

   tttGCATACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGAt                                                                                          <  1:1463600/100‑1 (MQ=255)

   tttGCATACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGAt                                                                                          <  2:2971356/100‑1 (MQ=255)

    ttGCATACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATa                                                                                         >  2:2863675/1‑100 (MQ=255)

     tGCATACTCCGCCATCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCCATGTGCTGGAAACGCTCGACACCGAt                                                                                          >  2:995267/1‑98 (MQ=255)

        aTACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAaaca                                                                                     <  1:2057186/100‑1 (MQ=255)

        aTACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACCTCCCGCGATGTGCTGGAAACGCTCTACACCGATAaaca                                                                                     <  2:1901088/100‑1 (MQ=255)

         tACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAaacaa                                                                                    >  2:2336751/1‑100 (MQ=255)

          aCTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAaacaac                                                                                   >  2:505781/1‑100 (MQ=255)

          aCTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAaacaac                                                                                   >  2:680211/1‑100 (MQ=255)

          aCTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAaacaac                                                                                   >  2:2107873/1‑100 (MQ=255)

           cTCCGCCCTCCAATTTTGCTTCCAGCATTCGCCTTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAAcc                                                                                  >  2:2454439/1‑100 (MQ=255)

           cTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCGGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAAcc                                                                                  >  1:202637/1‑100 (MQ=255)

               gCCCTCCAATTTTGGTTCCAGCAGTCGTCTTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGGTGGAAACGCTCGACACCGATAaacaac                                                                                   <  2:2797958/95‑1 (MQ=255)

                 ccTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGa                                                                            >  2:954309/1‑100 (MQ=255)

                 ccTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGa                                                                            >  2:959821/1‑100 (MQ=255)

                  cTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGAc                                                                           <  2:2717776/100‑1 (MQ=255)

                     cAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGt                                                                         >  1:2412239/1‑99 (MQ=255)

                     cAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTa                                                                        >  2:322994/1‑100 (MQ=255)

                         tttGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATtgg                                                                    <  1:777848/100‑1 (MQ=255)

                          ttGCTTCCAGCATTCGCTTTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATtggn                                                                   <  2:1695476/100‑2 (MQ=255)

                          ttGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATtggt                                                                   <  1:1561628/100‑1 (MQ=255)

                          ttGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATtggt                                                                   <  1:1237618/100‑1 (MQ=255)

                             cTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATtgg                                                                    >  2:1295826/1‑96 (MQ=255)

                               tCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGctc                                                              >  2:3090554/1‑100 (MQ=255)

                               tCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGctc                                                              <  2:807707/100‑1 (MQ=255)

                                ccAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCATCTGACGTATTGGTGGctc                                                              >  2:1800498/1‑99 (MQ=255)

                                       tCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTAt                                                      <  1:1790473/100‑1 (MQ=255)

                                        cGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATg                                                     <  2:641320/100‑1 (MQ=255)

                                        cGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATg                                                     >  1:2744140/1‑100 (MQ=255)

                                        cGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATg                                                     <  2:1247994/100‑1 (MQ=255)

                                        cGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATg                                                     >  2:1227411/1‑100 (MQ=255)

                                         gCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGc                                                    >  1:726542/1‑100 (MQ=255)

                                         gCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGc                                                    <  1:1005192/100‑1 (MQ=255)

                                         gCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGc                                                    >  1:2584059/1‑100 (MQ=255)

                                           aTTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGcct                                                  >  2:2902704/1‑100 (MQ=255)

                                            tttCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGcctc                                                 >  1:2769076/1‑100 (MQ=255)

                                            tttCAGACAACACGGTAATGACAACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACAGTCTGACGTATTGGTGGCACTGGGCTATGccgc                                                 >  2:1152137/1‑98 (MQ=255)

                                             ttCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGcctcc                                                >  1:2896891/1‑100 (MQ=255)

                                              tCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGCTTATGcc                                                   >  2:2243586/1‑96 (MQ=255)

                                               cAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTg                                              >  2:1829557/1‑100 (MQ=255)

                                               cAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTg                                              <  1:2720456/100‑1 (MQ=255)

                                                aGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTgg                                             <  2:1165800/100‑1 (MQ=255)

                                                 gACAACACGGTAATGACCCCTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTggg                                            <  2:3087418/100‑1 (MQ=255)

                                                 gACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTggg                                            <  2:2373388/100‑1 (MQ=255)

                                                 gACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTgg                                             <  2:1503293/99‑1 (MQ=255)

                                                   caacaCGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGgag                                          <  1:2336774/100‑1 (MQ=255)

                                                     acacGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGaa                                        >  2:2045526/1‑100 (MQ=255)

                                                     acacGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGaa                                        >  1:1646761/1‑100 (MQ=255)

                                                     acacGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGaa                                        >  1:2643438/1‑100 (MQ=255)

                                                      cacGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGaaa                                       >  2:2995126/1‑100 (MQ=255)

                                                      cacGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGaaa                                       <  2:559238/100‑1 (MQ=255)

                                                       acGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAg                                      <  1:302343/100‑1 (MQ=255)

                                                        ccgTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGCTAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAgg                                     <  2:802750/98‑1 (MQ=255)

                                                        cGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAgg                                     <  1:2296165/100‑1 (MQ=255)

                                                        cGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAgg                                     <  2:3085448/100‑1 (MQ=255)

                                                        cGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAg                                      <  1:524382/99‑1 (MQ=255)

                                                         ggTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGt                                    <  1:1800517/100‑1 (MQ=255)

                                                         ggTAATCACCGCTTCCCGCGATTTGCTTGAAACGCGGGACACCGATAAACAACCGTCTGACGTATGGGTGGCTCTGGGTTATGCCTCCTGGGAGaaa                                       <  2:2665274/97‑1 (MQ=255)

                                                          gTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTc                                   >  1:2440487/1‑100 (MQ=255)

                                                          gTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTc                                   >  1:1727338/1‑100 (MQ=255)

                                                              tGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACt                               >  2:764972/1‑100 (MQ=255)

                                                              tGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACt                               <  2:1331609/100‑1 (MQ=255)

                                                              tGACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACt                               <  1:97949/100‑1 (MQ=255)

                                                               gACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGCAAGGTCAACTg                              >  2:335477/1‑100 (MQ=255)

                                                               gACCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGGGGCTCTGGGTTATGCCTCCTGGGCGCAACGTCAACTg                              >  2:62983/1‑100 (MQ=255)

                                                                aCCACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATCGGTGGCTCTGGGTTATGCCTCCAGGGAGAATGGTCAACTgg                             >  2:602401/1‑100 (MQ=255)

                                                                 ccACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGa                            <  1:111785/100‑1 (MQ=255)

                                                                  cACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGaa                           >  1:680340/1‑100 (MQ=255)

                                                                  cACTTCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGaa                           >  2:1289326/1‑100 (MQ=255)

                                                                     ttCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACaa                        <  2:2769049/100‑1 (MQ=255)

                                                                      tCCCGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAg                       >  1:2730942/1‑100 (MQ=255)

                                                                       cccGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAg                       <  2:2412216/99‑1 (MQ=255)

                                                                       cccGCGATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGa                      >  2:2929942/1‑100 (MQ=255)

                                                                       cccGCGATGTGCTGGAAACGCTCGACACAGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAg                       <  1:938511/99‑1 (MQ=255)

                                                                          gcgATGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAAt                   <  1:3083367/100‑1 (MQ=255)

                                                                             aTGTGCTGGAAACGCTCGACACCGATATACAACCGTCTGACGTATTGGAGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATtct                <  2:818699/100‑1 (MQ=255)

                                                                             aTGTGCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATtct                <  2:2916787/100‑1 (MQ=255)

                                                                              tgtgCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATtctc               >  1:1581241/1‑100 (MQ=255)

                                                                              tgtgCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATtctc               <  1:554642/100‑1 (MQ=255)

                                                                               gtgCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCg              <  2:773916/100‑1 (MQ=255)

                                                                               gtgCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCg              <  2:521902/100‑1 (MQ=255)

                                                                                tgCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGa             <  2:2475538/100‑1 (MQ=255)

                                                                                tgCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGa             >  1:2454423/1‑100 (MQ=255)

                                                                                 gCTGGAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGAt            <  2:2523310/100‑1 (MQ=255)

                                                                                    ggAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAAc         <  1:2512526/100‑1 (MQ=255)

                                                                                     gAAACGCTCGACACCGATAAACAACCGTCTGACGTATTTGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAAcg        <  2:1681267/100‑1 (MQ=255)

                                                                                     gAAACGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAAcg        >  1:562091/1‑100 (MQ=255)

                                                                                      aaaCGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAAcgc       <  1:2488675/100‑1 (MQ=255)

                                                                                       aaCGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAAcgcg      >  2:2561794/1‑100 (MQ=255)

                                                                                        aCGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAACGCGt     <  2:2903704/100‑1 (MQ=255)

                                                                                         cGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAACGCGTg    <  1:1767204/100‑1 (MQ=255)

                                                                                         cGCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGACAACGCGTg    >  1:1761330/1‑100 (MQ=255)

                                                                                          gCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAACGCGTgg   >  1:3080220/1‑100 (MQ=255)

                                                                                          gCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAACGCGTgg   <  1:2392151/100‑1 (MQ=255)

                                                                                          gCTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAACGCGTgg   >  2:2829490/1‑100 (MQ=255)

                                                                                           cTCGACACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAACGCGTGGc  <  1:2105583/100‑1 (MQ=255)

TATTTTGCATACTCCGCCCTCCAATTTTGCTTCCAGCATTCGCATTTCAGACAACACGGTAATGACCACTTCCCGCGATGTGCTGGAAACGCTCGGCACCGATAAACAACCGTCTGACGTATTGGTGGCTCTGGGTTATGCCTCCTGGGAGAAAGGTCAACTGGAACAAGAAATTCTCGATAACGCGTGGC  >  NC_000913/3093188‑3093378

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