BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,000,785	G→A	5.7%	G113G (GGC→GGT)	ycaC‑IV ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,000,785	0	G	A	5.7%	1101.6 / 33.6	348	G113G (GGC→GGT)	ycaC‑IV	putative hydrolase
Reads supporting (aligned to +/- strand): ref base G (162/165); new base A (9/11); total (171/177)
Fisher's exact test for biased strand distribution p-value = 8.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 351 total aligned reads displayed.

AGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTGGGGCC  >  NC_002947/6000636‑6000934

        aCCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCgat                                                                                                                                               >  1:6004836/1‑150 (MQ=255)

        aCCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCgat                                                                                                                                               >  2:2643281/1‑150 (MQ=255)

        aCCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCgat                                                                                                                                               >  1:5536943/1‑150 (MQ=255)

          ctgcctgGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCgatga                                                                                                                                             <  2:7212839/150‑1 (MQ=255)

             cctgGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCAACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCa                                                                                                                                          <  2:3361600/150‑1 (MQ=255)

              ctgGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAg                                                                                                                                         <  1:5891228/150‑1 (MQ=255)

                   tCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGct                                                                                                                                    >  2:8025814/1‑150 (MQ=255)

                     aTGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGcttc                                                                                                                                  >  1:4731175/1‑150 (MQ=255)

                           tCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCa                                                                                                                            >  1:8975899/1‑150 (MQ=255)

                              tGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTg                                                                                                                         <  2:2152548/150‑1 (MQ=255)

                              ggggCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGGTGCTTCTTGCCAGTg                                                                                                                         <  2:4155297/149‑1 (MQ=255)

                                      gTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCAc                                                                                                                 >  2:3155426/1‑150 (MQ=255)

                                      gTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCAc                                                                                                                 >  1:7199985/1‑150 (MQ=255)

                                       tCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACt                                                                                                                <  1:2504142/150‑1 (MQ=255)

                                            ggTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCtt                                                                                                           <  2:847225/150‑1 (MQ=255)

                                               cATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGAc                                                                                                        <  1:8406592/150‑1 (MQ=255)

                                                 tGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGa                                                                                                      >  1:5004311/1‑150 (MQ=255)

                                                  ttcattGAAGGTTCCCGACGCATCTGTGACGACAAAGACATCGAAACCGGCTTCGATAGCGGACAATGCCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGaa                                                                                                     >  2:4533426/5‑150 (MQ=255)

                                                  ttcattGAAGGTTCCCGACGCATCTGTGACGACAAAGACATCGAAACCGGCTTCGATAGCGGACAATGCCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGaa                                                                                                     >  1:1030459/5‑150 (MQ=255)

                                                  ttcattGAAGGTTCCCGACGCATCTGTGACGACAAAGACATCGAAACCGGCTTCGATAGCGGACAATGCCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGaa                                                                                                     >  2:642452/5‑150 (MQ=255)

                                                  ggCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGaa                                                                                                     <  2:1511740/150‑1 (MQ=255)

                                                  ggCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGaa                                                                                                     <  1:2050050/150‑1 (MQ=255)

                                                  ggCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGaa                                                                                                     <  1:3679610/150‑1 (MQ=255)

                                                     gTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTc                                                                                                  <  2:3477283/150‑1 (MQ=255)

                                                     attGAAGGTTCCCGACGCATCTGTGACGACAAAGACATCGAAACCGGCTTCGATAGCGGACAATGCCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTc                                                                                                  <  2:5675025/149‑1 (MQ=255)

                                                     attGAAGGTTCCCGACGCATCTGTGACGACAAAGACATCGAAACCGGCTTCGATAGCGGACAATGCCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTc                                                                                                  >  1:2147562/2‑150 (MQ=255)

                                                     attGAAGGTTCCCGACGCATCTGTGACGACAAAGACATCGAAACCGGCTTCGATAGCGGACAATGCCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTc                                                                                                  <  2:6315253/149‑1 (MQ=255)

                                                        gAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTc                                                                                               <  2:686740/150‑1 (MQ=255)

                                                        gAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTc                                                                                               <  2:9110949/150‑1 (MQ=255)

                                                        gAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTc                                                                                               <  1:3212456/150‑1 (MQ=255)

                                                           ggTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGtt                                                                                            <  2:24068/150‑1 (MQ=255)

                                                              aCCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTc                                                                                         <  2:2879724/150‑1 (MQ=255)

                                                               ccGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTcc                                                                                        >  2:8298013/1‑150 (MQ=255)

                                                               ccGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTcc                                                                                        >  2:8955348/1‑150 (MQ=255)

                                                                   aGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAgg                                                                                    >  2:4946160/1‑150 (MQ=255)

                                                                    ggCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGc                                                                                   <  2:407729/150‑1 (MQ=255)

                                                                    cgCATCTGTGACGACAAAGACATCGAAACCGGCTTCGATAGCGGACAATGCCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGc                                                                                   <  1:6730346/149‑1 (MQ=255)

                                                                       aTCAGTGACCACGAACACGTCGAACTCTTGTTCGAGGGCCGAGAGGGACTTTAAAGCTACGCTCCCTTGCGTCACAACTCCGGCGATGATGAGGTGGGTCTTGCGACGTGACTTCTCTTCCTTGACGTAGTCTTCGTTGTCCCAGGCGtt                                                                                <  2:7644962/150‑1 (MQ=255)

                                                                        tCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTg                                                                               >  1:3084010/1‑150 (MQ=255)

                                                                         ctgTGACCCCGAACACGTCGATCTCTTCTCCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCCGTCACAACGCCGGCGATGATCCGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGa                                                                              <  1:6366672/148‑1 (MQ=25)

                                                                          aGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGAt                                                                             <  1:1429672/150‑1 (MQ=255)

                                                                               ccACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGc                                                                        >  2:1009105/1‑150 (MQ=255)

                                                                               ccACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGc                                                                        >  1:9128141/1‑150 (MQ=255)

                                                                                cACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGTCCTTCACTGCCTTGACGAAGTCTTCGTTGTCACAGTCGTTGATCTGGcc                                                                       >  2:6540841/1‑150 (MQ=255)

                                                                                 aCGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCa                                                                      >  1:1516623/1‑150 (MQ=255)

                                                                                  cGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAg                                                                     >  1:5922846/1‑150 (MQ=255)

                                                                                    aaCACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAggg                                                                   <  1:1570470/150‑1 (MQ=255)

                                                                                     ccctggtgAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAgggc                                                                  <  1:6979376/143‑1 (MQ=255)

                                                                                     acacGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAgggc                                                                  >  2:370380/1‑150 (MQ=255)

                                                                                     acacGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAgggc                                                                  >  1:5746350/1‑150 (MQ=255)

                                                                                       acGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAgggcgg                                                                >  1:8022806/1‑150 (MQ=255)

                                                                                        cGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAgggcggg                                                               <  1:8789459/150‑1 (MQ=255)

                                                                                          tCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTc                                                                                               >  1:1134809/1‑116 (MQ=255)

                                                                                            gAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAAt                                                           <  2:2489128/150‑1 (MQ=255)

                                                                                             aaCTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATg                                                          >  2:6238858/1‑150 (MQ=255)

                                                                                              aCTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGt                                                         <  2:2069347/150‑1 (MQ=255)

                                                                                                tcttctTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAg                                                                                                    >  2:5842641/1‑105 (MQ=255)

                                                                                                tcttctTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGt                                                                                                   <  1:5492511/106‑1 (MQ=255)

                                                                                                tcttctTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTAc                                                       >  1:647151/1‑150 (MQ=255)

                                                                                                acttcttCCATGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTAc                                                       <  2:9145419/149‑1 (MQ=255)

                                                                                                 cttcttCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACg                                                      >  1:1741333/1‑150 (MQ=255)

                                                                                                  ttcttcCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACgg                                                     >  1:8245136/1‑150 (MQ=255)

                                                                                                  ttcttcCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACgg                                                     <  2:3775094/150‑1 (MQ=255)

                                                                                                  ttcttcCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACgg                                                     >  2:3241865/1‑150 (MQ=255)

                                                                                                   tcttcCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGGTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACggg                                                    <  2:1182464/150‑1 (MQ=255)

                                                                                                      tcCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCg                                                 >  2:5398902/1‑150 (MQ=255)

                                                                                                        cAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTc                                               >  2:4589104/1‑150 (MQ=255)

                                                                                                         aGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTcc                                              >  2:1403357/1‑150 (MQ=255)

                                                                                                               gACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAAc                                        >  1:4634116/1‑150 (MQ=255)

                                                                                                                tcAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCATCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGTGGCGTCCTGGAACa                                       <  2:8878854/149‑1 (MQ=255)

                                                                                                                   gggCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGTTCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAggg                                    <  1:7702963/150‑1 (MQ=255)

                                                                                                                    ggCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGc                                   <  1:2023206/150‑1 (MQ=255)

                                                                                                                    ggCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGc                                   <  1:398151/150‑1 (MQ=255)

                                                                                                                     gCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCt                                  >  1:4783930/1‑150 (MQ=255)

                                                                                                                       cGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCttt                                <  2:3064948/150‑1 (MQ=255)

                                                                                                                         ggAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCa                              <  2:4428099/150‑1 (MQ=255)

                                                                                                                          gAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAg                             >  1:6106850/1‑150 (MQ=255)

                                                                                                                              gCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCa                                                                                                                  <  2:8804389/61‑1 (MQ=255)

                                                                                                                              gCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCt                         <  1:1218996/150‑1 (MQ=255)

                                                                                                                               cTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTg                        >  1:8163328/1‑150 (MQ=255)

                                                                                                                               cTACGCACACTTCAGTCACAACGCCGGCGATGATCAGGTGCTTCGTGCCAGTGGCCTGCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTg                        <  2:7022847/150‑1 (MQ=255)

                                                                                                                                tACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTgg                       >  1:5787113/1‑150 (MQ=255)

                                                                                                                                tACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTgg                       <  1:4485728/150‑1 (MQ=255)

                                                                                                                                tACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTgg                       <  2:7320625/150‑1 (MQ=255)

                                                                                                                                tACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTgg                       >  1:4335519/1‑150 (MQ=255)

                                                                                                                                   gCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCAc                    >  2:1621203/1‑150 (MQ=255)

                                                                                                                                     acacTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCa                  <  2:1823800/150‑1 (MQ=255)

                                                                                                                                     acacTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCa                  >  2:4358127/1‑150 (MQ=255)

                                                                                                                                         ttCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGaa                                                                                                     >  1:3435190/1‑63 (MQ=255)

                                                                                                                                         ttCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTTATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTgt              <  2:6004836/150‑2 (MQ=255)

                                                                                                                                         ttCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTgg              >  2:5542100/1‑150 (MQ=255)

                                                                                                                                         ttCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTgg              <  1:4936672/150‑1 (MQ=255)

                                                                                                                                         ttCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTgg              <  1:4102445/150‑1 (MQ=255)

                                                                                                                                               cacaACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGtt        <  2:5092448/150‑1 (MQ=255)

                                                                                                                                                acaACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTg       >  2:1121319/1‑150 (MQ=255)

                                                                                                                                                acaACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTg       <  2:2739549/150‑1 (MQ=255)

                                                                                                                                                 caACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTgg      <  1:7914868/150‑1 (MQ=255)

                                                                                                                                                 caACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTgg      <  2:4983999/150‑1 (MQ=255)

                                                                                                                                                   aCGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTgggg    >  1:868838/1‑150 (MQ=255)

                                                                                                                                                     gCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTGGGGcc  >  1:7536948/1‑150 (MQ=255)

AGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTGGGGCC  >  NC_002947/6000636‑6000934

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