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	97.7%	549.9 / 5.1	171	coding (91/834 nt)	mcrA	McrA
Reads supporting (aligned to +/- strand): ref base . (2/2); new base G (89/78); total (91/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 193 total aligned reads displayed.

TTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCG  >  NC_000913/1209567‑1209746

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

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

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

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

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

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

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

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

            ttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGt                                                                              <  2:2060836/93‑1 (MQ=255)

            ttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGt                                                                              <  2:1859696/93‑1 (MQ=255)

            ttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGt                                                                              <  2:1357163/93‑1 (MQ=255)

             tGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTg                                                                             <  1:1042036/93‑1 (MQ=255)

              gATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGGGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGa                                                                            >  1:2772339/1‑93 (MQ=255)

               ataataATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGAc                                                                           >  1:4157508/1‑93 (MQ=255)

                taataaTGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACa                                                                          >  1:1955900/1‑93 (MQ=255)

                taataaTGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACa                                                                          >  2:1487490/1‑93 (MQ=255)

                 aataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAg                                                                         >  2:1988169/1‑93 (MQ=255)

                  ataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGa                                                                        <  1:4136822/93‑1 (MQ=255)

                  ataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑GGG‑GGCCGGGTGACAGaa                                                                       >  1:6467491/1‑93 (MQ=255)

                   taatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTGTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGaa                                                                       >  1:2510660/1‑93 (MQ=255)

                   taatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGaa                                                                       <  2:3206554/93‑1 (MQ=255)

                      tGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACa                                                                    >  1:3767973/1‑93 (MQ=255)

                      tGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACa                                                                    <  1:2423634/93‑1 (MQ=255)

                      tGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACa                                                                    >  2:4862019/1‑93 (MQ=255)

                       ggAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACaa                                                                   >  2:6618547/1‑93 (MQ=255)

                          aTTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGa                                                                >  2:5148289/1‑93 (MQ=255)

                          aTTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGa                                                                >  1:4903185/1‑93 (MQ=255)

                            tGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGAtt                                                              <  2:822313/93‑1 (MQ=255)

                            tGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGAtt                                                              >  1:1632901/1‑93 (MQ=255)

                             gAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTa                                                             >  1:4234768/1‑93 (MQ=255)

                               aCTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACa                                                           >  2:6630871/1‑93 (MQ=255)

                               aCTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACa                                                           >  1:6074215/1‑93 (MQ=255)

                                cTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACaa                                                          <  1:3959239/93‑1 (MQ=255)

                                 tGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAat                                                         >  1:4348015/1‑93 (MQ=255)

                                   aaaGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAatat                                                       >  2:5339185/1‑93 (MQ=255)

                                    aaGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATc                                                      >  1:2595503/1‑93 (MQ=255)

                                    aaGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATc                                                      <  1:5390940/93‑1 (MQ=255)

                                      gCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGc                                                    >  2:5499649/1‑93 (MQ=255)

                                      gCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGc                                                    <  1:831734/93‑1 (MQ=255)

                                        tGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGctt                                                  >  2:5365781/1‑92 (MQ=255)

                                        tGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGctc                                                  >  1:1769738/1‑93 (MQ=255)

                                        tGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGctc                                                  <  2:2800131/93‑1 (MQ=255)

                                        tGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGctc                                                  <  1:4539595/93‑1 (MQ=255)

                                           gtgtTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTg                                               >  2:5781210/1‑93 (MQ=255)

                                            tgtTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGa                                              <  2:2983656/93‑1 (MQ=255)

                                             gtTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGAt                                             <  1:2907516/93‑1 (MQ=255)

                                             gtTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGAt                                             <  1:1746537/93‑1 (MQ=255)

                                              ttCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGAtt                                            <  1:734108/93‑1 (MQ=255)

                                              ttCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑GGGCGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGAtt                                            <  2:4693883/93‑1 (MQ=255)

                                               tCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATta                                           <  2:2197025/93‑1 (MQ=255)

                                                cGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATtat                                          >  2:5572330/1‑93 (MQ=255)

                                                 gATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATtata                                         <  2:5957122/93‑1 (MQ=255)

                                                 gATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATtata                                         >  2:188487/1‑93 (MQ=255)

                                                 gATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATtata                                         <  1:4562027/93‑1 (MQ=255)

                                                 gATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATtata                                         >  2:1834534/1‑93 (MQ=255)

                                                  aTAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTatat                                        <  2:2654629/93‑1 (MQ=255)

                                                    aGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCa                                      <  1:5620023/93‑1 (MQ=255)

                                                       tGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTg                                   <  2:3849340/93‑1 (MQ=255)

                                                        gAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGa                                  >  1:2172099/1‑93 (MQ=255)

                                                         aaGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGaa                                 >  1:1355710/1‑93 (MQ=255)

                                                           gagGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACg                               >  2:6094326/1‑93 (MQ=255)

                                                           gagGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACg                               <  2:3522276/93‑1 (MQ=255)

                                                             ggATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTCCAATATCGCTCTTGATTATATCATTGAACggt                             <  2:5068911/93‑1 (MQ=255)

                                                             ggATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggt                             >  2:4383975/1‑93 (MQ=255)

                                                             ggATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggt                             <  1:4192508/93‑1 (MQ=255)

                                                             ggATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggt                             >  1:6184112/1‑93 (MQ=255)

                                                             ggATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggt                             >  1:4108454/1‑93 (MQ=255)

                                                             ggATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggt                             >  1:2206335/1‑93 (MQ=255)

                                                             ggATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggt                             <  2:2636530/93‑1 (MQ=255)

                                                              gATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggtt                            <  2:6708350/93‑1 (MQ=255)

                                                              gATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggtt                            <  1:5729311/93‑1 (MQ=255)

                                                              gATGGTCTTGATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggtt                            <  2:2413305/93‑1 (MQ=255)

                                                               aTGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttg                           <  1:2436989/93‑1 (MQ=255)

                                                                tGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttgg                          >  2:4136781/1‑93 (MQ=255)

                                                                tGGTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttgg                          >  2:3346234/1‑93 (MQ=255)

                                                                tGGTGTTTATGGTCTAATCCTTGAGTCGT‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttggt                         >  1:786878/1‑93 (MQ=255)

                                                                 ggTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttggt                         <  2:4278664/93‑1 (MQ=255)

                                                                 ggTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttggt                         <  1:1453851/93‑1 (MQ=255)

                                                                 ggTGTTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttggt                         <  1:6646509/93‑1 (MQ=255)

                                                                 ggTGTTTAGGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttggt                         <  1:979043/93‑1 (MQ=255)

                                                                  gggtTTATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGCTTACAATATCGCTCTTGATTATATCATTGAACggttggtt                        >  1:5259288/3‑93 (MQ=255)

                                                                     tttATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGAt                     <  2:1664382/93‑1 (MQ=255)

                                                                      ttATGGTCTAATCCTTGAGTCGTGGGG‑GGTCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGAtt                    <  1:4781482/93‑1 (MQ=255)

                                                                      ttATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGAtt                    >  2:2776508/1‑93 (MQ=255)

                                                                       tATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTc                   <  2:6329387/93‑1 (MQ=255)

                                                                       tATGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTc                   <  2:3480700/93‑1 (MQ=255)

                                                                        aTGGTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCt                  >  1:4535161/1‑93 (MQ=255)

                                                                           gTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGgt               <  1:3166097/93‑1 (MQ=255)

                                                                           gTCTAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGgt               <  1:6651518/93‑1 (MQ=255)

                                                                              tAATCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTa            >  2:2276497/1‑93 (MQ=255)

                                                                               aaTCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTTAACGGTTGGTTGATTCTGGTGTAt           >  2:6419921/1‑93 (MQ=255)

                                                                               aaTCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTAt           <  1:2061285/93‑1 (MQ=255)

                                                                                aTCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATc          >  1:2256070/1‑93 (MQ=255)

                                                                                aTCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATc          <  1:154277/93‑1 (MQ=255)

                                                                                aTCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATc          >  2:4523786/1‑93 (MQ=255)

                                                                                 tCCTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATcc         >  2:4217141/1‑93 (MQ=255)

                                                                                  ccTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATccc        <  2:5398482/93‑1 (MQ=255)

                                                                                   cTTGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCa       >  2:5426859/1‑93 (MQ=255)

                                                                                    ctGAGTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGATGTATCCCaa      <  2:3669592/92‑1 (MQ=255)

                                                                                        gTCGTGGGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCg  >  2:3758406/1‑93 (MQ=255)

TTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑GGG‑GGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCG  >  NC_000913/1209567‑1209746

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