BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	625,985	G→A	6.2%	A367A (GCG→GCA)	entC →	isochorismate synthase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	625,985	0	G	A	6.2%	300.9 / 6.6	96	A367A (GCG→GCA)	entC	isochorismate synthase 1
Reads supporting (aligned to +/- strand): ref base G (42/48); new base A (3/3); total (45/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01

Only 100 of 103 total aligned reads displayed.

TGGAACCGTTCGACCGCGAACTGTTTGGCGGCATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGCATTCCATTCACCCGCTGGCCGGAAGAGTTTGCCCGTCGCTATCGGGAAAAAGGCTACT  >  NC_000913/625840‑626133

tGGAACCGTTCGACCGCGAACTGTTTGGCGGCATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGcgtcg                                                                                                                                                   >  5:343477/1‑149 (MQ=255)

tGGAACCGTTCGACCGCGAACTGTTTGGCGGCATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGcgtcg                                                                                                                                                   >  7:29869/1‑149 (MQ=255)

 ggAACCGTTCGACCGCGAACTGTTTGGCGGCATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGcgtcgt                                                                                                                                                  >  6:414304/1‑149 (MQ=255)

  gAACCGTTCGACCGCGAACTGTTTGGCGGCATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGcgtcgt                                                                                                                                                  >  2:481205/1‑148 (MQ=255)

  gAACCGTTCGACCGCGAACTGTTTGGCGGCATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGcgtcgt                                                                                                                                                  >  6:1575/1‑148 (MQ=255)

     ccGTTCGACCGCGCACTGTTTGGCGGCATTGTGGGGTGGTGTGACAGCGAAGGTAACGGCGAAGGGGTGGGGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCAcc                                                                                                                                              <  5:242137/149‑1 (MQ=255)

     ccGTTCGACCGCGAACTGTTTGGCGGCATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTTGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCTGGGATTGTGCCTGCGTCGTCAcc                                                                                                                                              <  2:341250/149‑1 (MQ=255)

      cGTTCGACCGCGAACTGTTTGGCGGCATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCg                                                                                                                                             >  4:218417/1‑149 (MQ=255)

        ttCGACCGCGAACTGTTTGGCGGCATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGtt                                                                                                                                           >  2:465990/1‑149 (MQ=255)

           gACCGCGAACTGTTTGGCGGCATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTggg                                                                                                                                        <  4:365759/149‑1 (MQ=255)

           gACCGCGAACTGTTTGGCGGCATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTggg                                                                                                                                        >  6:264655/1‑149 (MQ=255)

            aCCGCGCACTGTTTGGCGGCATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGCCCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGt                                                                                                                                       >  4:10257/1‑149 (MQ=255)

             ccGCGAACTGTTTGGCGGCATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGt                                                                                                                                       >  7:114125/1‑148 (MQ=255)

              cgcgAACTGTTTGGCGGCATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGt                                                                                                                                       >  1:404979/1‑147 (MQ=255)

               gcgAACTGTTTGGCGGCATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAg                                                                                                                                    <  4:141400/149‑1 (MQ=255)

               gcgAACTGTTTGGCGGCATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAg                                                                                                                                    >  2:454131/1‑149 (MQ=255)

                  aaCTGTTTGGCGGCATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTg                                                                                                                                          >  2:287829/1‑140 (MQ=255)

                  aaCTGTTTGGCGGCATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTg                                                                                                                                          <  1:287829/140‑1 (MQ=255)

                  aaCTGTTTGGCGGCATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTgg                                                                                                                                 >  4:71561/1‑149 (MQ=255)

                       tttGGCGGCATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGgcgc                                                                                                                              <  4:1830/147‑1 (MQ=255)

                          ggcggcATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGaaa                                                                                                                          <  1:465990/148‑1 (MQ=255)

                             ggcATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTgg                                                                                                                                         <  5:347280/130‑1 (MQ=255)

                             ggcATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTgg                                                                                                                                         >  6:347280/1‑130 (MQ=255)

                                 ttGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGt                                                                                                                   >  6:267662/1‑148 (MQ=255)

                                 ttGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCATCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGt                                                                                                                   >  1:334409/1‑148 (MQ=255)

                                 ttGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCATCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGt                                                                                                                   >  3:209071/1‑148 (MQ=255)

                                 ttGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCATCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGt                                                                                                                   >  3:202059/1‑148 (MQ=255)

                                    tgGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCaaa                                                                                                               >  6:316419/1‑149 (MQ=255)

                                      ggTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACt                                                                                                             <  3:324988/149‑1 (MQ=255)

                                           gtgtGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTa                                                                                                        <  1:454131/149‑1 (MQ=255)

                                             gtgACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTAcc                                                                                                      <  5:316419/149‑1 (MQ=255)

                                              tgACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGtt                                                                                                                                           <  7:389747/111‑1 (MQ=255)

                                              tgACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGtt                                                                                                                                           >  8:389747/1‑111 (MQ=255)

                                              tgACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTAcc                                                                                                      <  5:414304/148‑1 (MQ=255)

                                                aCAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCCACCATg                                                                                                   <  1:479103/149‑1 (MQ=255)

                                                 cAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGt                                                                                                  <  3:218417/149‑1 (MQ=255)

                                                  aGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGc                                                                                                                     <  2:408491/129‑1 (MQ=255)

                                                  aGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGc                                                                                                                     >  1:408491/1‑129 (MQ=255)

                                                  aGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGGTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACAATGtt                                                                                                 >  1:61557/1‑149 (MQ=255)

                                                  aGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGCAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGtt                                                                                                 <  8:132158/149‑1 (MQ=255)

                                                   gCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGtt                                                                                                 >  5:113940/1‑148 (MQ=255)

                                                    cGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGa                                                                                               <  5:267662/149‑1 (MQ=255)

                                                     gAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGaa                                                                                              >  2:83237/1‑149 (MQ=255)

                                                        ggTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGACGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGCACGt                                                                                           >  3:119807/1‑149 (MQ=255)

                                                          gaaCGGCGAAGGGGTGGTGACCATCCGCTGCGCGCAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGAGGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGttt                                                                                         <  6:237048/148‑1 (MQ=255)

                                                            aCGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGttttt                                                                                       <  5:264655/149‑1 (MQ=255)

                                                            aCGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGttttt                                                                                       <  5:485545/149‑1 (MQ=255)

                                                              ggCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTg                                                                                      >  1:461726/1‑148 (MQ=255)

                                                               gCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGa                                                                                    <  1:481205/149‑1 (MQ=255)

                                                               gCGAATGGGTGGTGAACATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGGCTGCGTCGTCACCGTTGTGTGAGTGGCGCGAAACAGGCGGCAAACTTTCTACCATGTTGAACGTTTTTGGa                                                                                    >  1:297287/1‑149 (MQ=255)

                                                                    tGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCa                                                                               <  7:476404/149‑1 (MQ=255)

                                                                     gggTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGGCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCAt                                                                              >  3:119519/1‑149 (MQ=255)

                                                                     gggTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCAt                                                                              >  2:33671/1‑149 (MQ=255)

                                                                      ggtggtGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTg                                                                                                                                      <  4:116704/92‑1 (MQ=255)

                                                                      ggtggtGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTg                                                                                                                                      >  3:116704/1‑92 (MQ=255)

                                                                        tggtgACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTaa                                                                           <  4:265496/149‑1 (MQ=255)

                                                                         ggtgACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCATCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAg                                                                          <  4:209071/149‑1 (MQ=255)

                                                                          gtgACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGgcgc                                                                                                                              >  7:459808/1‑96 (MQ=255)

                                                                          gtgACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCATCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAg                                                                          <  4:202059/148‑1 (MQ=255)

                                                                          gtgACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCATCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAg                                                                          <  2:334409/148‑1 (MQ=255)

                                                                            gACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAATCAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCAt                                                                              <  3:315113/142‑1 (MQ=255)

                                                                            gACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAg                                                                       >  4:334156/1‑149 (MQ=255)

                                                                            gACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAg                                                                       >  1:460123/1‑149 (MQ=255)

                                                                            gACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTAGAACGTTTTTGGATAGCAt                                                                              >  4:315113/1‑142 (MQ=255)

                                                                             aCCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGc                                                                      >  3:35687/1‑149 (MQ=255)

                                                                             aCCATCCGCTGCGCGAAGCTGCGCGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACACGCGTCAAAATTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGc                                                                      >  6:397419/1‑149 (MQ=255)

                                                                               cATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCCCCGTTGGGTGAGTGGCGCGAAACAGGCGTCACACTTTCACCCATGGTGAACGGTTTGGGAGTGCAGTAAGGAGCgc                                                                    >  2:231091/1‑148 (MQ=255)

                                                                               cATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGc                                                                      >  1:344169/1‑147 (MQ=255)

                                                                               cATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGGCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGa                                                                    <  4:35687/149‑1 (MQ=255)

                                                                                   cGCTGCGCGAAGCTGCGGGAAAATCAGGTGGGTCTGAGTGCCTGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGGGAGTGGCGCGAAACAGGCGGCAAACATTCTACCATGTTGAACGGTTTTGGATTGCATTAAGGAGCGAGGAt                                                                <  4:119807/149‑1 (MQ=255)

                                                                                   cGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGGCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAAGGTTTTTGGATTGCATTAAGGAGCGAGGAt                                                                <  5:397419/149‑1 (MQ=255)

                                                                                     cTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGa                                                              <  2:461726/149‑1 (MQ=255)

                                                                                         gcgAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAgg                                                                                                                      <  6:458524/89‑1 (MQ=255)

                                                                                         gcgAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAgg                                                                                                                      >  5:458524/1‑89 (MQ=255)

                                                                                         gcgAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGcat                                                          >  7:115012/1‑149 (MQ=255)

                                                                                           gAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGcattc                                                        >  8:281091/1‑149 (MQ=255)

                                                                                             aGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGcattcca                                                      <  3:71561/149‑1 (MQ=255)

                                                                                              gCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGcattccat                                                     <  2:460123/149‑1 (MQ=255)

                                                                                               cTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGcattccatt                                                    <  1:83237/149‑1 (MQ=255)

                                                                                                   gggAAAATCAGGTGAGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGCATTCCATTCAcc                                                >  8:124254/1‑149 (MQ=255)

                                                                                                          tCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGAt                                                                <  4:155555/126‑1 (MQ=255)

                                                                                                          tCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGAt                                                                >  3:155555/1‑126 (MQ=255)

                                                                                                                 cggcTGATTGCCGGAGCGGGGACTGTGCCTGCGTCGTCACCGTTGTTTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAATTAGCGAGGATGAGCATTCCATTCACCCGCTGGCCGGAaga                                  <  4:119519/146‑1 (MQ=255)

                                                                                                                 cGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGCATTCCATTCACCCGCTGGCCGGAaga                                  <  5:1575/149‑1 (MQ=255)

                                                                                                                   tCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGCATTCCATTCACCCGc                                             >  1:228988/1‑136 (MQ=255)

                                                                                                                   tCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGCATTCCATTCACCCGc                                             <  2:228988/136‑1 (MQ=255)

                                                                                                                       tttGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGCATTCCATTCACCCGCTGGCCGGAAGAGTTTGc                            <  6:113940/149‑1 (MQ=255)

                                                                                                                       tttGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGCATTCCATTCACCCGCTGGCCGGAAGAGTTTGc                            <  1:33671/149‑1 (MQ=255)

                                                                                                                       tttGCCGGAGCGGGGATTGTGCCTGAGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGCATTCCATTCACCCGCTGGCCGGAAGAGttt                              <  6:343477/147‑1 (MQ=255)

                                                                                                                           caggAGCGTGGATTGTGCCTGCGTAGGCACCGTTGGGGGAGGGGCGCGAAACAGGCGTCAAACTTTCTACCATTGTGAACGTTTTTGGAGTGCATTAAGTAGCGAGGATGAGCATTCCATTCACCCGCTGGACTTATGATTTTGCCCGt                        <  2:297287/147‑1 (MQ=255)

                                                                                                                             ggAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGCATTCCATTCACCCGCTGGCCGGAAGAGTTTGCCCGTCg                      >  5:188222/1‑149 (MQ=255)

                                                                                                                                  ggggATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGCATTCCATTCACCCGCTGGCCGGAAGAGTTTGCCCGTCGCTATc                 >  3:249745/1‑149 (MQ=255)

                                                                                                                                   gggAGTGTGCCTGCGTCGTCACCGTGGGGTGAGTGGCGCGAAACAGGCGTCAAACTTGCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGCATTCCATTCACCCGCTGGCCGGAAGAGTTTGCCCGTCGCTATCg                <  2:395631/149‑1 (MQ=255)

                                                                                                                                    ggATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGCATTCCATTCACCCGCTGGCCGGAAGAGTTTGCCCGTCGCTATCgg               >  5:312918/1‑149 (MQ=255)

                                                                                                                                          tgCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGATGATGAGCATTCCATTCACCCGCTGGCCGGAAGAGTTTGCCCGTCGCTATCGGGaaaaa         <  7:124254/149‑1 (MQ=255)

                                                                                                                                           gCCTGCGTCGTCACCGTTGGGTGAGTGTCGCGTAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGAGTGCATTAAGGAGCGAGGATGAGCATTCCATGCCCCCGCTGGCCGGAAGAGTTTGACCGTCGCTATCGGGAAAAAg        <  2:61557/149‑1 (MQ=255)

                                                                                                                                           gCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGCATTCCATTCACCCGCTGGCCGGAAGAGTTTGCCCGTCGCTATCGGGAAAAAg        >  1:499048/1‑149 (MQ=255)

                                                                                                                                             cTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGCATTCCATTCACCCGCTGGCCGGAAGAGTTTGCCCGTCGCTATCGGGAAAAAGGc      >  7:316621/1‑149 (MQ=255)

                                                                                                                                                cgtcgtCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGCATTCCATTCACCCGCTGGCCGGAAGAGTTTGCCCGTCGCTATCGGGAAAAAGGCTAc   <  2:344169/149‑1 (MQ=255)

                                                                                                                                                 gtcgtcACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGCATTCCATTCACCCACTGGCCGGAAGAGTTTGCCCGTCGCTATCGGGAAAAAGGCTACt  >  3:180022/1‑149 (MQ=255)

TGGAACCGTTCGACCGCGAACTGTTTGGCGGCATTGTGGGTTGGTGTGACAGCGAAGGTAACGGCGAATGGGTGGTGACCATCCGCTGCGCGAAGCTGCGGGAAAATCAGGTGCGTCTGTTTGCCGGAGCGGGGATTGTGCCTGCGTCGTCACCGTTGGGTGAGTGGCGCGAAACAGGCGTCAAACTTTCTACCATGTTGAACGTTTTTGGATTGCATTAAGGAGCGAGGATGAGCATTCCATTCACCCGCTGGCCGGAAGAGTTTGCCCGTCGCTATCGGGAAAAAGGCTACT  >  NC_000913/625840‑626133

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