BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,100,365	C→G	5.2%	E333Q (GAA→CAA)	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,365	0	C	G	5.2%	1043.6 / 33.0	307	E333Q (GAA→CAA)	acsA‑I	acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand): ref base C (132/159); new base G (4/12); total (136/171)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 308 total aligned reads displayed.

ATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTTCATCA  >  NC_002947/5100217‑5100514

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

  gTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATg                                                                                                                                                    >  2:2403157/1‑150 (MQ=255)

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

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

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

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

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

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

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

                     ccTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATgg                                                                                                                                                   <  1:246854/132‑1 (MQ=255)

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

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

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

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

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

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

                        tCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAgg                                                                                                                              <  1:8859036/150‑1 (MQ=255)

                        tCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTGCATGGGTCAGCGCGGCGTAAACCAgg                                                                                                                              <  2:1730557/150‑1 (MQ=255)

                              aGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCg                                                                                                                        >  2:6786599/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                        cGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCgg                                                                                              <  1:926845/150‑1 (MQ=255)

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

                                                           gCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCtt                                                                                           <  2:1194499/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                                                                   tgtgACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAg                                                                   >  1:6256425/1‑150 (MQ=255)

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

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

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

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

                                                                                      gACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGAt                                                                <  2:1251601/150‑1 (MQ=255)

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

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

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

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

                                                                                                     cgcggTACACCAGTAAACATCGCCGTCGTGATAGTCGAAGACATATTGATGCGTCATCGAAGCATAGACGAGATAGCCGCCGGTGGTGTGAAGCACGCCCTTTGGCTTGCCGGTCGAGCCGGAGGTATAAAGGATGAACAGCGGGTCTTc                                                 <  1:7059864/149‑1 (MQ=2)

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

                                                                                                      gcggTACACCAGTAAACATCGCCGTCGTGATAGTCGAAGACATATTGATGCGTCATCGAAGCATAGACGAGATAGCCGCCGGTGGTGTGAAGCACGCCCTTTGGCTTGCCGGTCGAGCCGGAGGTATAAAGGATGAACAgcgg                                                       <  2:4786666/143‑2 (MQ=255)

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

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

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

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

                                                                                                             aCCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCa                                         <  1:7824849/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                                                             cGGACGGTAAACGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCa                         >  1:9100728/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                                                                                 gTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTtca     <  1:5479504/150‑1 (MQ=255)

                                                                                                                                                 gTTCATGGGTCAGCGCGGCGAAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTtca     <  2:4133907/150‑1 (MQ=255)

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

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

                                                                                                                                                  ttCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTtcat    <  1:1462505/150‑1 (MQ=255)

                                                                                                                                                  ttCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTtcat    >  1:7144011/1‑150 (MQ=255)

                                                                                                                                                   tCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTtcatc   <  1:178328/150‑1 (MQ=255)

ATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGGCCTGTGACCCAACCGACGTCGGCGGTGCACCAGTACACCTCGCCCGGACGGTAATCGAACACCCGTTCATGGGTCAGCGCGGCGTAAACCAGGTAACCGCCTGTGGTGTGCAGCACGCCCTTCGGCTTGCCGGTGGAGCCGGAGGTATAAAGGATGAACAGCGCCTCTTCAGCACCCATTTCTTTGGGGGCGCAGTGGCTGGAGGCCACTTTCATCA  >  NC_002947/5100217‑5100514

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