BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,181,294	G→C	15.1%	D292H (GAT→CAT)	ygiC →	ATP‑Grasp family ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,181,294	0	G	C	15.1%	501.5 / 75.3	180	D292H (GAT→CAT)	ygiC	ATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand): ref base G (82/70); new base C (15/12); total (97/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.39e-01

Only 100 of 194 total aligned reads displayed.

GCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGG  >  NC_000913/3181195‑3181393

gCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgatg                                                                                                     <  1:711907/100‑1 (MQ=255)

 ngTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgatga                                                                                                    >  1:3433776/2‑100 (MQ=255)

   tGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatc                                                                                                  <  2:2356313/100‑1 (MQ=255)

   tGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatc                                                                                                  <  2:1415350/100‑1 (MQ=255)

    ggAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatca                                                                                                 <  1:3698509/100‑1 (MQ=255)

     gAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatcat                                                                                                >  1:4093271/1‑100 (MQ=255)

      aaGAGCATTATCTCCAACACGGCACTTGTACCGATCCTGTTTGTGACGTCACCGAACCCCATGAACATGCAGCACGCTTATTTTGCAGAAGATGatcat                                                                                                >  2:3620276/1‑99 (MQ=255)

          gCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCa                                                                                           <  1:116477/100‑1 (MQ=255)

          gCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCa                                                                                           >  1:529606/1‑100 (MQ=255)

          gCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTACTGCCCGCTTATTTTGCGGATGATCATCTTCCGCa                                                                                           <  1:522742/100‑1 (MQ=255)

           cATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCaa                                                                                          <  2:3897338/100‑1 (MQ=255)

           cATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCaa                                                                                          >  1:3770246/1‑100 (MQ=255)

            aTTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaa                                                                                         <  1:2350479/100‑1 (MQ=255)

            aTTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCaaa                                                                                         <  1:467548/100‑1 (MQ=255)

            aTTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCaaa                                                                                         >  2:1951731/1‑100 (MQ=255)

            aTTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCaaa                                                                                         >  1:1422127/1‑100 (MQ=255)

              tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaat                                                                                       <  1:1148263/100‑3 (MQ=38)

              tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATg                                                                                       <  1:501604/100‑1 (MQ=255)

              tATCGCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaat                                                                                       <  2:1624553/100‑3 (MQ=38)

               aTCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaata                                                                                      <  1:3724567/100‑4 (MQ=38)

                 ctcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaa                                                                                    >  2:3262213/1‑100 (MQ=255)

                  tcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaaa                                                                                   >  2:3446468/1‑100 (MQ=255)

                  tcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaaa                                                                                   <  2:696996/100‑1 (MQ=255)

                  tcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGCTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagc                                                                                   <  1:2922394/100‑7 (MQ=37)

                   ccAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaaaa                                                                                  >  2:993107/1‑100 (MQ=255)

                   ccAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaaaa                                                                                  >  2:1877279/1‑100 (MQ=255)

                   ccAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTCGCGGAAGATGATCATCCGCAAATGGaaaa                                                                                  >  2:787024/1‑100 (MQ=255)

                      acaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAata                                                                               >  1:3226083/1‑100 (MQ=255)

                       caaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAatat                                                                              >  1:1764080/1‑100 (MQ=255)

                          ggCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtgtg                                                                           <  1:453520/100‑1 (MQ=255)

                            cACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGt                                                                         <  1:3547627/100‑1 (MQ=255)

                             aCTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGtt                                                                        >  1:2145633/1‑100 (MQ=255)

                               ttCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTaa                                                                      >  2:531507/1‑100 (MQ=255)

                                  tACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAcc                                                                   <  1:2324695/100‑1 (MQ=255)

                                      gCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtc                                                               <  1:3040073/100‑1 (MQ=255)

                                       cTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtct                                                              <  1:1427812/100‑1 (MQ=255)

                                        tACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtctt                                                             >  1:416020/1‑100 (MQ=255)

                                            gtgGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTcc                                                         <  1:572679/100‑1 (MQ=255)

                                                gagaTGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTg                                                     <  2:3813729/100‑1 (MQ=255)

                                                gagaTGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTg                                                     >  1:3630532/1‑100 (MQ=255)

                                                gagaTGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTg                                                     >  1:3064925/1‑100 (MQ=255)

                                                gagaTGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATGTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTg                                                     <  2:944297/100‑1 (MQ=255)

                                                 tattGTTCCAGAATCACCCGAACCTGCTGCCCGCTTATTTTCCGGAAAATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGt                                                      <  2:129644/95‑1 (MQ=255)

                                                    tGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAACTGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAg                                                  >  2:1735448/1‑99 (MQ=255)

                                                     gTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgc                                                >  1:2278220/1‑100 (MQ=255)

                                                      ttCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgcg                                               <  1:1945160/100‑1 (MQ=255)

                                                        cccGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCa                                             >  1:943946/1‑100 (MQ=255)

                                                        cccGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCa                                             >  1:253657/1‑100 (MQ=255)

                                                        cccGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCa                                             >  2:3560431/1‑100 (MQ=255)

                                                        cccGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGAACATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCa                                             <  2:1334399/100‑1 (MQ=255)

                                                          cGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCaaa                                           >  2:3142764/1‑100 (MQ=255)

                                                           gAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAAc                                          <  1:2192548/100‑1 (MQ=255)

                                                            aaTCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACg                                         <  1:86178/100‑1 (MQ=255)

                                                            aaTCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACg                                         >  1:3199295/1‑100 (MQ=255)

                                                            aaTCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACg                                         >  2:3221358/1‑100 (MQ=255)

                                                             aTCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACgt                                        <  1:1874412/100‑1 (MQ=255)

                                                              tCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACgtg                                       >  1:1476029/1‑100 (MQ=255)

                                                               cACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACgtgt                                      >  2:2906995/1‑100 (MQ=255)

                                                                aCCCGAACCTGCTGCCCGCTTATTTTGCGGTAGCTGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTc                                     >  2:2738059/1‑100 (MQ=255)

                                                                aCCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTc                                     <  2:861561/100‑1 (MQ=255)

                                                                aCCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTc                                     >  2:2649485/1‑100 (MQ=255)

                                                                aCCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTc                                     <  2:2697271/100‑1 (MQ=255)

                                                                aCCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTc                                     >  2:3182164/1‑100 (MQ=255)

                                                                aCCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTc                                     >  2:626445/1‑100 (MQ=255)

                                                                 cccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCg                                    >  2:1623674/1‑100 (MQ=255)

                                                                        ctgctgCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTg                             <  1:3313513/100‑1 (MQ=255)

                                                                         tgctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTgn                            <  1:307855/100‑2 (MQ=255)

                                                                          gctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTgag                           >  1:579732/1‑100 (MQ=255)

                                                                           ctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTgaga                          >  2:1394699/1‑100 (MQ=255)

                                                                           ctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTgaga                          >  2:2292039/1‑100 (MQ=255)

                                                                             gcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAAc                        >  1:1829114/1‑100 (MQ=255)

                                                                             gcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAAc                        >  1:3127864/1‑100 (MQ=255)

                                                                              cccGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACg                       >  1:1676491/1‑100 (MQ=255)

                                                                               ncGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGATCgg                      >  1:2354441/2‑100 (MQ=255)

                                                                               ccGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACgg                      <  2:73619/100‑1 (MQ=255)

                                                                                cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                     >  2:3308289/1‑100 (MQ=255)

                                                                                cGCTTATTTTGCGGAAGATGATCAGCCGCAAATGGAAAAATATGTGGTTAAACCGATCTGCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                     <  1:2437194/100‑1 (MQ=255)

                                                                                     catttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc                >  1:865506/3‑100 (MQ=38)

                                                                                     catttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc                >  1:672353/3‑100 (MQ=38)

                                                                                     catttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc                >  1:1276348/3‑100 (MQ=38)

                                                                                     catttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc                >  1:61106/3‑100 (MQ=38)

                                                                                     catttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc                >  1:161357/3‑100 (MQ=38)

                                                                                     catttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc                >  1:1910562/3‑100 (MQ=38)

                                                                                     aTTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCaaaa                 >  2:2977247/1‑99 (MQ=255)

                                                                                      ttttGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAcc               <  2:4093239/100‑1 (MQ=255)

                                                                                      atttGCGGATGATCGTCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAcc               >  1:3571267/2‑100 (MQ=255)

                                                                                      atttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAcc               >  1:1494658/2‑100 (MQ=39)

                                                                                         tGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCAtt            <  1:132740/100‑1 (MQ=255)

                                                                                         tGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCAtt            <  1:1099050/100‑1 (MQ=255)

                                                                                           cGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGa          >  1:3012124/1‑100 (MQ=255)

                                                                                              aaGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAagc       <  1:4000011/100‑1 (MQ=255)

                                                                                              aaGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAagc       <  2:2511674/100‑1 (MQ=255)

                                                                                               aGATGATCATCCGCAAATGGCAAAATATGTGGTTAAACCGATCTTCTCCCGGGAAGGCGCAAACGTGTCGATCATTGAGCACGGCAAAACCATTg           >  2:2439393/1‑95 (MQ=255)

                                                                                               aGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAagca      >  2:3853833/1‑100 (MQ=255)

                                                                                                gatgatCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAagcag     <  1:2034004/100‑1 (MQ=255)

                                                                                                 atgatCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAagcagc    <  1:2302584/100‑1 (MQ=255)

                                                                                                  tgatCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCg   <  1:2647342/100‑1 (MQ=255)

                                                                                                  tgatCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGGCGATCATTGAGAACGGCAAAACCATTGAAGCAGCg   >  1:1963780/1‑100 (MQ=255)

                                                                                                  tgatCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGAAGCAAACGAGACAAAAATTGAGAACGGCAAAACAATTGAAGCAGCg   >  1:2119753/1‑100 (MQ=255)

                                                                                                   gatCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCgg  >  1:1459363/1‑100 (MQ=255)

GCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGG  >  NC_000913/3181195‑3181393

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