BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,503,169	(G)_5→4	100%	coding (263/1329 nt)	insG ←	IS4 transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,503,165	0	G	.	79.8%	262.3 / 83.8	99	coding (267/1329 nt)	insG	IS4 transposase
Reads supporting (aligned to +/- strand): ref base G (13/7); new base . (36/43); total (49/50)
Fisher's exact test for biased strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

Only 100 of 111 total aligned reads displayed.

GATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCCAACAATAC  >  NC_000913/4503068‑4503262

gaggcggcgTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTggg                                                                                                 <  2:211366/97‑1 (MQ=255)

gATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTggg                                                                                                 <  1:2118280/100‑1 (MQ=255)

 atctggcgtggcgtTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTgggg                                                                                                <  2:1822148/96‑1 (MQ=255)

 aTGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTgggg                                                                                                <  1:2188558/100‑1 (MQ=255)

 aTGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTgggg                                                                                                <  1:2091466/100‑1 (MQ=255)

 aTGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTgggg                                                                                                >  1:2082556/1‑100 (MQ=255)

  tGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTggggg                                                                                               <  1:382761/100‑1 (MQ=255)

  tGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTggggc                                                                                               <  1:296930/100‑2 (MQ=255)

   gcggcgTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTggggca                                                                                              >  2:1078962/1‑98 (MQ=255)

     ggcgtggcgtTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTggggcaac                                                                                            <  1:407382/100‑5 (MQ=255)

     ggcgtggcgtTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTggggcaac                                                                                            <  1:2345763/100‑5 (MQ=255)

     ggcgtggcgtTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTggggagat                                                                                            >  1:1822138/1‑96 (MQ=255)

     ggcgtggcgtTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCaa                                                                                            >  2:1455126/1‑100 (MQ=255)

      gcgtggcgtTATGCCAGAGCTGCGCTGTTTTCGTGAAAACGAGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGAGCAAc                                                                                           >  1:2027289/1‑100 (MQ=255)

       cgtggcgtTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACg                                                                                          <  1:1341674/100‑1 (MQ=255)

            ngtTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAggg                                                                                     >  2:308070/2‑100 (MQ=255)

               tATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGa                                                                                 >  2:2339150/1‑100 (MQ=255)

                aTGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGAt                                                                                <  2:1737065/100‑1 (MQ=255)

                   ccAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGc                                                                             >  2:96862/1‑100 (MQ=255)

                    cAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGcc                                                                            <  1:24003/100‑1 (MQ=255)

                        gCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGgc                                                                        <  2:2416459/100‑1 (MQ=255)

                        gCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGgc                                                                        <  2:2416516/100‑1 (MQ=255)

                        gCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGgc                                                                        <  1:629205/100‑1 (MQ=255)

                          tGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGgcag                                                                      <  1:1834889/100‑1 (MQ=255)

                          tGCGCTGTTTTCATGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGTCACT‑GGGGCAACGAAGGGGCGATTGCCCGgcag                                                                      >  1:1061666/1‑100 (MQ=255)

                            cgcTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGgcagc                                                                     >  1:984570/1‑100 (MQ=255)

                             gcTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAATGGGCGATTGCCCGGCAGCat                                                                   <  1:1705679/100‑1 (MQ=255)

                                gTTTTCGGGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCatgat                                                                <  2:2184616/100‑1 (MQ=255)

                                 ttttCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCatgatg                                                               >  1:437725/1‑100 (MQ=255)

                                   ttCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTc                                                             >  1:1730275/1‑100 (MQ=255)

                                   ttCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTc                                                             >  1:56132/1‑100 (MQ=255)

                                   ttCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTc                                                             >  1:1127411/1‑100 (MQ=255)

                                     cGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCa                                                           >  2:487019/1‑100 (MQ=255)

                                        gAACACGCGGCGTACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGc                                                        >  1:2343935/1‑100 (MQ=255)

                                        gAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAgg                                                         >  2:2317574/1‑100 (MQ=255)

                                        gAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGc                                                        <  2:5977/100‑1 (MQ=255)

                                        gAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGc                                                        <  2:1071615/100‑1 (MQ=255)

                                         aaCACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCg                                                       <  1:1039014/100‑1 (MQ=255)

                                         aaCACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCg                                                       <  1:2114198/100‑1 (MQ=255)

                                            acGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGTTGTCCAGGCGAt                                                     <  2:1696774/100‑1 (MQ=255)

                                             cgcgGCGGACAGCCTCACTTCCCAGGAGCTGGCGGGCCGGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGAtt                                                    >  1:1232509/1‑99 (MQ=255)

                                              gcggcgGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTca                                                  <  1:1936873/100‑1 (MQ=255)

                                                ggcggACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTcaca                                                <  2:1334112/100‑1 (MQ=255)

                                                 gcggACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACaa                                               <  2:1852429/100‑1 (MQ=255)

                                                    gACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGCGCGATTCCCCGGCAGCATGATGTCCAGGCGATTCACaagc                                             >  1:2236819/1‑98 (MQ=255)

                                                    gACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCt                                            <  1:589236/100‑1 (MQ=255)

                                                     aCAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTg                                           >  2:1313428/1‑100 (MQ=255)

                                                     aCAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTg                                           >  2:1683452/1‑100 (MQ=255)

                                                      cAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTgg                                          >  2:2318809/1‑100 (MQ=255)

                                                       aGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGt                                         >  2:1915653/1‑100 (MQ=255)

                                                        gCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGt                                         >  2:719542/1‑100 (MQ=255)

                                                        gCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGt                                         >  1:1662706/1‑100 (MQ=255)

                                                        gCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTg                                        >  1:1046950/1‑100 (MQ=255)

                                                        gCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTg                                        <  1:2203861/100‑1 (MQ=255)

                                                         ccTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGCTTGCCCGGCAGGATGATGTCCAGGCGATTCACAATCTGGTGa                                       >  1:2211620/1‑100 (MQ=255)

                                                         ccTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGa                                       >  1:2211205/1‑100 (MQ=255)

                                                           tCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGaa                                      <  1:92445/100‑1 (MQ=255)

                                                           tCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAag                                     <  2:500989/100‑1 (MQ=255)

                                                           tCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAag                                     <  1:2490511/100‑1 (MQ=255)

                                                            cACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAaga                                    <  1:990399/100‑1 (MQ=255)

                                                             aCTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCACCGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAagag                                   <  1:1812563/100‑1 (MQ=255)

                                                             aCTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAagag                                   <  2:1738037/100‑1 (MQ=255)

                                                             aCTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAagag                                   <  2:1032476/100‑1 (MQ=255)

                                                              cTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAgg                                  >  2:2361327/1‑100 (MQ=255)

                                                              cTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAgg                                  <  2:1787569/100‑1 (MQ=255)

                                                                  ccAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCt                              >  2:1243497/1‑100 (MQ=255)

                                                                   cAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCtt                             >  2:105695/1‑100 (MQ=255)

                                                                   cAGGCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACAAAGGGGCGATTGCCCGGCAGCATGATGTCTAGGCGATTCACAATCTGGTTAAGAGGTtttt                             >  2:259946/1‑100 (MQ=255)

                                                                     ggCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACGATCTGGTGAAGAGGTGCTTTa                           >  2:249767/1‑100 (MQ=255)

                                                                     ggCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTa                           <  2:2438182/100‑1 (MQ=255)

                                                                     ggCGCTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTa                           <  2:263798/100‑1 (MQ=255)

                                                                      gcgcTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTAc                          <  1:706157/100‑1 (MQ=255)

                                                                        gcTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTAcac                        <  1:454509/100‑3 (MQ=255)

                                                                         cTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGc                        <  1:2574299/99‑1 (MQ=255)

                                                                         cTGGCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCt                       <  2:1096223/100‑1 (MQ=255)

                                                                           ggCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTc                      >  2:631505/1‑100 (MQ=255)

                                                                           ggCGGGCCTGAATAACGGCACTGGGGG‑AACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCa                     >  2:1580761/1‑100 (MQ=255)

                                                                           ggCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCa                     >  2:225839/1‑100 (MQ=255)

                                                                           ggCGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCa                     >  2:1582169/1‑100 (MQ=255)

                                                                            gCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCa                     <  2:495442/100‑1 (MQ=255)

                                                                             cGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCATGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAg                   >  1:271824/1‑100 (MQ=255)

                                                                             cGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAg                   >  1:1719757/1‑100 (MQ=255)

                                                                             cGGGCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAg                   >  2:1648664/1‑100 (MQ=255)

                                                                              gggCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAg                   <  2:218770/100‑1 (MQ=255)

                                                                               ggCCTGAATAACGGCACT‑GGGGCAAGGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAgcg                 <  1:1247040/100‑1 (MQ=255)

                                                                                gCCTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAgcgc                >  2:789475/1‑100 (MQ=255)

                                                                                 ccTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGcc               <  2:166127/100‑1 (MQ=255)

                                                                                 ccTGAATAACGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGcc               >  2:1810937/1‑100 (MQ=255)

                                                                                    gAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCAt             <  2:2494709/100‑1 (MQ=255)

                                                                                        aaCGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCcaa        >  2:2279695/1‑100 (MQ=255)

                                                                                         aCGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCcaac       <  2:274480/100‑1 (MQ=255)

                                                                                         aCGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCcaac       <  1:101273/100‑1 (MQ=255)

                                                                                         aCGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCcaac       <  1:789472/100‑1 (MQ=255)

                                                                                          cGGCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTCAAGAGGTTCTTTACGCTCAAGCGCCATGCcaaca      >  1:1357127/1‑100 (MQ=255)

                                                                                           ggCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGCTTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCcaacaa     <  2:1156804/100‑1 (MQ=255)

                                                                                           ggCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCcaacaa     <  1:1617740/100‑1 (MQ=255)

                                                                                           ggCACT‑GGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCcaacaa     >  2:665371/1‑100 (MQ=255)

                                                                                            gCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCcaacaa     <  2:1846636/100‑1 (MQ=255)

                                                                                             cACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCCAACAAt    >  1:1288085/1‑100 (MQ=255)

                                                                                               cTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCCAACAATac  >  2:1736524/1‑100 (MQ=255)

GATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCCAACAATAC  >  NC_000913/4503068‑4503262

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