BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,379,402:1	+A	100%	coding (596/1128 nt)	zapE ←	divisome ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,379,402	1	.	A	96.8%	430.7 / 5.8	126	coding (596/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand): ref base . (3/1); new base A (67/55); total (70/56)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 134 total aligned reads displayed.

GTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑T‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAGGGTAA  >  NC_000913/3379303‑3379495

gTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGt‑t                                                                                                >  2:1943442/1‑100 (MQ=255)

gTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGt‑t                                                                                                <  2:1382807/100‑1 (MQ=255)

 tttCATGGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGt‑t                                                                                                <  1:2503947/99‑1 (MQ=255)

 tttCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTTATTACATCACAATGCTGt‑ta                                                                                               <  1:2257498/100‑2 (MQ=255)

  ttCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGt‑ta                                                                                               <  1:1039267/99‑2 (MQ=255)

     aTCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGt‑tataat                                                                                           >  1:2836597/1‑95 (MQ=255)

       cGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATg                                                                                          <  1:3195845/99‑1 (MQ=255)

        gTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGc                                                                                        <  2:3180327/100‑1 (MQ=255)

         tGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCa                                                                                       <  1:2714321/100‑1 (MQ=255)

          gAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCAt                                                                                      >  2:437772/1‑100 (MQ=255)

          gAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCAt                                                                                      >  2:2157052/1‑100 (MQ=255)

           aaGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATc                                                                                     <  1:813705/100‑1 (MQ=255)

             gTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGa                                                                                   <  2:212010/100‑1 (MQ=255)

              tGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGAt                                                                                  <  2:3307114/100‑1 (MQ=255)

               ggCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGAtt                                                                                 >  2:1086755/1‑100 (MQ=255)

               ggCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGAtt                                                                                 >  2:3361461/1‑100 (MQ=255)

               ggCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGAt                                                                                  >  1:1987504/1‑99 (MQ=255)

                 cGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTgc                                                                               >  1:1471967/1‑100 (MQ=255)

                 cGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTgc                                                                               <  1:2675279/100‑1 (MQ=255)

                 cGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTgc                                                                               >  1:2084587/1‑100 (MQ=255)

                  gAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTgca                                                                              <  1:3163874/100‑1 (MQ=255)

                  gAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTgca                                                                              <  2:296251/100‑1 (MQ=255)

                  gAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTgca                                                                              >  1:2994704/1‑100 (MQ=255)

                  gAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTgca                                                                              <  1:2857465/100‑1 (MQ=255)

                   ttcgCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTgcag                                                                             <  1:1737051/97‑1 (MQ=255)

                   aaaGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTgcag                                                                             <  1:949609/100‑1 (MQ=255)

                    aaGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTgcagg                                                                            >  2:507107/1‑100 (MQ=255)

                     aGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTgcaggn                                                                           <  1:1770501/100‑2 (MQ=255)

                     aGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTgcaggc                                                                           <  2:1825355/100‑1 (MQ=255)

                     aGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑T‑TAATGGCATCGATTgcaggca                                                                          <  2:3525149/100‑1 (MQ=255)

                      gCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTgcaggca                                                                          >  2:2033196/1‑100 (MQ=255)

                        cacaGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGa                                                                        >  2:3105937/1‑100 (MQ=255)

                            gATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTGCACGTTCATTACATCACAATGCTGTCA‑TAATGGCATCGATTGCAGGCAGAAAAc                                                                    >  1:3573943/1‑100 (MQ=255)

                            gATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAAc                                                                    >  1:3064002/1‑100 (MQ=255)

                             aTGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCTTCGATTGCAGGCAGAAAAcg                                                                   <  1:2263238/100‑1 (MQ=255)

                             aTGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAAcg                                                                   <  1:3582561/100‑1 (MQ=255)

                              tGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAATcg                                                                   <  2:159256/99‑1 (MQ=255)

                              tGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAAcgc                                                                  >  1:2830829/1‑100 (MQ=255)

                                ngcCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACgcgc                                                                >  1:3223403/2‑100 (MQ=255)

                                cgcCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACgcgc                                                                >  1:2461984/1‑100 (MQ=255)

                                cgcCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACgcgc                                                                >  2:2080398/1‑100 (MQ=255)

                                 gcCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACgcgc                                                                >  1:1460888/1‑99 (MQ=255)

                                  ccTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCAc                                                              >  1:2522551/1‑100 (MQ=255)

                                    tgagtgagTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGt                                                            <  1:497628/100‑1 (MQ=255)

                                      agtgagtgTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTg                                                          >  1:3134333/1‑100 (MQ=255)

                                           gtgtACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGc                                                     >  1:1885697/1‑100 (MQ=255)

                                           gtgtACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGc                                                     <  1:2718028/100‑1 (MQ=255)

                                             gtACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCa                                                   >  1:3345689/1‑100 (MQ=255)

                                             gtACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCa                                                   >  2:1010743/1‑100 (MQ=255)

                                              tACGCAGACGATAAGCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCAt                                                  <  2:1717167/100‑1 (MQ=255)

                                               aCGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCAtt                                                 <  2:3182773/100‑1 (MQ=255)

                                                cGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATACTGGCATCGACTGCAGGCAGAAAACGCGCACGTTGCAGGCCAttt                                                >  1:1305257/1‑100 (MQ=255)

                                                cGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCAttt                                                >  2:744508/1‑100 (MQ=255)

                                                  catccGATAACCAACACCGGCGTCCACGTTCATTACATCACAATGCTTT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCg                                              <  1:1213633/96‑1 (MQ=255)

                                                  cAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCg                                              >  2:2149758/1‑100 (MQ=255)

                                                  cAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCg                                              >  2:176312/1‑100 (MQ=255)

                                                   aGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGa                                             >  2:2547987/1‑100 (MQ=255)

                                                    gACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGAt                                            <  2:2182600/100‑1 (MQ=255)

                                                     aCGATAATCAACATCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATATTGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGAt                                            <  1:1461390/99‑1 (MQ=255)

                                                      cGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATaa                                          >  1:1038496/1‑100 (MQ=255)

                                                        aTAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAg                                        <  1:619554/100‑1 (MQ=255)

                                                         tAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGt                                       >  2:2327786/1‑100 (MQ=255)

                                                          aaTCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGtt                                      >  2:2916067/1‑100 (MQ=255)

                                                           aTCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTc                                     >  2:1688885/1‑100 (MQ=255)

                                                           aTCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTc                                     >  2:330123/1‑100 (MQ=255)

                                                           aTCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTc                                     >  2:2727125/1‑100 (MQ=255)

                                                            tCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCa                                    >  2:2435778/1‑100 (MQ=255)

                                                            tCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCa                                    >  2:2436291/1‑100 (MQ=255)

                                                             cAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCAt                                   <  2:3265304/100‑1 (MQ=255)

                                                             cAACACAAACGTCCACGTTCATTACATCACAATGATGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCa                                    <  1:444331/99‑1 (MQ=255)

                                                               acacCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATcc                                 >  1:2460505/1‑100 (MQ=255)

                                                               acacCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATcc                                 <  1:2541096/100‑1 (MQ=255)

                                                                cacCAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCcg                                >  1:2030988/1‑100 (MQ=255)

                                                                 acCAGCGTCCACGTTCATTACATCACACTTCTGT‑TATAATGCCATCGATTGCAGGCAGAAAACGCGCCCGTTGCAGGCCATTCCGATAAAGTTCATc                                  >  1:643556/1‑97 (MQ=255)

                                                                  ccAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCcggc                              >  2:2266429/1‑100 (MQ=255)

                                                                   cAGCGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCcggcg                             <  2:1450909/100‑1 (MQ=255)

                                                                      cGTCCACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGaa                          >  2:2406193/1‑100 (MQ=255)

                                                                         ccACGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAAAGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGACGGAatat                       >  1:482846/1‑100 (MQ=255)

                                                                            cGTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTg                    >  2:461806/1‑100 (MQ=255)

                                                                             gTTCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGa                   <  1:2547973/100‑1 (MQ=255)

                                                                              ttCATTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTGGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATGTGAc                  <  1:2972861/100‑1 (MQ=255)

                                                                                 atgaCATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGTCCACTTCGATAAACAGCATCCGGCGGAATATTTGAc                  >  1:926644/4‑97 (MQ=255)

                                                                                 aTTACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTc               <  1:330123/100‑1 (MQ=255)

                                                                                  ttACATCACAATGCTGT‑TATAATGGCATCGATTGCAGGAAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCg              <  1:1632742/100‑1 (MQ=255)

                                                                                    aCATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCATGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCt            <  1:1169242/100‑1 (MQ=255)

                                                                                    aCATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCt            >  1:1647804/1‑100 (MQ=255)

                                                                                    aCATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGCCGCt            >  2:822401/1‑100 (MQ=255)

                                                                                    aCATCACAATGCTGT‑T‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAGGTTAATCCGGCGGATTATTTGACGTCGc             >  1:1327832/1‑98 (MQ=255)

                                                                                     cATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTa           >  2:165726/1‑100 (MQ=255)

                                                                                     cATCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTa           >  2:1516790/1‑100 (MQ=255)

                                                                                       tCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTAcc         >  2:1237944/1‑100 (MQ=255)

                                                                                       tCACAATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTAcc         >  1:1414828/1‑100 (MQ=255)

                                                                                         acaATGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAg       <  2:2416918/100‑1 (MQ=255)

                                                                                           aaTGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAggg     <  2:1700133/100‑1 (MQ=255)

                                                                                           aaTGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAggg     <  2:3424938/100‑1 (MQ=255)

                                                                                           aaTGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAggg     >  1:3677042/1‑100 (MQ=255)

                                                                                             tGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAGGGTa   <  2:2608366/100‑1 (MQ=255)

                                                                                             tGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAGGGTa   <  2:3467075/100‑1 (MQ=255)

                                                                                             tGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAGGGTa   <  1:709790/100‑1 (MQ=255)

                                                                                             tGCTGT‑TATAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAGGGTa   <  1:2406180/100‑1 (MQ=255)

GTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGT‑T‑TAATGGCATCGATTGCAGGCAGAAAACGCGCACGTTGCAGGCCATTTCGATAAAGTTCATCCGGCGGAATATTTGACGTCGCTACCAGGGTAA  >  NC_000913/3379303‑3379495

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