BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,100,381	G→A	5.3%	Y327Y (TAC→TAT)	acsA‑I ←	acetyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,100,381	0	G	A	5.3%	958.4 / 26.1	299	Y327Y (TAC→TAT)	acsA‑I	acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand): ref base G (129/153); new base A (4/12); total (133/166)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 304 total aligned reads displayed.

GGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTTCATCAGGTCTTCGTACCAGAT  >  NC_002947/5100232‑5100530

ggCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCgcggcg                                                                                                                                                       <  2:6160611/150‑1 (MQ=255)

 gCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTCCATGGGTCAGCGCGGCGt                                                                                                                                                      >  1:5152939/1‑150 (MQ=255)

       cTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCa                                                                                                                                                <  2:5741042/150‑1 (MQ=255)

       cTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGCCAGCGCGGCGTAAACCa                                                                                                                                                <  1:743268/150‑1 (MQ=255)

        ttCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAg                                                                                                                                               >  1:1042766/1‑150 (MQ=255)

        ttCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAg                                                                                                                                               >  2:1558239/1‑150 (MQ=255)

        ttCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAg                                                                                                                                               <  1:3922727/150‑1 (MQ=255)

        ttCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTCCATGGGTCAGCGCGGCGTAAACCAg                                                                                                                                               <  2:5152939/150‑1 (MQ=255)

           gAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTa                                                                                                                                            >  1:6186301/1‑150 (MQ=255)

                 cAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGcc                                                                                                                                      >  2:7450538/1‑150 (MQ=255)

                 cAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGcc                                                                                                                                      >  1:3138531/1‑150 (MQ=255)

                    ggtggtGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCtgt                                                                                                                                   <  2:3502833/150‑1 (MQ=255)

                    ggtggtGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCtgt                                                                                                                                   <  1:5621812/150‑1 (MQ=255)

                     gtggtgGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCtgtg                                                                                                                                  >  1:1544555/1‑150 (MQ=255)

                     gtggtgGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCtgtg                                                                                                                                  <  2:1922563/150‑1 (MQ=255)

                      tggtggCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTgtgg                                                                                                                                 >  2:3805846/1‑150 (MQ=255)

                       ggtggCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTgtggt                                                                                                                                <  2:1664123/150‑1 (MQ=255)

                          ggCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGtgtg                                                                                                                             <  2:8147832/150‑1 (MQ=255)

                          ggCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGtgtg                                                                                                                             <  1:3420960/150‑1 (MQ=255)

                           gcgcCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTgc                                                                                                                            >  1:1986897/1‑150 (MQ=255)

                            cgcCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTgca                                                                                                                           >  1:5187950/1‑150 (MQ=255)

                              ccGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTgcagc                                                                                                                         >  1:5315746/1‑150 (MQ=255)

                               cGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGAACCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAAACAGGTAACCGCCTGTGGTGTgcagca                                                                                                                        >  2:1134101/1‑150 (MQ=255)

                                  tGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGc                                                                                                                     >  1:7563552/1‑150 (MQ=255)

                                       aGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTc                                                                                                                >  2:622391/1‑150 (MQ=255)

                                       aGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTc                                                                                                                >  2:7558313/1‑150 (MQ=255)

                                        gCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCg                                                                                                               >  1:8121633/1‑150 (MQ=255)

                                             ccGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTg                                                                                                          >  1:911357/1‑150 (MQ=255)

                                              cGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGc                                                                                                         <  2:8276774/150‑1 (MQ=255)

                                                   aCGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTg                                                                                                    <  1:6379960/150‑1 (MQ=255)

                                                     gATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGa                                                                                                  <  2:2004991/150‑1 (MQ=255)

                                                            cTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAg                                                                                           >  2:4463089/1‑150 (MQ=255)

                                                                ggCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGtat                                                                                       <  1:1812768/150‑1 (MQ=255)

                                                                  ccTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATaa                                                                                     >  1:7335204/1‑150 (MQ=255)

                                                                   cTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATaaa                                                                                    <  2:1330735/150‑1 (MQ=255)

                                                                   cTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATaaa                                                                                    >  2:8941145/1‑150 (MQ=255)

                                                                   cTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATaaa                                                                                    >  2:7285590/1‑150 (MQ=255)

                                                                    tgtgACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAg                                                                                   <  2:5808946/150‑1 (MQ=255)

                                                                     gtgACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAgg                                                                                  <  2:1777158/150‑1 (MQ=255)

                                                                     gtgACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAgg                                                                                  >  1:9074602/1‑150 (MQ=255)

                                                                      tgACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGa                                                                                 <  1:8264979/150‑1 (MQ=255)

                                                                       caCCCAGCCCACATCCGCGGTACACCAGTAAACATCGCCGTCGTGATAGTCGAAGACATATTGATGCGTCATCGAAGCATAGACGAGATAGCCGCCGGTGGTGTGAAGCACGCCCTTTGGCTTGCCGGTCGAGCCGGAGGTATAAAGGAt                                                                                <  1:228677/149‑1 (MQ=255)

                                                                            aaCCGACGTCGGCGGTGCTCCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACa                                                                           <  1:4080134/150‑1 (MQ=255)

                                                                               cGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAgcg                                                                        >  1:3497943/1‑150 (MQ=255)

                                                                                gACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAgcgc                                                                       <  1:7051626/150‑1 (MQ=255)

                                                                                    tCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCctct                                                                   >  2:4699044/1‑150 (MQ=255)

                                                                                     cggcggTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCtt                                                                  <  2:8350752/150‑1 (MQ=255)

                                                                                       gcggTACACCAGTAAACATCGCCGTCGTGATAGTCGAAGACATATTGATGCGTCATCGAAGCATAGACGAGATAGCCGCCGGTGGTGTGAAGCACGCCCTTTGGCTTGCCGGTCGAGCCGGAGGTATAAAGGATGAACAgcgg                                                                       >  1:4786666/1‑142 (MQ=255)

                                                                                       gcggTACACCAGTAAACATCGCCGTCGTGATAGTCGAAGACATATTGATGCGTCATCGAAGCATAGACGAGATAGCCGCCGGTGGTGTGAAGCACGCCCTTTGGCTTGCCGGTCGAGCCGGAGGTATAAAGGATGAACAGCGGGTCTTcc                                                                <  1:8035112/150‑2 (MQ=255)

                                                                                       gcggTACACCAGTAAACATCGCCGTCGTGATAGTCGAAGACATATTGATGCGTCATCGAAGCATAGACGAGATAGCCGCCGGTGGTGTGAAGCACGCCCTTTGGCTTGCCGGTCGAGCCGGAGGTATAAAGGATGAACAGCGGGTCTTcc                                                                >  2:1019407/1‑149 (MQ=255)

                                                                                       gcggTACACCAGTAAACATCGCCGTCGTGATAGTCGAAGACATATTGATGCGTCATCGAAGCATAGACGAGATAGCCGCCGGTGGTGTGAAGCACGCCCTTTGGCTTGCCGGTCGAGCCGGAGGTATAAAGGATGAACAGCGGGTCTTcc                                                                <  2:3918133/150‑2 (MQ=255)

                                                                                       gcggTACACCAGTAAACATCGCCGTCGTGATAGTCGAAGACATATTGATGCGTCATCGAAGCATAGACGAGATAGCCGCCGGTGGTGTGAAGCACGCCCTTTGGCTTGCCGGTCGAGCCGGAGGTATAAAGGATGAACAGCGGGTCTTcc                                                                <  2:3289115/150‑2 (MQ=255)

                                                                                        cggTTCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCTGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAg                                                               <  2:4562476/150‑1 (MQ=255)

                                                                                              aCCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCa                                                         >  1:2026679/1‑150 (MQ=255)

                                                                                               ccAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCAt                                                        <  1:6052555/150‑1 (MQ=255)

                                                                                                cAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCAtt                                                       >  1:2767946/1‑150 (MQ=255)

                                                                                                     cacCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCttt                                                  <  2:4818050/150‑1 (MQ=255)

                                                                                                     cacCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCttt                                                  <  2:4122289/150‑1 (MQ=255)

                                                                                                     cacCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCttt                                                  <  1:7627014/150‑1 (MQ=255)

                                                                                                     cacCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCttt                                                  <  1:6769843/150‑1 (MQ=255)

                                                                                                      acCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTg                                                                                                    >  2:772680/1‑99 (MQ=255)

                                                                                                      acCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTg                                                 <  1:3606273/150‑1 (MQ=255)

                                                                                                       ccTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTgg                                                <  1:8367346/150‑1 (MQ=255)

                                                                                                        cTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTggg                                               <  1:881595/150‑1 (MQ=255)

                                                                                                             ccGGACGGTTATCGAACACCCTTTCATGGGTCAGCGCGGCGTATACCTGGTATCCGCCTGTGGTGTGCAGCACGGCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGACCAGCGCCTCTTCAGCACCCATTTCTTTGGGGgcgc                                          <  2:6798107/150‑1 (MQ=255)

                                                                                                              cGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACACATTTCTTTGGGGGCGCa                                         >  2:5656669/1‑150 (MQ=255)

                                                                                                               ggACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAg                                        >  1:2378322/1‑150 (MQ=255)

                                                                                                                  cGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTgg                                     <  2:7999666/150‑1 (MQ=255)

                                                                                                                   ggTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCAAGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGc                                    >  2:6055105/1‑150 (MQ=255)

                                                                                                                     tAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTg                                  <  1:7746919/150‑1 (MQ=255)

                                                                                                                      aaTCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTgg                                 >  2:8591073/1‑150 (MQ=255)

                                                                                                                      aaTCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTgg                                 <  1:2403157/150‑1 (MQ=255)

                                                                                                                      aaTCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTgg                                 >  2:7198723/1‑150 (MQ=255)

                                                                                                                       aTCTAACACCCGTTCATTGTTCATCGGGGCTTAAACAATGTAAGCTCCTGTGGTTTGCATCTCTCCCTTCGGCTTTCCGGTGGTGAGGGAGGTATAATGGTTGATCAGCGCCTCTCGATCACCCATTTCTTTGGGGGCGCAGTGGCTgta                                <  1:4677743/150‑3 (MQ=255)

                                                                                                                       aTCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGTGGGCGCAGTGGCTGGa                                <  1:1057766/150‑1 (MQ=255)

                                                                                                                       aTCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGa                                >  1:295769/1‑150 (MQ=255)

                                                                                                                        tCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAg                               >  1:4326748/1‑150 (MQ=255)

                                                                                                                        tCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAg                               >  2:6635625/1‑150 (MQ=255)

                                                                                                                        tCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCAGTTCTTTGGGGGCGCAGTGGCTGGAg                               >  2:1201265/1‑150 (MQ=255)

                                                                                                                          gAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGc                             <  1:2707828/150‑1 (MQ=255)

                                                                                                                          gAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGc                             <  1:7083280/150‑1 (MQ=255)

                                                                                                                           aaCACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGcc                            >  2:4156167/1‑150 (MQ=255)

                                                                                                                           aaCACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGcc                            >  2:4046921/1‑150 (MQ=255)

                                                                                                                            acacCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCa                           <  2:2853755/150‑1 (MQ=255)

                                                                                                                            acacCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCa                           >  2:590061/1‑150 (MQ=255)

                                                                                                                                ccGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACttt                       >  2:5239936/1‑150 (MQ=255)

                                                                                                                                ccGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACttt                       >  1:1314730/1‑150 (MQ=255)

                                                                                                                                ccGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACttt                       >  2:6892979/1‑150 (MQ=255)

                                                                                                                                   ttCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCAttt                                                      >  2:2320419/1‑116 (MQ=255)

                                                                                                                                   gccgtcagTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTCCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTtcat                    <  2:7223268/143‑1 (MQ=255)

                                                                                                                                     cATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTtcatca                  <  1:1625607/150‑1 (MQ=255)

                                                                                                                                      aTGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTTCATCAg                 >  2:2153417/1‑150 (MQ=255)

                                                                                                                                      aTGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTTCATCAg                 >  2:5373058/1‑150 (MQ=255)

                                                                                                                                           tCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTTCATCAGGTCtt            >  1:4609814/1‑150 (MQ=255)

                                                                                                                                            cAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTTCATCAGGTCTTc           <  2:3092694/150‑1 (MQ=255)

                                                                                                                                              gcgcGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTTCATCAGGTCTTCGt         <  2:3364445/150‑1 (MQ=255)

                                                                                                                                               cgcgGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTTCATCAGGTCTTCGTa        >  1:658337/1‑150 (MQ=255)

                                                                                                                                               cgcgCCGTAAACCAGGAAACCGCCTGTGGTGTGCAGCACGCCCTTCCGCATGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTTCATCAGGTCTTCGTa        >  1:1389150/1‑150 (MQ=255)

                                                                                                                                                  ggcgTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTTCATCAGGTCTTCGTACCa     <  1:8509426/150‑1 (MQ=255)

                                                                                                                                                     gTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTTCATCAGGTCTTCGTACCAGag  >  2:7688466/1‑149 (MQ=255)

GGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTTCATCAGGTCTTCGTACCAGAT  >  NC_002947/5100232‑5100530

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