BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,064,161	G→A	15.4%	L42L (CTG→CTA) ‡	yccE →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,064,161	0	G	A	15.4%	385.2 / 45.1	159	L42L (CTG→CTA)	yccE	uncharacterized protein
Reads supporting (aligned to +/- strand): ref base G (57/76); new base A (17/8); total (75/84)
Fisher's exact test for biased strand distribution p-value = 2.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 168 total aligned reads displayed.

TTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAA  >  NC_000913/1064062‑1064260

ttGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTg                                                                                                     >  2:3067904/1‑100 (MQ=255)

 tGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTgg                                                                                                    <  2:1865854/100‑1 (MQ=255)

 tGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTgg                                                                                                    >  2:4167755/1‑100 (MQ=255)

  gCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGgt                                                                                                   >  2:1169617/1‑100 (MQ=255)

    aCTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGgtg                                                                                                  <  2:931593/99‑1 (MQ=255)

     cTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGgtgtg                                                                                                <  1:1709507/100‑1 (MQ=255)

     cTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAgtgtg                                                                                                <  2:1988901/100‑1 (MQ=255)

       gCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGtt                                                                                              >  2:2639386/1‑100 (MQ=255)

          aataatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTa                                                                                           <  2:2154974/100‑1 (MQ=255)

          aataatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTa                                                                                           <  2:2121936/100‑1 (MQ=255)

           ataatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTAt                                                                                          <  2:747955/100‑1 (MQ=255)

           ataatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTAt                                                                                          <  1:573680/100‑1 (MQ=255)

           ataatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAGTGTGttttca                                                                                          <  2:3633363/100‑3 (MQ=255)

           ataatTATTTGAAACAGGCCTGGAACGATATAAAAAATAAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTAt                                                                                          <  1:3812291/100‑1 (MQ=255)

            taatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATg                                                                                         >  2:4052387/1‑100 (MQ=255)

             aatTATTTGAAACAGGCCTAAAAAGAAAAAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGc                                                                                        <  1:3967240/100‑1 (MQ=255)

              attattTGAAACAGGCCTGGAACGATATAAAAAATGGGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCg                                                                                       <  2:3123565/100‑1 (MQ=255)

              attattGGAAACAAGCCTGTAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCACTCACGCTTTTTGAAAACACACCAGTGTGttttca                                                                                          <  1:2470561/97‑3 (MQ=255)

              agtattGGAGACAGGACTGGGAAGATATAAGAAATGTGGACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAGTGTGttttcaa                                                                                         <  1:3883408/96‑4 (MQ=255)

               ttattTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCgg                                                                                      <  2:3831106/100‑1 (MQ=255)

               ttattTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCgg                                                                                      <  2:2800406/100‑1 (MQ=255)

               ttattTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCgg                                                                                      <  2:3885724/100‑1 (MQ=255)

                tattTGAAACAGGCCTGGAAGGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGt                                                                                     <  1:1609205/100‑1 (MQ=255)

                tattTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGt                                                                                     >  2:3229623/1‑100 (MQ=255)

                tattTGAAACAGGCCTCGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAGTGTGttttcaaaaag                                                                                     <  1:2585228/100‑8 (MQ=38)

                 attTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGtt                                                                                    <  1:4094066/100‑1 (MQ=255)

                 attTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGtt                                                                                    >  2:56212/1‑100 (MQ=255)

                 attTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGt                                                                                     <  1:839681/99‑1 (MQ=255)

                  tttGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTta                                                                                   <  1:4138377/100‑1 (MQ=255)

                  tttGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTta                                                                                   >  2:1342968/1‑100 (MQ=255)

                         caagCCTTGAACAATATAAAAAATGAGTACAAAAAATATTAAATGGTTGCGATCTGTCTTTTTTAAAACACATTGTTGTGTTTTATGCGGTTAGACAATg                                                                            <  1:1666993/97‑1 (MQ=255)

                          aGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGa                                                                           >  1:4358024/1‑100 (MQ=255)

                           ggCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGaa                                                                          >  2:280841/1‑100 (MQ=255)

                           ggCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACAGGTGTGTTTTATGCGGTTATACAATGaa                                                                          >  1:2083948/1‑100 (MQ=255)

                            gCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAAc                                                                         <  2:2686559/100‑1 (MQ=255)

                            gCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAAc                                                                         >  1:3201790/1‑100 (MQ=255)

                            gCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAAc                                                                         <  2:3165363/100‑1 (MQ=255)

                               tGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCa                                                                      <  2:4137070/100‑1 (MQ=255)

                                ggAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAg                                                                     >  2:849142/1‑100 (MQ=255)

                                 gAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAg                                                                     <  1:1982944/99‑1 (MQ=255)

                                 gAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGa                                                                    >  2:2352583/1‑100 (MQ=255)

                                  aaCGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGAc                                                                   <  2:609143/100‑1 (MQ=255)

                                  aaCGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGAc                                                                   <  2:1749205/100‑1 (MQ=255)

                                   aCGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACg                                                                  <  2:1659461/100‑1 (MQ=255)

                                   aCGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACg                                                                  <  1:1392360/100‑1 (MQ=255)

                                     gATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgta                                                                <  1:4414172/100‑1 (MQ=255)

                                     gATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgta                                                                <  1:2713808/100‑1 (MQ=255)

                                        ataAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAATACACTGATGTGTATTATGCGGTTATACAATGAACTAAGACgtaaag                                                             >  2:404981/1‑100 (MQ=255)

                                         taAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaaagt                                                            <  2:4298404/100‑1 (MQ=255)

                                           aaaaaTGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTAAAAACACACTGGTGTGTTTTATGCGGTTATACAATGACCTCAGACgtcacgta                                                           <  1:3514972/99‑1 (MQ=255)

                                              aaTGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATg                                                       <  2:351508/100‑1 (MQ=255)

                                                tGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGaa                                                     <  2:596627/100‑1 (MQ=255)

                                                  aGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAaga                                                   <  2:4298616/100‑1 (MQ=255)

                                                  aGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAaga                                                   <  2:2937245/100‑1 (MQ=255)

                                                  aGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAaga                                                   <  2:2039345/100‑1 (MQ=255)

                                                   gTACGAAAAAAATCAAACATATTCAATCACGCTTTTTTTTTTTTTTTTTTTGTGTTTTATGCGGTTATACAATGAACTTAGACGTAAAGTAAATGAagag                                                  >  1:2302976/1‑100 (MQ=255)

                                                   gTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAagag                                                  <  1:2918986/100‑1 (MQ=255)

                                                      cGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACAGTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGGAAATGAAGAGGAt                                               >  1:363746/1‑100 (MQ=255)

                                                      cGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGAt                                               >  2:1044400/1‑100 (MQ=255)

                                                      cGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGAt                                               >  2:335826/1‑100 (MQ=255)

                                                         aaaaaaTCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACt                                            >  2:3260583/1‑100 (MQ=255)

                                                          aaaaaTCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTc                                           <  2:860874/100‑1 (MQ=255)

                                                              aTCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATg                                       >  2:3016925/1‑100 (MQ=255)

                                                                 aaaCATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTcagatcgg                                    >  1:1858361/1‑93 (MQ=255)

                                                                  aaCATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTg                                   >  1:756147/1‑100 (MQ=255)

                                                                      tatTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAAt                               >  1:3484531/1‑100 (MQ=255)

                                                                        ttCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAAtgt                             >  2:3336678/1‑100 (MQ=255)

                                                                        ttCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAAtgt                             >  1:10405/1‑100 (MQ=255)

                                                                         tCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAAtgtg                            <  1:2345042/100‑1 (MQ=255)

                                                                           aaTCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGaa                          <  1:24392/100‑1 (MQ=255)

                                                                              cACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCa                       >  2:2303521/1‑100 (MQ=255)

                                                                              cACGCCTTTTGAAAGCACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCa                       >  2:2494982/1‑100 (MQ=255)

                                                                               ataaccgcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCAc                      >  2:3342356/11‑100 (MQ=38)

                                                                                cGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACt                     >  2:18361/1‑100 (MQ=255)

                                                                                cGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACt                     <  2:3728349/100‑1 (MQ=255)

                                                                                 gCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTa                    <  1:131038/100‑1 (MQ=255)

                                                                                 gCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTa                    >  1:3848967/1‑100 (MQ=255)

                                                                                 aaccgcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGCCGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTa                    >  2:1488569/9‑100 (MQ=38)

                                                                                  cTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAg                   <  2:4414264/100‑1 (MQ=255)

                                                                                  cTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAg                   >  2:2615918/1‑100 (MQ=255)

                                                                                  cTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAg                   >  2:1343850/1‑100 (MQ=255)

                                                                                  accgcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTAAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAg                   >  1:4031838/8‑100 (MQ=38)

                                                                                   tttttGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGa                  >  1:1530433/1‑100 (MQ=255)

                                                                                   tttttGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGa                  <  2:363747/100‑1 (MQ=255)

                                                                                   tttttGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGa                  <  1:4010827/100‑1 (MQ=255)

                                                                                    cgcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaa                 >  2:1727706/6‑100 (MQ=255)

                                                                                    cgcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaa                 >  2:2217291/6‑100 (MQ=255)

                                                                                     gcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaa                >  1:992501/5‑100 (MQ=255)

                                                                                     gcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaa                >  1:1238446/5‑100 (MQ=255)

                                                                                      ttGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaaa               <  1:1343845/100‑1 (MQ=255)

                                                                                      cataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaaa               >  1:697824/4‑100 (MQ=255)

                                                                                       tGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaaaa              >  1:1188622/1‑100 (MQ=255)

                                                                                        gAAAAAACACTAGTGTGTTTTATGCGGTTATACAAAGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAg                   <  1:3677733/94‑1 (MQ=255)

                                                                                            acacacGGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACCCCATGTCTGGAATGTGAATCACTAGAAAAAGAAt          <  1:2374386/99‑1 (MQ=255)

                                                                                              acacTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTg       <  2:2649459/100‑1 (MQ=255)

                                                                                                acTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAg     <  1:1748515/100‑1 (MQ=255)

                                                                                                 cTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAgg    >  2:1822954/1‑100 (MQ=255)

                                                                                                 cTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAgg    >  2:1634583/1‑100 (MQ=255)

                                                                                                 cTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAg     >  1:3409941/1‑99 (MQ=255)

                                                                                                   ggTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGaa  >  1:2474564/1‑100 (MQ=255)

TTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAA  >  NC_000913/1064062‑1064260

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