BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,100,362	G→A	5.4%	R334W (CGG→TGG) ‡	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,362	0	G	A	5.4%	956.1 / 25.6	298	R334W (CGG→TGG) ‡	acsA‑I	acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand): ref base G (131/150); new base A (4/12); total (135/163)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 303 total aligned reads displayed.

GGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTTCA  >  NC_002947/5100213‑5100511

ggTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCg                                                                                                                                                       >  1:5808946/1‑150 (MQ=255)

     tGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCAt                                                                                                                                                  >  2:76213/1‑150 (MQ=255)

     tGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCAt                                                                                                                                                  <  1:5394136/150‑1 (MQ=255)

      gTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCACGGACGGTAATCGAACACCCGTTCATg                                                                                                                                                 >  2:1001339/1‑150 (MQ=255)

       tCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATgg                                                                                                                                                <  2:2752759/150‑1 (MQ=255)

       tCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATgg                                                                                                                                                <  2:1723954/150‑1 (MQ=255)

       tCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATgg                                                                                                                                                >  1:7999666/1‑150 (MQ=255)

        ccGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATggg                                                                                                                                               <  1:70856/150‑1 (MQ=255)

         cGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGt                                                                                                                                              <  1:6949352/150‑1 (MQ=255)

           ggTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCa                                                                                                                                            <  2:9072970/150‑1 (MQ=255)

                ttGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGcgcg                                                                                                                                       >  2:6513570/1‑150 (MQ=255)

                  gggCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCgcggc                                                                                                                                     >  1:7485870/1‑150 (MQ=255)

                    gCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGt                                                                                                                                   >  2:4140409/1‑150 (MQ=255)

                          cTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCa                                                                                                                             >  2:1625607/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

                                        gtggtgGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCACCCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCggtg                                                                                                               >  2:8509426/1‑150 (MQ=255)

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

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

                                               cgcCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTgca                                                                                                        <  2:4636303/150‑1 (MQ=255)

                                                gcCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTgcag                                                                                                       >  2:6961068/1‑150 (MQ=255)

                                                gcCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTgcag                                                                                                       >  2:8107541/1‑150 (MQ=255)

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

                                                  cGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTgcagca                                                                                                     >  2:5680342/1‑150 (MQ=255)

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

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

                                                          aGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTc                                                                                             <  2:2397267/150‑1 (MQ=255)

                                                             gggCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGc                                                                                          <  2:2231354/150‑1 (MQ=255)

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

                                                                     gACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGt                                                                                  <  1:7374498/150‑1 (MQ=255)

                                                                     gACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGt                                                                                  <  2:618431/150‑1 (MQ=255)

                                                                     gACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGt                                                                                  >  1:5167475/1‑150 (MQ=255)

                                                                      aCGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTg                                                                                 >  1:8419889/1‑150 (MQ=255)

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

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

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

                                                                               cTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCGCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTTAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCGGGAg                                                                        <  2:3974868/150‑1 (MQ=255)

                                                                               cTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAg                                                                        <  2:8155786/150‑1 (MQ=255)

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

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

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

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

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

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

                                                                                         tcaCCCAGCCCACATCCGCGGTACACCAGTAAACATCGCCGTCGTGATAGTCGAAGACATATTGATGCGTCATCGAAGCATAGACGAGATAGCCGCCGGTGGTGTGAAGCACGCCCTTTGGCTTGCCGGTCGAGCCGGAGGTATAAAGGa                                                              <  2:8815713/148‑1 (MQ=255)

                                                                                         tcaCCCAGCCCACATCCGCGGTACACCAGTAAACATCGCCGTCGTGATAGTCGAAGACATATTGATGCGTCATCGAAGCATAGACGAGATAGCCGCCGGTGGTGTGAAGCACGCCCTTTGGCTTGCCGGTCGAGCCGGAGGTATAAAGGa                                                              <  1:4202281/148‑1 (MQ=255)

                                                                                         tcaCCCAGCCCACATCCGCGGTACACCAGTAAACATCGCCGTCGTGATAGTCGAAGACATATTGATGCGTCATCGAAGCATAGACGAGATAGCCGCCGGTGGTGTGAAGCACGCCCTTTGGCTTGCCGGTCGAGCCGGAGGTATAAAGGa                                                              <  2:5859896/148‑1 (MQ=255)

                                                                                             ccAGCCCACATCCGCGGTACACCAGTAAACATCGCCGTCGTGATAGTCGAAGACATATTGATGCGTCATCGAAGCATAGACGAGATAGCCGCCGGTGGTGTGAAGCACGCCCTTTGGCTTGCCGGTCGAGCCGGAGGTATAAAGGATGaa                                                          >  2:668923/1‑150 (MQ=255)

                                                                                             ccAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGaa                                                          >  2:578933/1‑150 (MQ=255)

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

                                                                                               aaCCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACa                                                        >  2:519624/1‑150 (MQ=255)

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

                                                                                                   gACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAgcgc                                                    >  2:5270854/1‑150 (MQ=255)

                                                                                                     cGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCct                                                  <  2:7410584/150‑1 (MQ=255)

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

                                                                                                         ggcggTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTc                                              <  1:8222855/150‑1 (MQ=255)

                                                                                                         ggcggTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTc                                              <  2:7960860/150‑1 (MQ=255)

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

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

                                                                                                           cggTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAg                                            <  2:4562461/150‑1 (MQ=255)

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

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

                                                                                                                  ccAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCAt                                     >  2:2941179/1‑150 (MQ=255)

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

                                                                                                                      tACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCt                                 >  1:8054702/1‑150 (MQ=255)

                                                                                                                      tACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCt                                 >  1:3494308/1‑150 (MQ=255)

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

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

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

                                                                                                                         acCTCGCCCGGACGGTCATCGAACACCCGTTCCTGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCGCTTCAGCACCCCTTTCtggg                              >  2:5908434/1‑147 (MQ=255)

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

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

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

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

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

                                                                                                                            tCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAg                                                                        >  1:2232232/1‑105 (MQ=255)

                                                                                                                                ccGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCttt                               >  1:1725366/1‑142 (MQ=255)

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

                                                                                                                                      ggTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATGAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGc                 <  1:3551984/150‑1 (MQ=255)

                                                                                                                                       gTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCt                <  2:2448974/150‑1 (MQ=255)

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

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

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

                                                                                                                                         aaTCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTgg              >  1:2475093/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)

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

                                                                                                                                            cGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTgc                                                                                                         <  1:8284413/56‑1 (MQ=255)

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

                                                                                                                                                 acCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACt      >  2:4276987/1‑150 (MQ=255)

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

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

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

GGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTTCA  >  NC_002947/5100213‑5100511

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