BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,209,664	(G)_5→6	coding (96/834 nt)	mcrA →	McrA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,209,659	1	.	G	96.9%	912.1 / 18.1	287	coding (91/834 nt)	mcrA	McrA
Reads supporting (aligned to +/- strand): ref base . (2/7); new base G (135/143); total (137/150)
Fisher's exact test for biased strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 351 total aligned reads displayed.

TTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAG  >  NC_000913/1209567‑1209749

ttatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGt                                                                                               <  2:2126129/93‑1 (MQ=255)

ttatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGt                                                                                               >  2:461719/1‑93 (MQ=255)

ttatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGt                                                                                               >  2:5394497/1‑93 (MQ=255)

 tatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGACGAGGATGGTGGTTATGGTCTAATCCTTGAGTAGT‑‑g                                                                                            >  2:236761/1‑93 (MQ=255)

 tatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑g                                                                                            <  2:3429252/93‑1 (MQ=255)

 tatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑g                                                                                            <  2:283490/93‑1 (MQ=255)

 tatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑g                                                                                            >  2:942017/1‑93 (MQ=255)

  atGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑gg                                                                                           <  1:4763590/93‑1 (MQ=255)

   tGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑ggg                                                                                          <  2:5297467/93‑1 (MQ=255)

   tGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑ggg                                                                                          >  2:5865066/1‑93 (MQ=255)

    gCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑ggg‑g                                                                                        <  1:790900/93‑1 (MQ=255)

    gCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑ggg‑g                                                                                        <  2:9576781/93‑1 (MQ=255)

     cATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGCGTCGT‑‑ggg‑gg                                                                                       >  1:4021270/1‑93 (MQ=255)

     cATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑ggg‑gg                                                                                       <  1:307588/93‑1 (MQ=255)

      aTGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑ggg‑ggg                                                                                      >  1:8511415/1‑92 (MQ=255)

       tGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑ggg‑gggg                                                                                     >  1:6622807/1‑91 (MQ=255)

       tGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑ggg‑gggg                                                                                     >  1:3550236/1‑91 (MQ=255)

       tGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑ggg‑gggc                                                                                     <  2:1372567/93‑3 (MQ=255)

       tGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑ggg‑gggc                                                                                     <  1:8047875/93‑3 (MQ=255)

       tGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑ggg‑gggc                                                                                     >  2:4882877/1‑91 (MQ=255)

        gTTTTGGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑ggg‑ggggc                                                                                    <  1:4577441/93‑4 (MQ=255)

         tttttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑ggg‑ggggcg                                                                                   >  1:4870318/1‑89 (MQ=255)

          ttttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑GGG‑GGCCggg                                                                                  <  1:7197129/93‑1 (MQ=255)

           tttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCggg                                                                                  <  2:6920640/93‑1 (MQ=255)

            ttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGt                                                                                 <  1:5433204/93‑1 (MQ=255)

                 aataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAg                                                                            >  2:4954649/1‑93 (MQ=255)

                 aataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAg                                                                            <  1:8418227/93‑1 (MQ=255)

                  ataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGCGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGa                                                                           >  1:1710667/1‑93 (MQ=255)

                  ataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGa                                                                           <  1:2791748/93‑1 (MQ=255)

                  ataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGa                                                                           <  2:4745034/93‑1 (MQ=255)

                   taatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGaa                                                                          >  2:983696/1‑93 (MQ=255)

                    aatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGGCGGGTGACAGaaa                                                                         <  2:734272/93‑1 (MQ=255)

                    aatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGaaa                                                                         >  2:5544340/1‑93 (MQ=255)

                     atGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAAc                                                                        >  1:4264148/1‑93 (MQ=255)

                     atGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAAc                                                                        <  2:923631/93‑1 (MQ=255)

                      tGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACa                                                                       >  1:4478053/1‑93 (MQ=255)

                       ggAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACaa                                                                      <  1:5708175/93‑1 (MQ=255)

                       ggAATAGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACaa                                                                      <  2:4474653/93‑1 (MQ=255)

                        gAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGGCGGGCGACAGAaaaaaa                                                                     >  2:7445530/1‑93 (MQ=255)

                        gAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACaaa                                                                     <  2:8508334/93‑1 (MQ=255)

                         aaTTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGGCGGGTGACAGAAACAAAg                                                                    <  2:5908641/93‑1 (MQ=255)

                         aaTTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAg                                                                    <  2:7083414/93‑1 (MQ=255)

                          aTTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGa                                                                   >  1:1926549/1‑93 (MQ=255)

                          aTTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑GGGTGGCCGGGTGACAGAAACAAAGa                                                                   <  2:2771767/93‑1 (MQ=255)

                             gAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTa                                                                <  1:5802295/93‑1 (MQ=255)

                                cTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGGCGGGTGACAGAAACAAAGATTACaa                                                             <  1:1445509/93‑1 (MQ=255)

                                cTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGGCGGGTGACAGAAACAAAGATAACaa                                                             >  1:2616129/1‑93 (MQ=255)

                                cTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACaa                                                             <  1:6064271/93‑1 (MQ=255)

                                 tGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAat                                                            <  1:2214671/93‑1 (MQ=255)

                                 cgAAAGCTGAGTGGTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAat                                                            <  1:7700772/92‑1 (MQ=255)

                                  gAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAata                                                           >  1:1215794/1‑93 (MQ=255)

                                   aaaGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAatat                                                          <  2:8336452/93‑1 (MQ=255)

                                     aGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCg                                                        >  1:6818195/1‑93 (MQ=255)

                                        tGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGctc                                                     >  2:8100542/1‑93 (MQ=255)

                                         gAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGGCGGGTGACAGAAACAAAGATTACAATATCGctct                                                    <  1:561873/93‑1 (MQ=255)

                                          aGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGGCGGGTGACAGAAACAAAGATTACAATATCGCTCtt                                                   <  2:6501764/93‑1 (MQ=255)

                                          aGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCtt                                                   <  2:7869352/93‑1 (MQ=255)

                                           gtgtTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTg                                                  <  2:1592817/93‑1 (MQ=255)

                                           gtgtTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTg                                                  >  1:751193/1‑93 (MQ=255)

                                            tgtTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGa                                                 >  2:1104446/1‑93 (MQ=255)

                                             gtTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGAtt                                               <  1:9976542/93‑1 (MQ=255)

                                              ttCGATAGGTTAAGAGGATTGTGTTTATGGTCTAATCCTTTAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGAtt                                               <  1:2548131/93‑1 (MQ=255)

                                              ttCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGAtt                                               >  2:8717333/1‑93 (MQ=255)

                                              ttCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGAtt                                               <  2:1042692/93‑1 (MQ=255)

                                               tCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGGCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATta                                              >  1:5886005/1‑93 (MQ=255)

                                               tCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATta                                              <  2:354513/93‑1 (MQ=255)

                                               tCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATta                                              >  1:1973548/1‑93 (MQ=255)

                                               tCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATta                                              <  1:9735591/93‑1 (MQ=255)

                                                cGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGGCGTG‑GGG‑GGGCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATtat                                             <  1:1155980/93‑1 (MQ=255)

                                                 gATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGGCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATtata                                            <  1:637572/93‑1 (MQ=255)

                                                 gATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATtata                                            <  2:9376262/93‑1 (MQ=255)

                                                  aTAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGGCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTatat                                           >  1:1022235/1‑93 (MQ=255)

                                                   tAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATc                                          <  2:1946012/93‑1 (MQ=255)

                                                    aGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGGCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCa                                         <  1:3361599/93‑1 (MQ=255)

                                                     ggTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑‑GG‑GGCCGCGGGACAGAAACAAAGAGTATAATATCTCTCTTTATTTTATCtttt                                      >  1:8200383/1‑92 (MQ=255)

                                                      gTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCAtt                                       >  2:886120/1‑93 (MQ=255)

                                                        gAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGa                                     >  2:3149628/1‑93 (MQ=255)

                                                         aaGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGaa                                    <  1:8627169/93‑1 (MQ=255)

                                                         aaGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGaa                                    <  1:5299582/93‑1 (MQ=255)

                                                           gagGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACg                                  >  2:7489793/1‑93 (MQ=255)

                                                            agGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACgg                                 <  2:2538277/93‑1 (MQ=255)

                                                             ggATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggt                                <  2:7246025/93‑1 (MQ=255)

                                                              gATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggtt                               >  2:606137/1‑93 (MQ=255)

                                                              gATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggtt                               <  2:5649485/93‑1 (MQ=255)

                                                               aTGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttg                              >  1:3938987/1‑93 (MQ=255)

                                                               aTGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttg                              >  1:9088572/1‑93 (MQ=255)

                                                               aTGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttg                              <  2:6355987/93‑1 (MQ=255)

                                                                tGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTTCAATATCGCTCTTGATTATATCATTGAACggttgg                             >  2:8800080/1‑93 (MQ=255)

                                                                 ggTGTTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttggt                            <  2:9027199/93‑1 (MQ=255)

                                                                    gtTTATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGa                         >  2:4001942/1‑93 (MQ=255)

                                                                     tttATGGTCTAATCCTTGAGTCGTG‑GGG‑GGGCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGAt                        >  1:9609560/1‑93 (MQ=255)

                                                                      ttATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGAtt                       >  1:3285897/1‑93 (MQ=255)

                                                                       tATGGTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTc                      <  1:2319269/93‑1 (MQ=255)

                                                                         tGGTCTAATCCTTGAGTCGTG‑GGG‑GGGCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTg                    <  2:3083378/93‑1 (MQ=255)

                                                                          ggTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTgg                   >  2:6187184/1‑93 (MQ=255)

                                                                           gTCTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGgt                  <  2:4310408/93‑1 (MQ=255)

                                                                             cTAATCCTTGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGgtgt                >  1:775350/1‑93 (MQ=255)

                                                                              tAATCCTTGAGTCGTG‑GGG‑GGGCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTa               >  1:8850782/1‑93 (MQ=255)

                                                                                aTCCTTGAGTCGTGGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTAt              >  1:1300114/1‑93 (MQ=255)

                                                                                    ttGAGTCGTG‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCaa         <  2:4418852/93‑1 (MQ=255)

TTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAG  >  NC_000913/1209567‑1209749

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